Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,058,473 G→A D23N (GAC→AAC)  bipA → 50S ribosomal subunit assembly factor BipA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,058,4730GA92.9% 32.2 / NA 14D23N (GAC→AAC) bipA50S ribosomal subunit assembly factor BipA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (6/7);  minor base C (1/0);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAA  >  NC_000913/4058258‑4058664
                                                                                                                                                                                                                       |                                                                                                                                                                                               
gACCTTTGGATCACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGAGCGATTAATCAGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCAACCTGGTAAACAAGCTGCTCCAACAATACGGTACGTTCGActct                                                                                                                                                              >  1:58619/1‑251 (MQ=255)
                    tATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAAGTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAgc                                                                                                                                           <  1:410944/250‑1 (MQ=255)
                              gAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGAc                                                                                                                                >  1:37482/1‑251 (MQ=255)
                                         caAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCGAATGCCTGAGGCGAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCCCCCTGGTAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGAttt                                                                                                                     <  1:433952/251‑1 (MQ=255)
                                                                                 gATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGGTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAGCTGCTCCACCAATCCGGTACGTTCTACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCCAAAAACAcc                                                                             <  1:3872/251‑1 (MQ=255)
                                                                                         tttGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCGCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGTAGAAAGAGCGTGGGATTACCATc                                                                                            <  1:149315/228‑1 (MQ=255)
                                                                                            ntAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATg                                                                  >  1:278229/2‑251 (MQ=255)
                                                                                                  cGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAACTTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAgctgct                                                                                                                                                                                       >  1:63028/1‑128 (MQ=255)
                                                                                                          ccGTGTACAATAACGCGCTATTTCTAATGCCTGAGTCAAAGTGGGGATCGAAAAAGTGCGTAAGATCGCCATCATCGCGCCCGTAGACCATGTTAAAACCACCCTGGAAAACAAGCTGCTCCAACAATCCGGTACTTTCGACTCTCGTGCCGAAACCCAAAAGCGCGTGATGGACTCAAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGc                                                                           <  1:249980/228‑1 (MQ=255)
                                                                                                           cGTGTACAATAAAGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGAAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGATATCAAATGGAATGATTACCGTAt                                                   >  1:320123/1‑251 (MQ=255)
                                                                                                                                 cTAATGCCTGAGGCAAAGTTGTGATCGAAAAAGTGCGTAATATCGCCATCATCGCGCATGTATACCATGGTAAAACCACCCTGGTAAACATGCTGCTCCAACAATCCGGTACGt                                                                                                                                                                      <  1:374899/114‑1 (MQ=255)
                                                                                                                                      gCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATAGCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACg                         <  1:353672/250‑1 (MQ=255)
                                                                                                                                                 aGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAAACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCATGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGAtt                                                                                                  >  1:200888/1‑166 (MQ=255)
                                                                                                                                                            aaaaaTTGCGTAATATCGCCATCATCGCGCACGTAGACAATCGTAAAACAACCCTGGCACCCAAGCGGCCCCCCCAATCCGGAAAGTTCGACTCTCGCGCCGCAACCCAAGAGCGCGTGATGGTCTCCACAGATATGGTGAAAGCGCGTGGGATTACCATCCTCGCGAACAACACCACTATTAAATGGCATGATTACCGTATCACCATCGATGATACCCCGGGGCACGCCCACTTCggcggccgaaacatg  >  1:216483/1‑238 (MQ=255)
                                                                                                                                                                                                             cnccNTGGTANACAANCTGCTCCAACAAGCCNGTACGTTCGACTCTCGTGCNNANTCCCANGAgcg                                                                                                                                          <  2:190166/66‑1 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                               
GACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAA  >  NC_000913/4058258‑4058664

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: