Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,431,197 A→G 11.4% S79S (TCA→TCG USA300HOU_2294 → transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,431,1970AG11.4% 87.8 / 3.4 35S79S (TCA→TCGUSA300HOU_2294transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (18/13);  new base G (2/2);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTTTAAGAGATATTGATGATTTAGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGAAAATTGCTTAAACGCATGGACTTTTATAT  >  CP000730/2431065‑2431337
                                                                                                                                    |                                                                                                                                            
aaTTTTAAGAGATATTGATGATTTAGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTAtt                                                                                                                                      >  1:248232/1‑141 (MQ=255)
        gagaTATTGATGATTTAGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGc                                                                                                                              <  2:424055/141‑1 (MQ=255)
        gagaTATTGATGATTTAGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGc                                                                                                                              <  1:182811/141‑1 (MQ=255)
          gaTATTGATGATTTAGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCGGCCTTATTTTTAGTGCtt                                                                                                                            <  2:155320/141‑1 (MQ=255)
           atatTGATGATTTAGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTa                                                                                                                           >  1:263854/1‑141 (MQ=255)
                        aGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAagtcagtc                                                                                                              <  1:358462/141‑1 (MQ=255)
                        aGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAagtcagtc                                                                                                              >  1:438575/1‑141 (MQ=255)
                         gAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCg                                                                                                             >  2:33389/1‑141 (MQ=255)
                               cAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCg                                                                                                       <  2:73535/141‑1 (MQ=255)
                                             ttATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATaa                                                                                         <  1:74922/141‑1 (MQ=255)
                                               aTGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATaaaa                                                                                       >  2:64565/1‑141 (MQ=255)
                                                 gcgcATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCGGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGc                                                                                     >  1:366866/1‑141 (MQ=255)
                                                 gcgcATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGc                                                                                     >  1:442304/1‑141 (MQ=255)
                                                 gcgcATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGc                                                                                     <  2:175035/141‑1 (MQ=255)
                                                       tATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATc                                                                               <  2:304218/141‑1 (MQ=255)
                                                       tATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATc                                                                               <  1:208432/141‑1 (MQ=255)
                                                          gggAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCGGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAAc                                                                            <  2:433563/141‑1 (MQ=255)
                                                            gAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGc                                                                          >  1:294509/1‑141 (MQ=255)
                                                              aaaaTGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCaa                                                                        >  1:192221/1‑141 (MQ=255)
                                                              aaaaTGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCaa                                                                        >  2:36792/1‑141 (MQ=255)
                                                                aaTGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAAtt                                                                      >  2:13733/1‑141 (MQ=255)
                                                                   ggtggtTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTCGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATtata                                                                   >  2:449967/1‑141 (MQ=255)
                                                                   ggtggtTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATtata                                                                   >  1:78287/1‑141 (MQ=255)
                                                                    gtggtTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATaa                                                                  >  1:267880/1‑141 (MQ=255)
                                                                     tggtTACCAAAAACAACAACCACAAACTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATaaa                                                                 <  1:172576/141‑1 (MQ=255)
                                                                         tACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACaa                                                             <  2:209384/141‑1 (MQ=255)
                                                                            cAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGt                                                          >  2:32852/1‑141 (MQ=255)
                                                                                          acaATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCcaca                                            >  1:322406/1‑141 (MQ=255)
                                                                                            aaTCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACaaa                                          <  1:12690/141‑1 (MQ=255)
                                                                                               cTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAacac                                       >  1:172039/1‑141 (MQ=255)
                                                                                                          ttAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGaaaa                            <  1:341317/141‑1 (MQ=255)
                                                                                                          ttAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGaaaa                            >  1:203461/1‑141 (MQ=255)
                                                                                                                       aaaCACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGAAAATTGCTTAAACGCa               <  1:138766/141‑1 (MQ=255)
                                                                                                                           acaATTATCGGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGAAAATTGCTTAAACGCATGGa           >  1:161520/1‑141 (MQ=255)
                                                                                                                               ttATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGAAAATTGCTTAAACGCATGGACttt       >  1:65450/1‑141 (MQ=255)
                                                                                                                                    aGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGAAAATTGCTTAAACGCATGGACTTTTatat  <  2:334632/141‑1 (MQ=255)
                                                                                                                                    aGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGAAAATTGCTTAAACGCATGGACTTTTatat  <  2:256420/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                            
AATTTTAAGAGATATTGATGATTTAGAAAATCAAGGTGTTAAAATTTATGCGCATTATGGGAAAAATGGTGGTTACCAAATACAACAAGCACAATCTAAAATTGCATTAAACTTATCTGAAACACAATTATCAGCCTTATTTTTAGTGCTTAATGAAAGTCAGTCGTACTCGACATTACCATATAAAAGCGAAATCAACGCAATTATAAAACAATGTTTAAGTCTTCCACAAACACGCTTAAGAAAATTGCTTAAACGCATGGACTTTTATAT  >  CP000730/2431065‑2431337

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: