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

CTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGT  >  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)
                                                                           gggAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑Ca                                                       >  2:286340/1‑150 (MQ=255)
                                                                                  tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCatgat                                                <  2:286012/150‑1 (MQ=255)
                                                                                     aCAGTTTTTTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATgatgat                                             <  2:251526/150‑1 (MQ=255)
                                                                                       aGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCa                                           <  1:249771/150‑1 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGAATGCATTTCAATGCACCACGATAGAGCAACAAACCTGACGTG‑CACCataccaatcat                                          >  1:162397/1‑141 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACCAACCTGGCGTG‑CACCATGATGATCAt                                          >  1:131753/1‑150 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATCTTTCATTAAAAATGCTCAAGTATCCATTTCAATGCACCACGATAGAGCAACAAACCTGGCATC‑CACCATAACGATCAt                                          >  1:120663/1‑150 (MQ=255)
                                                                                        ggTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAAGATCCCGCCCCGTTAGACGTATTATTTTTCATGAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCAt                                          <  1:9404/150‑1 (MQ=255)
                                                                                               tACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCa                                   >  2:167720/1‑150 (MQ=255)
                                                                                                acacTTATCGTGCCATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CAAAatgaaaaggaaaacattag                                  <  1:182078/150‑16 (MQ=255)
                                                                                                         gTGCCATAAAATTGACCATCAACGCAAGATATCCTGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACACCCCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCtt                         <  2:62412/150‑1 (MQ=255)
                                                                                                            ccATAAAATTGACCATCAAGGCAAGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTccc                      >  1:136826/1‑150 (MQ=255)
                                                                                                             gaTAAAATTGACCATCAAGGCAAAGTATCCCGCCCCGTTTGACGTATGATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGGTAGTGCAACAAACCTGGCGAG‑CACCATGACAACCATCCCATCAGGTACACCTTCCCa                     >  1:268146/2‑150 (MQ=255)
                                                                                                             cATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTCCCCa                     >  2:106688/1‑150 (MQ=255)
                                                                                                               tAAAATTGACCATCAAGGCAGGATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTCCCAGc                   >  2:260685/1‑150 (MQ=255)
                                                                                                                 aaaTTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGa                 >  2:276086/1‑150 (MQ=255)
                                                                                                                  aaTTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGAc                >  2:138174/1‑150 (MQ=255)
                                                                                                                         cATCAAGGCAAAATTTCCCGCCCCGTTAGACGTATTATTGTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAg         <  2:20292/150‑1 (MQ=255)
                                                                                                                          aTCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGGATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGt        <  2:206688/150‑1 (MQ=255)
                                                                                                                                gCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATATAAATACAACACGATAGTGCAACAAACCTGGCGTGACAACTTCATG‑TCATTCCTTCAGGGCCTGATTCCCAGCGACGTATCAGTCACCGt  >  2:315970/1‑150 (MQ=255)
                                                                                                                                    |                                                                                                                                                  
CTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTG‑CACCATGATGATCATTCCATCAGGTACAGCTTCCCAGCGACGTATCAGTCACCGT  >  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: ‑