Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 66,157 T→G 19.2% V447V (GTT→GTG U712_00295 → putative membrane protein yabM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00688166,1570TG19.2% 92.7 / 4.3 48V447V (GTT→GTGU712_00295putative membrane protein yabM
Reads supporting (aligned to +/- strand):  ref base T (9/29);  new base G (8/2);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 1.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTAAACACACTTTTGGTTCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGG  >  CP006881/66021‑66294
                                                                                                                                        |                                                                                                                                          
tgtTAAACACACTTTTGGTTCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGtttt                                                                                                                                          >  3:379875/1‑139 (MQ=255)
   tAAACACACTTTTGGTTCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGAt                                                                                                                                       >  2:125059/1‑139 (MQ=255)
   tAAACACACTTTTGGTTCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGAt                                                                                                                                       >  2:118463/1‑139 (MQ=255)
          acTTTTGGTTCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGatt                                                                                                                                >  2:356441/1‑139 (MQ=255)
              ttGGTTCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGattat                                                                                                                             <  6:133129/137‑1 (MQ=255)
                  ttCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTTTTTGATc                                                                                                                       >  4:340290/1‑139 (MQ=255)
                     cgctCTATGGGATTGTAGGAGCATCTCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCTGCAAAATGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTg                                                                                                                     <  5:356259/138‑1 (MQ=255)
                     cgcgCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTg                                                                                                                     <  6:406891/138‑1 (MQ=255)
                     cgcgCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTg                                                                                                                     <  1:95993/138‑1 (MQ=255)
                     cgagcTATGGGATTGAAGGCGCATCGCTTGCGTCAGCGGCTTCCTTTGCGGCTGGCGCCTGCCTTTAATTGTTTCATTTTCTGCATTAGGAATGGCTCGATAAGCTTTGAGGCGTGTTGACTCCGATTATTGGATCTg                                                                                                                     <  4:427650/135‑1 (MQ=255)
                                    aaGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGTTTAGGCTTATTTGATCTGCTCTTTTGAtttgagg                                                                                                     >  1:391309/1‑133 (MQ=255)
                                    aaGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGGTTGGGGTTATTGGTTCTGCTCTTTTGATGTGAgg                                                                                                     >  7:388671/1‑138 (MQ=255)
                                    aaGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGATTACGATTATTGGATCTGGTCTTTTGATGTCAGc                                                                                                     >  6:211503/1‑139 (MQ=255)
                                           aTCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTGCGATTTTTTGATATGGTCTTTTGATGTCAGCTGTTttg                                                                                              >  8:394062/1‑139 (MQ=255)
                                            tCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTttgc                                                                                             <  7:38131/139‑1 (MQ=255)
                                            tCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTttgc                                                                                             <  3:145475/139‑1 (MQ=255)
                                            tCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTttgc                                                                                             <  7:65763/139‑1 (MQ=255)
                                                  gCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGCTTGCGATTATTGGGTGTGGTCGTTTGATGTCAGGTGGTTTGGTTGGGt                                                                                       >  1:289668/1‑139 (MQ=255)
                                                       aGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATAcg                                                                                  <  5:364896/139‑1 (MQ=255)
                                                       aGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATAcg                                                                                  <  7:394062/139‑1 (MQ=255)
                                                       aGCGGCGTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTGCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATAcg                                                                                  <  1:149302/139‑1 (MQ=255)
                                                       aGCGGCGTCCTTTGCGGCTGTCGCCGGCCGTAACTTGTATCCGCTGCGGCAACAGGAATGGCTCGCTAAGCTTAGAGGCGTTTTGATTACGATTATTGGATCTGCTCTTTTGATGGCAGCTGTTTTGCGGGCGTATAcg                                                                                  <  4:226457/139‑1 (MQ=255)
                                                              tCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGTTTGCGTATTTTTGTTGTGGTCTTTTGTTGTGAGCGGGTTTGGTTGGGGATAggcgggggt                                                                           >  3:295892/1‑134 (MQ=255)
                                                                   tGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGa                                                                        <  5:211503/137‑1 (MQ=255)
                                                                   tGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGa                                                                        <  6:202068/137‑1 (MQ=255)
                                                                   tGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGa                                                                        <  1:356441/137‑1 (MQ=255)
                                                                   tGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGa                                                                        <  8:90239/137‑1 (MQ=255)
                                                                   tGCGGCGGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGa                                                                        <  6:337320/137‑1 (MQ=255)
                                                                   tGCGGCGGGCGCCGGCCTTTACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGGTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGGCTCTGGa                                                                        <  7:111503/137‑1 (MQ=255)
                                                                          gTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGTTTGCGATTATTGGTTCTGGGCGTTTGATGTGAGCTGTTTTTCTTTGGTATTAGCGTCTCTGGGCTTtgttgt                                                               >  7:298100/1‑139 (MQ=255)
                                                                               cGGCCTTAACTTGTATCAGCGTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTGCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTGGCGTATACGCGTCTCTGGACTTTCTTGTATCCt                                                          <  3:103578/139‑1 (MQ=255)
                                                                                  ccTTAACTTTTATCAGCTTCGTCAATAGGAATGGCTCGATATGATTAGAGGCGTGTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTTTTTTGCTTGCTTTTACGCGTCTCTGGACTTTCTTGTTTGCTGCa                                                       <  2:8832/139‑1 (MQ=255)
                                                                                  ccTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCa                                                       <  7:407844/139‑1 (MQ=255)
                                                                                           gTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGATGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAag                                               <  3:37509/138‑1 (MQ=255)
                                                                                            tATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTTATTACGATTATTTTTTGTGGGATTTTTTTTTGATGTGGTTTTGTTTCGTATACCCGTCTATGTACTTTTATTTTTTATTCAACAGGAagag                                             >  4:64598/1‑139 (MQ=255)
                                                                                            tATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGTTGTTTGCGATTTTTTGGTATGCTGTTTTGATGTGAGCTGTTTTGCTTGGGGATACGCGTGTCTGGGCGTTGTTGTTTCCTGCGAAAGGAagag                                             >  5:128029/1‑139 (MQ=255)
                                                                                            tATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGGTGATTAAGATTATTTGGATTGCTCTTTTGATGTGAGGAGTTTTTCTTGCGTATTCTCTGGGCTGGTGTGTCTTGTTTCATGCAACGGGAagag                                             >  8:240600/1‑139 (MQ=255)
                                                                                                  cTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTGGTTTTTACGATTATTTGATCGGGTGTTTTTTTGTAAGCGGGTTTGGTTTGGTTTACTCGGCGCGGGTCTTTGTTTTTTGCTGGAACAGGAAGAGGcgcgc                                       >  6:279626/1‑138 (MQ=255)
                                                                                                      ggCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTTGAt                                                                                                                        >  6:178556/1‑54 (MQ=255)
                                                                                                      ggCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGAt                                                                                                                        <  5:178556/54‑1 (MQ=255)
                                                                                                               aaTGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGagag                          <  3:177293/139‑1 (MQ=255)
                                                                                                                aTGTCTCTATAATCTTAGAGGCGTTTTGATTCCTTTTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGGGTATCAGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGagag                          <  3:64598/138‑1 (MQ=255)
                                                                                                                    cTCGATAAGCTTAGAGGCGTTTTGATTCGGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTTCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCttt                     <  3:325022/139‑1 (MQ=255)
                                                                                                                       gATACGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTctg                  <  1:118463/139‑1 (MQ=255)
                                                                                                                       gATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCGGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGGGATTGAGAGTCTTTctg                  <  1:338105/139‑1 (MQ=255)
                                                                                                                       gATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGGCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTctg                  <  1:125059/139‑1 (MQ=255)
                                                                                                                         ttaGCTTAGAGGCGTTTTGATTCCGATTATTGGATATGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTctgct                <  7:240600/137‑1 (MQ=255)
                                                                                                                         tAAGCTTAGAGGCGTTTTGATTGGGATTTTTGGAGGTGCTGTTTTGTTTTAAGCGGTTTTTGTTGGGTATACGCGGGGGTGGACTTTCTTGTTTGCTGGGACAGGGAGTGGGGGGGCGGGGGTTGAGAGGCTTTctggg                >  6:105913/1‑137 (MQ=255)
                                                                                                                         tAAGATTAGAGGCGTTTTGATTCCGATTATTGTATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGAATTTCTTGTTTCCTGCAACAGCAAGAGGAGCGGCTGTGATTGAGAGTCTTTctgct                <  1:14828/139‑1 (MQ=255)
                                                                                                                                       ttttGATTCCGATTATTTGTTCTGCTGTTTTTTTTTGAGCGGATTTGCTTGCGTATAAGCGTCTCTGGGCATTATTGTTTACTGCAACACGAATAGGAGCGGCGGTGATTGAGAGTCTTTCTGCTGTTGCGATtggtgg  >  8:105571/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TGTTAAACACACTTTTGGTTCCGCGCTATGGGATTGAAGGCGCATCGCTTGCGACAGCGGCTTCCTTTGCGGCTGTCGCCGGCCTTAACTTGTATCAGCTTCGGCAAAAGGAATGGCTCGATAAGCTTAGAGGCGTTTTGATTCCGATTATTGGATCTGCTCTTTTGATGTCAGCTGTTTTGCTTGCGTATACGCGTCTCTGGACTTTCTTGTTTCCTGCAACAGGAAGAGGAGCGGCTGTGATTGAGAGTCTTTCTGCTGTTGCCATTGGTGG  >  CP006881/66021‑66294

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: