Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,420,328 G→A 52.9% intergenic (‑61/+160) 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,3280GA52.9% 3.8 / 40.2 36intergenic (‑61/+160)hisJ/argThistidine/lysine/arginine/ornithine transporter subunit/lysine/arginine/ornithine transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base A (8/11);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

TTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACA  >  CP009273/2420188‑2420471
                                                                                                                                            |                                                                                                                                               
ttttGCGGAAACACATCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGATCCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaa                                                                                                                                         <  1:139308‑M1/149‑3 (MQ=255)
             cagcaAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGGTAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATAGTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatc                                                                                                                            >  1:275081‑M1/1‑134 (MQ=255)
               gcaAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaa                                                                                                                          <  1:314357‑M1/149‑18 (MQ=255)
                  aaCGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCCCCAGTTGTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGcgaca                                                                                                                       >  1:79105/1‑145 (MQ=255)
                     gccgcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATCATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGCTTGCTACACTTATCGTAAAATAAaattaattgaccatcga                                                                                                                          >  2:449931‑M1/1‑126 (MQ=255)
                     gccgcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgagcatcaa                                                                                                                          >  2:124444‑M1/1‑126 (MQ=255)
                       cgcAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACCCTTATCGTGCCATAAAATTGCCCATCAAGGCaca                                                                                                                        >  2:390491/1‑141 (MQ=255)
                               cgcTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgcc                                                                                                          >  1:305916‑M1/1‑116 (MQ=255)
                                gcTGGGAAAGGTCAGCACCAGAGAGAGCTATACCACCAGTTTTTGCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattaaccatcaaggcaaaatacc                                                                                                               <  2:411349‑M1/143‑29 (MQ=255)
                                 cTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGcc                                                                                                              >  2:444633/1‑143 (MQ=255)
                                  tGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAg                                                                                                       >  1:406440/1‑149 (MQ=255)
                                   ggAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGa                                                                                                      >  1:413540/1‑149 (MQ=255)
                                   ggAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgt                                                                                                      >  1:525408‑M1/1‑112 (MQ=255)
                                        aGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCTACTTACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGtat                                                                                                 <  1:30790/149‑1 (MQ=255)
                                                aCCAGAGAGAGCGATATCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatt                                                                                         <  1:78480‑M1/149‑51 (MQ=255)
                                                 ccAGAGAGAGCGTTAGCACCAGTTTTTTCTTTATGTCCTCAACTGGCAGTCTTTTTATAACGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattattt                                                                                        <  1:259877‑M1/149‑52 (MQ=255)
                                                 ccAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGtattat                                                                                              >  2:86711/1‑143 (MQ=255)
                                                       agagCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCa                                                                                        <  2:17882/143‑1 (MQ=255)
                                                       agacCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTACCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCa                                                                                        <  2:339781/143‑1 (MQ=255)
                                                               aGCACCATTTTTTTCATAATGCCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTaa                                                                          <  1:397120/149‑1 (MQ=255)
                                                                             cATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTCTGCATTTGAAt                                                            >  1:129427/1‑149 (MQ=255)
                                                                             cATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCACAAAATTGACCATCAATGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATACTTAAGTATGCATTTGAAt                                                            >  1:318997/1‑149 (MQ=255)
                                                                                    tCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAACGCTTAAGTATGCATTTAAATg                                                           >  2:217060/1‑143 (MQ=255)
                                                                                             gACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcgaaatatcccgccccgttagacgtattatttttcattaaacatgcttacgtatgaatttccatgcccaa                                                  >  2:129255‑M1/1‑54 (MQ=255)
                                                                                                        tataGGGAATATTTACAGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatcgtgcaa                                       >  2:454441‑M1/1‑43 (MQ=255)
                                                                                                             ggAATATTTACAGGTTAGTTGCACTTATCGTAACATAAaattaattgaccatcaaggcaaagtatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaac                                  >  2:333787‑M1/1‑38 (MQ=255)
                                                                                                              gAATATTTACAGGTTTGCTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGc                                 <  2:294257/143‑1 (MQ=255)
                                                                                                                   tttACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCa                            >  2:377532/1‑143 (MQ=255)
                                                                                                                       cAGGTTTGTGACCCTTATCGTAACATAAaattaatagaccatcaaggcaaaatctcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatga                  <  1:397017‑M1/149‑122 (MQ=255)
                                                                                                                        aGGTTTGTTACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttctttaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctgccgtgcaccatgat                 <  1:361877‑M1/149‑123 (MQ=255)
                                                                                                                          gTTTCTCACACTTATCGTAACATAAaattaatggcccatcaaggcaaaatatcccgccccgttagacgtattatttttcattacaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcacca                     <  2:426756‑M1/143‑119 (MQ=255)
                                                                                                                           tttGTTACACTTATCAGGCCATAAAAGTGACCATCAATGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCAtt              <  1:302645/149‑1 (MQ=255)
                                                                                                                               ttACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatg                <  2:270891‑M1/143‑124 (MQ=255)
                                                                                                                                tACACTTATCGTAACATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatga               <  2:475695‑M1/143‑125 (MQ=255)
                                                                                                                                     ttATCGTAATATAAaattaattgaccatcaaggcaaaatatcccgccccgttagacgtattatttttcattaaaaatgcttaagtatgcatttcaatgcaccacgatagtgcaacaaacctggcgtgcaccatgatgatcatt          <  2:213782‑M1/143‑130 (MQ=255)
                                                                                                                                       ctCGTGCCATAAATATTACCATCAAGGCAAAATTTCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATCCATTTCAATGCACTACGATAGTGCAACAAACCTGGCGAGCACCATGATGATCATTCCATCAGGTACa  <  1:322059/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
TTTTGCGGAATCGCAGCAAACGCCGCAGTTGCGCTGGAGAAGGCCAGAACCAGAGAGAGCGATAGCACCAGTTTTTTCATAATGTCCTCAACTGACAGTCTTTTTATAGGGAATATTTACAGGTTTGTTACACTTATCGTGCCATAAAATTGACCATCAAGGCAAAATATCCCGCCCCGTTAGACGTATTATTTTTCATTAAAAATGCTTAAGTATGCATTTCAATGCACCACGATAGTGCAACAAACCTGGCGTGCACCATGATGATCATTCCATCAGGTACA  >  CP009273/2420188‑2420471

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: