Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,982,842 G→T intergenic (+31/‑116) proM → / → aslB tRNA‑Pro/putative anaerobic sulfatase maturation enzyme AslB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,982,8420GT81.2% 41.3 / 4.4 16intergenic (+31/‑116)proM/aslBtRNA‑Pro/putative anaerobic sulfatase maturation enzyme AslB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (2/11);  minor base . (2/1);  total (4/12)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

TTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAAC  >  NC_000913/3982706‑3982967
                                                                                                                                        |                                                                                                                             
ttGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGtttttttttt                                                                                                                            <  2:50194/140‑1 (MQ=255)
           aTTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCAttt                                                                                                                 <  1:266944/140‑1 (MQ=255)
                     cAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTCTGTGCATTTCGTCAccctcc                                                                                                      <  2:389935/140‑1 (MQ=255)
                                ggAGTAGCGCAACTTGGGAGCGCAAATGGTTTGGGACCAGTGGGTCGGGGGTTCCAATCCTCTCTCGCCGGCCAATTTTGAGCCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATa                                                                                            <  2:394953/139‑1 (MQ=255)
                                     aGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTTTCTGCATTTTGTCACCCTCCCTTCCCCATAAACGc                                                                                       >  2:428622/1‑140 (MQ=38)
                                      gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGcc                                                                                      <  2:341572/140‑1 (MQ=255)
                                      gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGcc                                                                                      <  1:428622/140‑1 (MQ=255)
                                                    cgcAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTg                                                                       >  1:126355/1‑140 (MQ=255)
                                                       aaCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGcc                                                                     <  2:58491/140‑1 (MQ=255)
                                                        aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCACAATAAACGCCCGTAATAACTCATTGccc                                                                    <  2:70196/140‑1 (MQ=255)
                                                          tGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCa                                                                  <  2:80267/140‑1 (MQ=255)
                                                               tGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTa                                                             <  2:239834/140‑1 (MQ=255)
                                                                                             ctctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTTCCCCTTAACGTATTGAAGGATGa                               >  2:190795/1‑140 (MQ=255)
                                                                                                     gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGc                       <  1:98613/140‑1 (MQ=255)
                                                                                                               gAACCCCGCTTCGGCGGGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGAc             <  2:252770/140‑1 (MQ=255)
                                                                                                                     cGCTTCGGCGGGGTTTTTTTTTTTCTGTGGATTTCTGCACCCCCCCTTCGCAATAAACGGCCGTAATAACTCCTTGCCCCACGGTTTTATTTCCCCCTTAAACGATTGAAGGATGACTTCAGGCAAGGAGCGACCaggtt       >  2:40677/1‑136 (MQ=255)
                                                                                                                     cGCTTCGGCGGGGTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAtgctg      >  2:445304/1‑140 (MQ=255)
                                                                                                                          cggcggGGTTTTTTTTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAAc  <  1:453777/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                             
TTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAAC  >  NC_000913/3982706‑3982967

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: