Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,301,992 A→T N271Y (AAT→TAT)  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,9920AT81.8% 44.2 / 2.6 22N271Y (AAT→TAT) oppAoligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (10/8);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

ATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAT  >  NC_000913/1301773‑1302238
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
aTCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGTAACAACGCGAAAACCGTTATTAAACAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAg                                                                                                                                                                                                                           <  1:244562/249‑1 (MQ=255)
                                                       ccTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAACAAGAGATCACGGACGAAGTTCACGTTGATCTCTACCTGTGCACTTAc                                                                                                                                                                  >  2:37630/1‑251 (MQ=255)
                                                                tAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACg                                                                                                                                                           <  1:310195/249‑1 (MQ=255)
                                                                                        gaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTc                                                                                                                                 >  1:241505/1‑251 (MQ=255)
                                                                                        aaaTCGTTCTTGAACGCAGCCCGCCCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTCCCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATCCCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTc                                                                                                                                 <  1:111000/250‑1 (MQ=255)
                                                                                                             cGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAACCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCTATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTAcc                                                                                                            <  1:37630/251‑1 (MQ=255)
                                                                                                              aaCCTACTGGAACACCGCTAAACCCGTTATTAACAAGGTAACCTTTTTGCCTATTGCTTCTGAAGTTACCGATTTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTACACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCg                                                                                                           <  1:251564/250‑1 (MQ=255)
                                                                                                                      gaaaaaacgctaaatccttgttaacccAGTTAACCTATTGGCCTATTGCTTTTGAATTTACCGATGTAAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACGGCATCCCGATCGAATTGTTCCAGAAGTGGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATTCCTGGTCACTTTCTATGCCGAATTTAACAACCAAAATCCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGaaa                                                                                                   <  2:341236/226‑1 (MQ=255)
                                                                                                                        aacaacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACt                                                                                                 >  1:282885/1‑251 (MQ=255)
                                                                                                                          caacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACt                                                                                                 >  2:17617/1‑249 (MQ=255)
                                                                                                                                  aacccGTTATTAACCAGGTAACCTATTTGCCTATTTCTTCTGAATTTTCCTAGGTCATCC‑GCTACCGTAGTGGGGAAATTGACATGACGTAAAACCGCATGCCGATCGAATTGTTCCATAAGCTGAAAAAAGAGATCCCGGACGAAGTCCACGTTGATCCATCCCTTTGCACTTTCTATTACGAAATTACCATCCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGAc                                                                                       <  1:215528/248‑1 (MQ=255)
                                                                                                                                         tattaACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATAc                                                                                                                                                                              <  1:8177/157‑1 (MQ=255)
                                                                                                                                         tattaACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATAc                                                                                                                                                                              >  2:8177/1‑157 (MQ=255)
                                                                                                                                              accagttaccttATTTGACTATTGCTTCCGAAGTTACCGATGTCATCC‑TATACCGTAGTGGTAAAACTGACTTGACGTATAACAGCATGCTGATCGAATTGTCCCAGAATCATAAAAAAGAGATCCCGGACGAAGTCCACGTTGATCCATACCTGTGCCTTTACTATTACGAAATTAACACCCAAATACCGCCAGTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCAtt                                                                           <  2:291749/240‑1 (MQ=255)
                                                                                                                                               ccAGATAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGT‑AACCATCTACCGTCGCGGTGAATTCTACATGAAGTATAACCGCATGCCGATCGCAATGTTCAAGAAGCTGAAAAAAGAGATCCCGGACGAAGGTAAAGTTGATCCATACCTGTGCACTTACTTTTACGAAATTAACAACAAGAAACCGCCATTAAACGATGTGCGAGGGCGTACCGAACTGAAACTAGGTAGGGACCGGGCTATCATTg                                                                          >  2:55404/1‑251 (MQ=255)
                                                                                                                                                                   ttGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGACCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCAGGACGCAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAACCTAGGTATGGACCGCGATCTCCTTGTTAATAAAGTGAAAtcgc                                                        >  2:21209/1‑249 (MQ=255)
                                                                                                                                                                          tGAAGTTACCGATGTCAACC‑GCTACCGTCGAGGTGACATCGAACTGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAAc                                                >  1:133599/1‑250 (MQ=255)
                                                                                                                                                                               ttACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGccc                                          >  1:173241/1‑251 (MQ=255)
                                                                                                                                                                                               gCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAAGGTTACACt                           >  1:100653/1‑251 (MQ=255)
                                                                                                                                                                                                           gtgAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACt               >  2:47854/1‑251 (MQ=255)
                                                                                                                                                                                                                    aCATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGgcgc       >  1:254468/1‑250 (MQ=255)
                                                                                                                                                                                                                        aaCGAATAACACAAGGCCGAGAGAAATGTCCCAGAAGCTAAAAAAAGAGTACCCTGACGAAGTCCACGTTGATCAATACCGGTGCACTTACCTTTACAAAAGTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTAAGGTTACACTCCGCCGTATACTGATGGCGCAAAAt  <  2:100653/250‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
ATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACC‑GCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAT  >  NC_000913/1301773‑1302238

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