Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,883 +C coding (667/687 nt) rph ← truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,8831.C94.4% 48.8 / ‑3.0 18coding (667/687 nt)rphtruncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (7/10);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

ATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCC  >  NC_000913/3815646‑3816115
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                        
aTTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTGCTACAGACAAAAAAAAGGCGACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCggg                                                                                                                                                                                                                                     <  1:171070/244‑1 (MQ=255)
                             ggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                                                                                                                                                                                  >  1:312992/1‑250 (MQ=255)
                                                   ngTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                                                                                           >  1:86852/2‑251 (MQ=255)
                                                   cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                                                                                                                                                                           >  1:179620/1‑251 (MQ=255)
                                                                                  ctAAGCGCAAATTCAAAAAACGGGCGCTGATATGGTATCATGCCCTCGCTCCTCATCTTACTTTGCTACAGAAAAAAAAAAGTCGACTCGTCAGTCGCCTGAAAAATCAGCTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCGAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                                                                                                                                                                 <  1:264052/214‑1 (MQ=255)
                                                                                     aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGtcatc                                                                                                                                         <  1:201962/251‑1 (MQ=255)
                                                                                                                 atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCg                                                                                                                                                                               <  1:271459/185‑1 (MQ=255)
                                                                                                                                    cctcATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTACCAATAAGTTTGCCCGCGACGCCTCCTGCGTCGATACAATGGATTCCATTCC‑CCCTCGGGCCCGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACt                                                                                                                                                                    >  1:175157/1‑177 (MQ=255)
                                                                                                                                      tcATCTTACTTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGACTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGATCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgatg                                                                                                                                                           <  1:53279/184‑1 (MQ=255)
                                                                                                                                            tACTTTGCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGGCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCg                                                                                      <  1:118115/247‑1 (MQ=255)
                                                                                                                                                     aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAGTCAGTTCGCCAGCGCCGCCTTCTGCGTCGCTTAAAGGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGATGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTCCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCATAGTCTTCAACGTATTCCAGATCGCAAACCGCTTc                                                                           <  1:28329/249‑1 (MQ=255)
                                                                                                                                                                     cGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCCTCTGCGTCGCTCCAATGGATTCGATTCCACC‑TCGGGCCAGAGCCACCAAGATGAGTAGCTCTTCCTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTACTGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGAc                                                         >  1:49795/1‑251 (MQ=255)
                                                                                                                                                                                              aGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcgcn                                                                                                                                                       <  1:86140/132‑2 (MQ=255)
                                                                                                                                                                                                         nccgccTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGCTGAGTAGCTATCCATAGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCACCGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGt                     >  1:69557/2‑251 (MQ=255)
                                                                                                                                                                                                            nccTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                  >  1:70234/2‑155 (MQ=255)
                                                                                                                                                                                                             ccTTCTGCGTCGCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCGCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTc                 <  1:176443/251‑1 (MQ=255)
                                                                                                                                                                                                                        gCTACAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGCTCGCAAACCGCTTCGCCGTTCACAATTCCGACAGACACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGtt      >  1:133975/1‑251 (MQ=255)
                                                                                                                                                                                                                           aCAATGGATTCGATTCC‑CCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAGCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTGCGCCGTTCACAATTC‑GACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTcc  <  1:356572/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                        
ATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCC‑CC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCC  >  NC_000913/3815646‑3816115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑