Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,302,204 Δ1 bp coding (1023/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,302,2040T.91.7% 90.5 / 3.8 24coding (1023/1632 nt)oppAoligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (11/11);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

AATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACT  >  NC_000913/1301979‑1302434
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
aaTCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATgg                                                                                                                                                                                                              >  1:221345/1‑251 (MQ=255)
         gtcgtaTAACACCAGGCCGACGGAATTGTTCCAGAAGCTGAAAAAAGAGTTCCCGGACGAAGTTCAGTTTGATCTATGCTTGTGCACTTACTATGACGAAATTAACTGCTATATTCCGCCATTCAACGATGTGGGTGTTCGTACCGCTCTGAAACTAGGTATTGACCGCGATATCATTGGTAATAATGTGAAATCGCAGGGTAACAGGCCCGCCTAGTTTACACTCCACCGTATACTGATGGCGCAAAAtt                                                                                                                                                                                                     <  1:14821/246‑1 (MQ=255)
            gtaTAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACACCCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGAc                                                                                                                                                                                                  <  1:1069/249‑1 (MQ=255)
            gtaTAACAGCATGCCGATCGAATTGTTCCAGAAGCAGAAAAAAGAGATCCCGGACGAAGTTCACGTGGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGAc                                                                                                                                                                                                  <  2:117679/249‑1 (MQ=255)
                      aTGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGa                                                                                                                                                                                         >  1:264600/1‑250 (MQ=255)
                               gAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGCTCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGc                                                                                                                                                                               >  1:93557/1‑251 (MQ=255)
                                 aTTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTg                                                                                                                                                                             <  1:104075/251‑1 (MQ=255)
                                             gCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGaaa                                                                                                                                                                  >  1:13727/1‑250 (MQ=255)
                                                      agagATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACg                                                                                                                                                         >  1:134800/1‑250 (MQ=255)
                                                                 aCGAAGGTCCCGTTGATCCATACATGTGCACTTCCTATTCCGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGCTGGCGCACACTTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGaa                                                                                                                                               >  2:238997/1‑250 (MQ=255)
                                                                             ttgatgaaTACCTGTGCACTTACTATTACGCCATTAACAACCAGAAACCGCAATTCAAAGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGCATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAACACTgctggctg                                                                                                                                 >  2:35557/8‑251 (MQ=255)
                                                                                      aCCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAACCTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGtt                                                                                                                        >  2:127035/1‑251 (MQ=255)
                                                                                                     aTTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAAc                                                                                                         <  1:161343/251‑1 (MQ=255)
                                                                                                       tgcgatagtaccacacAGACACCGCCAGTCAACGATGGTCGTGTGCGTACCGCACTGAAACTATGTATGGACCGCGATATCATTTTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCCAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACaagcc                                                                                                        <  2:61383/235‑4 (MQ=255)
                                                                                                           gaaTTACCAACCCAGAATCCGCCATTCAACGATGTGCGGTTGCGTACCGCACTGAACCGAGCTATGGACCGCTATATCTTGGTTAATGAAGTGAAAGCGCAGGGCAAAATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGTCGCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAATATCTGCTGGCTGAAGCGGGTTATCCCGCAGACAAACCGTTg                                                                                                    <  2:170168/250‑1 (MQ=255)
                                                                                                                   aacCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAAc                                                                                           >  1:41125/1‑251 (MQ=255)
                                                                                                                                    tCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACAcc                                                                            <  2:129443/250‑1 (MQ=255)
                                                                                                                                                    gcgtACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATaaaa                                                            <  1:138596/250‑1 (MQ=255)
                                                                                                                                                          cGCACTGAAACTAGGAATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTgtc                                                     <  2:70254/251‑3 (MQ=255)
                                                                                                                                                                  aaCTAGGTATGGACCGCGCTATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGGTtgctg                                             >  2:124672/1‑251 (MQ=255)
                                                                                                                                                                            ggACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCAtt                                   <  1:127035/251‑1 (MQ=255)
                                                                                                                                                                                   gATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGga                              >  1:226186/1‑249 (MQ=255)
                                                                                                                                                                                       tCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGaaa                         >  1:243167/1‑250 (MQ=255)
                                                                                                                                                                                                              cagggcaaattcaccGCCTATGGGTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACt  <  2:78922/238‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
AATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAATTGACTCAGCCAGAATGGTTTGGCTGGAGCCAGGAAAAACGTAACGAAGAAGCGAAAAAACTGCTGGCTGAAGCGGGTTATACCGCAGACAAACCGTTGACCATCAACCTGTTGTATAACACCTCCGATCTGCATAAAAAGCTGGCGATTGCTGCCTCTTCATTGTGGAAGAAAAACATTGGTGTAAACGTCAAACT  >  NC_000913/1301979‑1302434

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