Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,353 A→G 44.8% intergenic (‑86/+135) hisJ ← / ← argT histidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,420,3530AG44.8% 3.3 / 30.4 29intergenic (‑86/+135)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (10/6);  new base G (4/9);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

CTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGT  >  CP009273/2420221‑2420498
                                                                                                                                    |                                                                                                                                                  
cTGGAGCAGGCCAGCACCAGGTAGAGTGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTGTTATAGGGAATATTTACAGGTTTGTTACACGTCTCGTTCCATAACATTGCCCATCAAGGCAAGATATCCCGCCCCGTTAg                                                                                                                                   <  1:180646/150‑1 (MQ=255)
cTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTGCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAg                                                                                                                                   <  1:68176/150‑1 (MQ=255)
           cAGAACCAGAGAGAGCGATAGCACCAGTTTTTGCATAATGTCCTCAACTGACAGTTTTTTAAAAACGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTAttt                                                                                                                        <  2:261391/150‑1 (MQ=255)
                                caacAGTTTTTTCATAATGGCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTa                                                                                                   <  2:23934/147‑1 (MQ=255)
                                cACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTa                                                                                                   <  2:23919/150‑1 (MQ=255)
                                  ncAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATg                                                                                                 >  2:222824/2‑150 (MQ=255)
                                             aTAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGc                                                                                      >  2:235986/1‑150 (MQ=255)
                                                             aCAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGcaac                                                                      <  1:273296/150‑1 (MQ=255)
                                                                 tCTTGTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAAc                                                                  <  2:166375/150‑1 (MQ=255)
                                                                           gggAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCa                                                       >  2:286340/1‑150 (MQ=255)
                                                                                  tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCatgat                                                <  2:286012/150‑1 (MQ=255)
                                                                                     aCAGTTTTTTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATgatgat                                             <  2:251526/150‑1 (MQ=255)
                                                                                       aGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCa                                           <  1:249771/150‑1 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGAATGCATTTCAATGCACCACGATAGAGCAACAAACCTGACGTGCACCataccaatcat                                          >  1:162397/1‑141 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACCAACCTGGCGTGCACCATGATGATCAt                                          >  1:131753/1‑150 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATCTTTCATTAAAAATGCTCAAGTATCCATTTCAATGCACCACGATAGAGCAACAAACCTGGCATCCACCATAACGATCAt                                          >  1:120663/1‑150 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAAGATCCCGCCCCGTTAGACGTATTATTTTTCATGAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCAt                                          <  1:9404/150‑1 (MQ=255)
                                                                                               tACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCa                                   >  2:167720/1‑150 (MQ=255)
                                                                                                acacTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCAAAatgaaaaggaaaacattag                                  <  1:182078/150‑16 (MQ=255)
                                                                                                         gTGCCATAAAATTGACCATCAACGCAAGATATCCTGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACACCCCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCtt                         <  2:62412/150‑1 (MQ=255)
                                                                                                            ccATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTccc                      >  1:136826/1‑150 (MQ=255)
                                                                                                             gaTAAAATTGACCATCAAGGCAAAGTATCCCGCCCCGTTTGACGTATGATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGGTAGTGCAACAAACCTGGCGAGCACCATGACAACCATCCCATCAGGTACACCTTCCCa                     >  1:268146/2‑150 (MQ=255)
                                                                                                             cATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTCCCCa                     >  2:106688/1‑150 (MQ=255)
                                                                                                               tAAAATTGACCATCAAGGCAGGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGc                   >  2:260685/1‑150 (MQ=255)
                                                                                                                 aaaTTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGa                 >  2:276086/1‑150 (MQ=255)
                                                                                                                  aaTTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGAc                >  2:138174/1‑150 (MQ=255)
                                                                                                                         cATCAAGGCAAAATTTCCCGCCCCGTTAGACGTATTATTGTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAg         <  2:20292/150‑1 (MQ=255)
                                                                                                                          aTCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGGATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGt        <  2:206688/150‑1 (MQ=255)
                                                                                                                                gCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATATAAATACAACACGATAGTGCAACAAACCTGGCGTGACAACTTCATGTCATTCCTTCAGGGCCTGATTCCCAGCGACGTATCAGTCACCGt  >  2:315970/1‑150 (MQ=255)
                                                                                                                                    |                                                                                                                                                  
CTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGT  >  CP009273/2420221‑2420498

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