Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,269,250 A→G I9T (ATT→ACT)  ldrA ← small toxic polypeptide LdrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,269,2500AG94.7% 37.9 / ‑3.4 19I9T (ATT→ACT) ldrAsmall toxic polypeptide LdrA
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (8/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01

GCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTT  >  NC_000913/1269012‑1269331
                                                                                                                                                                                                                                              |                                                                                 
gCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACt                                                                       >  2:52655/1‑251 (MQ=255)
gCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACt                                                                       >  1:393449/1‑251 (MQ=255)
      aaTCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCAGGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGGTTCGGACGCATAGTGCCATCACGGCAGCCTTACGCCCGTTTTCACCTTACTTGCGGTTACGCCACCAGCTGACAAACGCTGCGGACATAATTCACACAAGGATCGGGGCTGCTAGGACGAGCCAGAAAGTGAAGGCAAACTGagcgg                                                                 >  2:300761/1‑250 (MQ=255)
             gttGCTTGAAAATTAAAGTCTTGGCGTCTTTCGTTAAGATTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGATGGCTGTTACCCCGGCGGCAGGTTGCGGGCGAATATTGCCATCCCGGCAGCCGGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGATGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCa                                                           <  2:311071/249‑1 (MQ=14)
                 cTTGAAAAATAAAGTATTAGCGTTCTGAGTTAAGACTTTTTTCATGGGAGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTTCGGCAGGTTTCGGGCGCATATTGCCATCACGGCAGCATGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTAACAAACGCTCCCGTAATAATTCCCGCCAGGTTCGGTGCTACCAGGTCGTGCCAGATAGTCATGGCAAACTGCGCGAGCGTAAtatac                                                      >  1:321768/1‑250 (MQ=17)
                      aaaaTAAATTATTAGCTTTCTGCGTTAAGTCTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAgccgcc                                                 <  1:462517/251‑1 (MQ=255)
                      aaaaTAAAGTATTAGCTTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAgccgcc                                                 <  1:141927/251‑1 (MQ=255)
                        aaaaaaggattgccgtgcTGCGTAATGACTTTTGGCAGGGGTGCTGTATTCAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCGGGTTTCGGCCGCATATTGCCATCCCGGCAGCCTGACGCCCGTTTTCACCTTCCTTCCGGTGAGGCCACCAGCGGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCtt                                               <  2:187681/234‑1 (MQ=11)
                        aaTAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGGGCGTCCTATAGCCGCCtt                                               >  2:116383/1‑251 (MQ=255)
                         aTAAAGGATGAGCGTTCTGTGTTAAGACTTTTTTCATGGGGGCGGGATACAAAAAAGTCCGCAGGCGGTTTCCCCGGCGGCGGGTTTGGGGCGCAGATGTCCACCACGGAAGCCAGACGCCGGTTGTCTCCTTACTTCCGGTGCCGCCACCGGTTGATAATCGCTGCGGTAATAATTCCCGCCAGTATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCttg                                              <  1:30599/251‑1 (MQ=9)
                           aaaggatgagcgtgcggggTTAAGTCTTTTTTCAGGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCTAATATTGCCATCCCGGCAGCCTGACGCCCGTTTTCCCCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATGCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCTTCATATAGCCGCCttgt                                             <  2:430477/232‑1 (MQ=14)
                                         cTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACAtt                              >  1:572521/1‑251 (MQ=255)
                                           tggttaaggtttttttCTGGGGTGCCGGATACAAAAAAGGCCGCAGGCTTTTCCCTCTGGGGCCGGTTTCGGGGGCATATTGCCATCACGGCAGCTGGACGCCCGTTTTCCCCTTACTCCCGGTTACGCCACCAGCTGAAAAGCGCTGCGGTATAAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCATAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACAt                               <  1:266110/238‑1 (MQ=11)
                                            ggttaaggcgttttccttGGGGGCCGGATACAAAAAAGGCCGCAGGCTGTTTCCCCTGCGGCCGGTTTCGGGCGAATATTGCCATCCCGGCAGCCTGACGCCCGTTTTCCCCTTCCTTCCGGTCACGCCACCAGTTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTAATGGCAAACGGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTg                           <  2:194822/234‑1 (MQ=11)
                                            cGTTAAGACTTTTTTCAGGGGTGCCGGATACAAAAAAGGCCGCAGCCTGTTACCCCTCCGGCCGGTTTCGGGCGCATATTGCCAGCACGGCAGCCTGACGCCCGTTTTCACCTTACTCCCGGTTACGCCACCAGCTGACAATCGCTGCGTTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTg                           <  1:194274/251‑1 (MQ=14)
                                             gTTAAGACTTTTTTCAGGGGTGCCGGATAAAAAAAAGGCCGCAGGCTGTTACCCCGGCGGCGGGTTTCGGGCGCATTGGGCCACCACGGCAGCCTGACGCCCGTTTCCCCCTTACTCCCGGTTACGCCACCAGCTGACAATCGCGGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGc                          <  2:572521/251‑1 (MQ=14)
                                               tatgtcTTTTTGCATGGTCGCCGGATAAAAAAAAGGCCGCAGGCTGTTCCCCCTGGGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCCCTTTCCTTCCGGTTACGCCACCAGCTGCCAATCGCGGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCATGTCGTTTCAGAAAGGCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCg                         <  1:362584/245‑1 (MQ=11)
                                                          tCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGa              >  1:505172/1‑250 (MQ=38)
                                                                     ggATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATAGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCAAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCAGGCTTGTTGTAATGACAANATTTTGCGGCTATTCTTggagtggctggttt  >  2:266260/1‑238 (MQ=11)
                                                                                                                                                                                                                                              |                                                                                 
GCAAGTAATCTTTTTTGCTTGAAAAATAAAGTATTAGCGTTCTGCGTTAAGACTTTTTTCATGGGTGCCGGATACAAAAAAGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGGTT  >  NC_000913/1269012‑1269331

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