Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 855,421 C→T M108I (ATG→ATA ldtB ← L,D‑transpeptidase LdtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913855,4210CT80.0% 16.1 / ‑5.0 10M108I (ATG→ATAldtBL,D‑transpeptidase LdtB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (3/5);  minor base A (2/0);  total (5/5)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACAGCCCCATCGGGTTATCCGGACCTGCCGGAACGACAGCCGGAAGCGGTTCGCCCGCAGCGCGGTACTCTGCGTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTGGAAAGCCCCATCTGGTACTCGGCGGCAAAATACTCCAGCGGCTGAGTGTTACCTTCAGGAATGGTGATCACCTGATTCTGACCAACCAGGCGACTCCCGTCGGTTGGCAGAGGA  >  NC_000913/855196‑855664
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
aTACAGCCCCATCGGGTGATCCGGACCGGCCGGAACGACAGCCGGAAGCGGTGCGTTCGCAGCGCGGTACTCTGCGTGCATTTTGGCGGTCTGCGTCCAGGTCGTTCCTGCTTTTTGACGCTCATCTGTGGTGGTGCAGTTGATAGGCGTATCTTGGTCTAACTGACCAATGTCGATCGTAAGCGCGATAACGTTGTTGGTCCCTTTCGGATAGTAATAAAGACGTATCTCAGCACTGTTAATGTCGATGc                                                                                                                                                                                                                            <  1:186972/251‑1 (MQ=255)
                     cGGACCTGCCGGAACGACAGCCGGAAGCGGTTCGCTCGCAGCGCGGTACTCTGCGTGCATTTTGGCGGTCGGTGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTAGCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACTGTTGTGGTCCCTTTCGGATAGTAATAAAGACGTATCTCAGc                                                                                                                                                                                                                                             <  1:222893/213‑1 (MQ=255)
                                          cGGAAGCGGTTCGCTCGCAGCGCGGTACTCTGCGTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGTATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGt                                                                                                                                                                                  >  1:8852/1‑251 (MQ=255)
                                                                                                                                tgttggtCCAGTGGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGTATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATcc                                                                                                                                                                                                           <  1:193918/137‑1 (MQ=255)
                                                                                                                                                            nccTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGTATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATcc                                                                                                                                               >  1:24948/2‑172 (MQ=255)
                                                                                                                                                             ccTAACTGACCAATGCCGATAGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATCGTAATAAAGACGAATCTCAGCACTGTTAATGACGATg                                                                                                                                                                                                                             >  1:132845/1‑93 (MQ=255)
                                                                                                                                                                                   ncagccCNNTAACGGTGTTGGTACCTTTCGGATAGTAAGAAAGACGTATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATTCGGCAGGAt                                                                                                                                                                                                   <  2:17665/96‑1 (MQ=255)
                                                                                                                                                                                                ggTGTNGGTCANTNNCGGATAGTAAAAAAGACGTNTCTCAGCACTGTTAANGACGNTGCCTNCANg                                                                                                                                                                                                                     >  2:97239/1‑66 (MQ=255)
                                                                                                                                                                                                                           aaGACNAATCANANNACTGTTAATGCCGATGCCTNCATGAACGGTATCCGNCAGGNTCAgcngc                                                                                                                                                                                            >  2:100774/1‑64 (MQ=255)
                                                                                                                                                                                                                           aaGACGTATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTGGAAAGCCGCATCTGGTACGCGGCGGCAAAATGCTCCAGCGGCTGAGTGTTACCTTGAGGAATGGTGATCACCTGATGCTTACCAACCAGGCGACTCCCGTCGGTTGGCAGAGGa  <  1:100774/250‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
ATACAGCCCCATCGGGTTATCCGGACCTGCCGGAACGACAGCCGGAAGCGGTTCGCCCGCAGCGCGGTACTCTGCGTGCATTTTGGCGGTCGGCGTCCAGGTCGGGCCTGCTTTTTTACGCTCAACTTTGGTGGTCCAGTTGATAGGCGTATCTTTGCCTAACTGACCAATGCCGATCGGCAGCACGATAACGGTGTTGGTCCCTTTCGGATAGTAATAAAGACGCATCTCAGCACTGTTAATGACGATGCCTTCATGAACGGTATCCGGCAGGATCAGCTGCTGCGGAATGTTCAGTACAGTACCGCCTTTCGGCAGGAAGGTATCCACACCCGGGTTCGCTTCCATCATATTGGAAAGCCCCATCTGGTACTCGGCGGCAAAATACTCCAGCGGCTGAGTGTTACCTTCAGGAATGGTGATCACCTGATTCTGACCAACCAGGCGACTCCCGTCGGTTGGCAGAGGA  >  NC_000913/855196‑855664

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