Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,306,736 T→G S325A (TCC→GCC)  oppF → murein tripeptide ABC transporter/oligopeptide ABC transporter ATP binding subunit OppF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,306,7360TG95.5% 57.5 / ‑6.1 22S325A (TCC→GCC) oppFmurein tripeptide ABC transporter/oligopeptide ABC transporter ATP binding subunit OppF
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (9/12);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGGTAGAACTGGGGACCTATGATGAGGTCTACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCA  >  NC_000913/1306490‑1306968
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
gCGGTAGAACGGGGGACCTATGATGAGGTCTACCCCAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTCCCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAACCACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGt                                                                                                                                                                                                                                      <  2:93083/251‑1 (MQ=255)
     aGAACTGGGGACCTATGATGAGGTCTACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCCGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCtt                                                                                                                                                                                                                                 >  1:29280/1‑251 (MQ=255)
               aCCTATGATGAGGTCTACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTTTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACTGCCTGTTCTGGAGGGGGGTTTCAGACACGCCGTTTCTTGCCTGAAAGt                                                                                                                                                                                                                       <  1:27564/251‑1 (MQ=255)
                                 cacaATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACa                                                                                                                                                                                                     <  1:436026/251‑1 (MQ=255)
                                 cacaATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACa                                                                                                                                                                                                     <  1:393876/251‑1 (MQ=255)
                                       ccACTACATCCTTACCCCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCACCCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAAACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAgg                                                                                                                                                                                               <  1:9579/251‑1 (MQ=255)
                                       ccACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAgg                                                                                                                                                                                               >  1:332025/1‑251 (MQ=255)
                                                ccTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGAGCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCGTTCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGt                                                                                                                                                                                      <  1:440474/251‑1 (MQ=255)
                                                                 gATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAACCCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGACGGTCCAGAGTGCGCCAACACACGTCCTGTTCTGGAGGGGAGTTTCAGACCCGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCAc                                                                                                                                                                     >  1:282144/1‑251 (MQ=255)
                                                                                 agaCCTGATCGGGATCTGGAGAAGAACAAAACCAGCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGGGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCGTGCCTGAAAGTTGATCCGCTGTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattatt                                                                                                                                                     <  2:29280/249‑1 (MQ=255)
                                                                                 aTACCTGATCCGGATCTGGAGAAGAACAAAACCATCCCGTTACTGGAAGGAGAATTACCGTCGCCGCTCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCACGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGCCACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGCGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGattatt                                                                                                                                                     >  1:381488/1‑251 (MQ=255)
                                                                                        aTCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTccc                                                                                                                                              >  2:84825/1‑251 (MQ=255)
                                                                                                                     cAGTTACTGGAAGGGGAATTACCGTCGCCGATCACCCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCt                                                                                                                 <  2:188025/251‑1 (MQ=255)
                                                                                                                            tGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTAAAAGTTGCCCGATTGCCGGTCCAGAGTGCGCCAACACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTcgc                                                                                                          >  1:280722/1‑251 (MQ=255)
                                                                                                                              gAAGGGGAATTACCGTCGCCGATCACCCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACCCGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCgg                                                                                                        <  1:246048/251‑1 (MQ=255)
                                                                                                                                      tttcccgggggcGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGCAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAACAGCGa                                                                                                                                                          >  2:29158/12‑193 (MQ=255)
                                                                                                                                                ccGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGACGGTGCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGGGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGAt                                                                                                                                                                                                                   >  1:458819/1‑126 (MQ=255)
                                                                                                                                                         ccGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCg                                                                             >  1:63798/1‑251 (MQ=255)
                                                                                                                                                                        gtTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCGGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATc                                                               <  1:118159/250‑1 (MQ=255)
                                                                                                                                                                          ttcccGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCtt                                                             <  1:34316/247‑1 (MQ=255)
                                                                                                                                                                          ttaccGTACCCGTTGCCCGCTTGCCGGTCCAGAGTGCGCCAATACATGTCCTGTTCTGGAGGGGAGTTTCCGACACGCCGTTTCTTGCCTGAAAGACGATCCGATTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTAACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTATGATCTTTTTAGCGGCATAACAAAATCCGGACAAAAACTTCGTTGATTTCACCATCTTc                                                            >  2:194373/4‑251 (MQ=255)
                                                                                                                                                                                                                                    gAGGGGAGTTTCAGACACGCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTGTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTTATTTCACCATCTTCACTGTCCGCCAGATCAATCACCACTTCGGCAAAAAATTCCGGGTTCAAATCGAAATCa  <  2:26559/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        
GCGGTAGAACTGGGGACCTATGATGAGGTCTACCACAATCCACTACATCCTTACACCAGGGCATTGATGTCGGCAGTCCCCATACCTGATCCGGATCTGGAGAAGAACAAAACCATCCAGTTACTGGAAGGGGAATTACCGTCGCCGATCAACCCGCCTTCCGGTTGTGTTTTCCGTACCCGTTGCCCGATTGCCGGTCCAGAGTGCGCCAAAACACGTCCTGTTCTGGAGGGGAGTTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCA  >  NC_000913/1306490‑1306968

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