Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 178,334 C→T T43T (ACG→ACA btuF ← vitamin B12 ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913178,3340CT90.0% 21.2 / ‑5.7 10T43T (ACG→ACAbtuFvitamin B12 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/5);  minor base G (1/0);  total (5/5)
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

AGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTGGAGGATAGTCGGAATAGCTGCTGACCCCAACCGGCGTGATCCCGGCGGCAAAGGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTG  >  NC_000913/178090‑178550
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
aGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAAGCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCGCAATCTGTTGCGCTTGTGGAGGATAGTCGGAATAGTTGCTGACCCCAACCGGTGTGATc                                                                                                                                                                                                                    <  1:94180/251‑1 (MQ=255)
                                                   cGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGACAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGAGCGGGTGTCAGCGCGACAAGGCTGTCCAGAGTCATCCCCGGCCAGGTGGAAACCTGCTCAATCTTTTGCGCATGTGGAGGATAGTCGGAATAGTTGCTGACCCCAACCGGTGTGATCCCGGCGGCAAAGGCAAGTTCAGTGTTGGCTGGAGAAAGCGTGATGAcgcgc                                                                                                                                                                 <  1:110686/251‑1 (MQ=255)
                                                              ngtnngaatNCCAGCGAGGCCTGCGGGTCAACCGNCCGCTCGGCATTACCTCCGCGCCAGGCAATCNCNAGATCGGGTTTCAGCGCGACAATGCGTTCCAGAGTCATCCCCTGCCAGNTGGAAACCTGNTCANNCTGTNGNGNTTGTGGATGATAGTCGGAATAGTTGCTGANCCCAACCGGTGTGATCCCGGCGGNNAANNNAANNTCAgc                                                                                                                                                                                             <  2:127601/204‑2 (MQ=255)
                                                                        ccAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGAAATCACCAGATCGGGTTTCAGCGCGCCAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTGGAGGATAGTCGGAATAGTTGCTGACCCCAACCGGTGTGATCCCGGCGGCAAAGGCACGTTCAgc                                                                                                                                                                                             >  1:127601/1‑201 (MQ=255)
                                                                                                                         gtcAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTGGAGGATAGTCGGAATAGTTGCTGACCCCAACCGGTGTGATCCCGGCGGCAAAGGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGcc                                                                                            >  1:161724/3‑250 (MQ=255)
                                                                                                                                aCCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTGGAGGATAGTCGGAATAGTTGCTGACCCCAACCGGTGTGATCCCGGCGGCAAAGGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCa                                                                                     >  1:60696/1‑250 (MQ=255)
                                                                                                                                                                                                                 ggaggaTAGGCGGAATAGTTGCTGACCCCAACCGGTGTGATCCCGGCGGCAAAGGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGTGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCATCAGCCAGGAACTCATCGAAGCTAAGATGAGACt   <  1:151413/251‑1 (MQ=255)
                                                                                                                                                                                                                  gaggaTAGTCGGAATAGTTGCTGACCCCAACCGGTGTGATCCCGGCGGCAAAGGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTg  >  1:118871/1‑251 (MQ=255)
                                                                                                                                                                                                                  gaggaTAGTCGGAATAGTTGCTGACCCCAACCGGGGAGATCCCGGCGGCAAAGGCAAGTTCAATGTTGGCGGGAGAAAGAGTGATGACGCGCGGCGCGGCTTTGAGCCACAGTGGCGCACGAAAAGACAGGACGACCAGCGCCCTGAACAGTGACTTAACCATGTGCAAGTTTCTGCACCAGTGACTCTACCATCAGGCTGGACTGATAAGCGGCAACAGCCAGGACCTCATCGACGATAAGATGAGAcgg  >  1:118866/1‑249 (MQ=255)
                                                                                                                                                                                                                                          ccccAACCGGTGTGATCCCGGCGGCAAAGGCAAGTTCATTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGTGGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACATTGACTTAGCCATGTGCAAGTGGCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGc                                         <  1:93602/188‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        
AGGGGGCCAGTTGACGTAACGCATTGGCAATTTGTTCAATGCTTGTCGCATCGACCCACATCACTTTTATTCCCAGCGAAGCCAGCTGGTCAACCTGCCGCTCGGCATTACCTCCACGCCAGGCAATCACCAGATCGGGTTTCAGCGCGACAATGCGTTCCAGATTCATCCCCTGCCAGGTGGAAACCTGCTCAATCTTTTGCGCTTGTGGAGGATAGTCGGAATAGCTGCTGACCCCAACCGGCGTGATCCCGGCGGCAAAGGCAAGTTCAGTGTTGGCGGGAGAAAGCGTGATGACGCGCGGCGCGGCGTTGAGCCACAGTGGCGCAAGAAAAGACAGGGCGACCAGCGCCCTGAACAGTGACTTAGCCATGTGCAAGTTTCTGCACCAGTGACTCAACCATCAGGCTGGACTGTTTAGCGGCAACAGCCAGGAACTCATCGAAGCTAAGATGAGACTG  >  NC_000913/178090‑178550

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