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

ATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAT  >  NC_000913/1301773‑1302238
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
aTCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGTAACAACGCGAAAACCGTTATTAAACAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAg                                                                                                                                                                                                                           <  1:244562/249‑1 (MQ=255)
                                                       ccTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAACAAGAGATCACGGACGAAGTTCACGTTGATCTCTACCTGTGCACTTAc                                                                                                                                                                  >  2:37630/1‑251 (MQ=255)
                                                                tAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACg                                                                                                                                                           <  1:310195/249‑1 (MQ=255)
                                                                                        gaaTCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTc                                                                                                                                 >  1:241505/1‑251 (MQ=255)
                                                                                        aaaTCGTTCTTGAACGCAGCCCGCCCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTCCCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATCCCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTc                                                                                                                                 <  1:111000/250‑1 (MQ=255)
                                                                                                             cGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAACCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCTATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTAcc                                                                                                            <  1:37630/251‑1 (MQ=255)
                                                                                                              aaCCTACTGGAACACCGCTAAACCCGTTATTAACAAGGTAACCTTTTTGCCTATTGCTTCTGAAGTTACCGATTTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTACACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCg                                                                                                           <  1:251564/250‑1 (MQ=255)
                                                                                                                      gaaaaaacgctaaatccttgttaacccAGTTAACCTATTGGCCTATTGCTTTTGAATTTACCGATGTAAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACGGCATCCCGATCGAATTGTTCCAGAAGTGGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATTCCTGGTCACTTTCTATGCCGAATTTAACAACCAAAATCCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGaaa                                                                                                   <  2:341236/226‑1 (MQ=255)
                                                                                                                        aacaacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACt                                                                                                 >  1:282885/1‑251 (MQ=255)
                                                                                                                          caacGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACt                                                                                                 >  2:17617/1‑249 (MQ=255)
                                                                                                                                  aacccGTTATTAACCAGGTAACCTATTTGCCTATTTCTTCTGAATTTTCCTAGGTCATCCGCTACCGTAGTGGGGAAATTGACATGACGTAAAACCGCATGCCGATCGAATTGTTCCATAAGCTGAAAAAAGAGATCCCGGACGAAGTCCACGTTGATCCATCCCTTTGCACTTTCTATTACGAAATTACCATCCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGAc                                                                                       <  1:215528/248‑1 (MQ=255)
                                                                                                                                         tattaACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATAc                                                                                                                                                                              <  1:8177/157‑1 (MQ=255)
                                                                                                                                         tattaACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATAc                                                                                                                                                                              >  2:8177/1‑157 (MQ=255)
                                                                                                                                              accagttaccttATTTGACTATTGCTTCCGAAGTTACCGATGTCATCCTATACCGTAGTGGTAAAACTGACTTGACGTATAACAGCATGCTGATCGAATTGTCCCAGAATCATAAAAAAGAGATCCCGGACGAAGTCCACGTTGATCCATACCTGTGCCTTTACTATTACGAAATTAACACCCAAATACCGCCAGTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCAtt                                                                           <  2:291749/240‑1 (MQ=255)
                                                                                                                                               ccAGATAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTAACCATCTACCGTCGCGGTGAATTCTACATGAAGTATAACCGCATGCCGATCGCAATGTTCAAGAAGCTGAAAAAAGAGATCCCGGACGAAGGTAAAGTTGATCCATACCTGTGCACTTACTTTTACGAAATTAACAACAAGAAACCGCCATTAAACGATGTGCGAGGGCGTACCGAACTGAAACTAGGTAGGGACCGGGCTATCATTg                                                                          >  2:55404/1‑251 (MQ=255)
                                                                                                                                                                   ttGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGACCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCAGGACGCAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAACCTAGGTATGGACCGCGATCTCCTTGTTAATAAAGTGAAAtcgc                                                        >  2:21209/1‑249 (MQ=255)
                                                                                                                                                                          tGAAGTTACCGATGTCAACCGCTACCGTCGAGGTGACATCGAACTGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAAc                                                >  1:133599/1‑250 (MQ=255)
                                                                                                                                                                               ttACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGTATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGccc                                          >  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)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      
ATCGAGAAATTCGGCGAAAAATGGACCCAGCCTGGTAATATCGTCACCAACGGTGCCTATACCTTAAAAGATTGGGTCGTAAACGAACGAATCGTTCTTGAACGCAGCCCGACCTACTGGAACAACGCGAAAACCGTTATTAACCAGGTAACCTATTTGCCTATTGCTTCTGAAGTTACCGATGTCAACCGCTACCGTAGTGGTGAAATCGACATGACGAATAACAGCATGCCGATCGAATTGTTCCAGAAGCTGAAAAAAGAGATCCCGGACGAAGTTCACGTTGATCCATACCTGTGCACTTACTATTACGAAATTAACAACCAGAAACCGCCATTCAACGATGTGCGTGTGCGTACCGCACTGAAACTAGGTATGGACCGCGATATCATTGTTAATAAAGTGAAAGCGCAGGGCAACATGCCCGCCTATGGTTACACTCCACCGTATACTGATGGCGCAAAAT  >  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: