Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,301,970 +GT coding (789/1632 nt) oppA → oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,301,9681.G92.0% 83.7 / 0.4 25coding (787/1632 nt)oppAoligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (10/13);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
*NC_0009131,301,9682.T92.0% 74.7 / 0.4 25coding (787/1632 nt)oppAoligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (10/13);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

ATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑‑G‑‑TAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTA  >  NC_000913/1301726‑1302203
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                             
ataAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTAGCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCGTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCGATTTGCCTATTGCTGCTGAAGTTACCGATGTCAACCGCTACC‑‑G‑‑tgta                                                                                                                                                                                                                                        <  1:239462/248‑4 (MQ=255)
                           cACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGc                                                                                                                                                                                                             >  1:422991/1‑250 (MQ=255)
                            aCCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACa                                                                                                                                                                                                               >  2:357038/1‑247 (MQ=255)
                            aCCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACa                                                                                                                                                                                                               <  1:357038/247‑1 (MQ=255)
                                       aaGCCGCTATCGAGACATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAAt                                                                                                                                                                                                >  1:99972/1‑251 (MQ=255)
                                       aaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGAATGCCGATCGAAt                                                                                                                                                                                                >  2:101720/1‑251 (MQ=255)
                                       aaGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTGACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAAt                                                                                                                                                                                                <  2:322005/251‑1 (MQ=255)
                                                                                        cGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATcc                                                                                                                                               >  1:329919/1‑251 (MQ=255)
                                                                                            cccAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAATCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATAc                                                                                                                                           <  1:174476/250‑1 (MQ=255)
                                                                                              cAACGGTGCCTATACCTTAAAAGTTTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAGCCAGGTAACCTATTTGCCTATTGCTTCTGAAGGTACCGATGTCAACCGCTACC‑‑GGGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTTTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCt                                                                                                                                         <  1:334883/251‑1 (MQ=255)
                                                                                                         ataCCGTAAAAGATTTGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACt                                                                                                                              <  2:296453/251‑1 (MQ=255)
                                                                                                            ccTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTAtt                                                                                                                           >  2:75429/1‑251 (MQ=255)
                                                                                                                       tGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCGATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTaa                                                                                                                 <  2:82778/250‑1 (MQ=255)
                                                                                                                               aaaCGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTTCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGGGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACAAAGTTCACGTTGATCCATACCTGTGCATTTACTATTTCGAAATGAACAACCAgc                                                                                                         <  2:111203/250‑2 (MQ=255)
                                                                                                                               aaaCGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGa                                                                                                         <  1:75429/250‑1 (MQ=255)
                                                                                                                                   gaacgaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCg                                                                                                    <  1:322047/251‑1 (MQ=255)
                                                                                                                                                                               gAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTAt                                                        >  1:47837/1‑251 (MQ=255)
                                                                                                                                                                                        tattaACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTTTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACcgcg                                                <  2:99972/250‑1 (MQ=255)
                                                                                                                                                                                                 ggTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACTGCGATATCATTGt                                      <  1:20555/251‑1 (MQ=255)
                                                                                                                                                                                                  gTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGCATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCACCGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGtt                                     >  2:79958/1‑251 (MQ=255)
                                                                                                                                                                                                             gcctttttcTTCTGAAGTTTCTGATGTCAACCGCAACCGTG‑‑TAGTGGGGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATCCCTGTGCACTAACTTTTTCGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGTGATATCATTGTTAATAAAGTGaa                          <  2:252845/243‑1 (MQ=255)
                                                                                                                                                                                                                      ttCTGAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAggg                 >  2:144649/1‑251 (MQ=255)
                                                                                                                                                                                                                          gAAGTTACCGATGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGATATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACACCCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCa               <  1:332021/249‑1 (MQ=255)
                                                                                                                                                                                                                                     tGTCAACCGCTACCGTG‑‑TAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCt   >  1:266667/1‑250 (MQ=255)
                                                                                                                                                                                                                                      gTCAACCGCTACCGTG‑‑TAGTGGTGCAATCGACATGAAGTACAACATCATGCTGATGGAATTTTTCCAGAAGCTGAAAAAAGAGACCCCGGACAAAGTTCACGTTGATCCATGCTGGTGCCTTTATTATTACAAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACGAGGTATGGACCGCGATATCTTTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTa  <  2:257821/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                             
ATAAATTACTTGTTCACCCATCAACTTCACCGGTGCCAAAAGCCGCTATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACC‑‑G‑‑TAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTA  >  NC_000913/1301726‑1302203

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