Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,477,878 T→A 33.3% intergenic (+20/+29) slyX → / ← slyD phi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,477,8780TA33.3% 8.7 / 9.0 12intergenic (+20/+29)slyX/slyDphi X174 lysis protein/FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base A (2/2);  total (6/6)
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.69e-02

TCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGT  >  NC_000913/3477777‑3478010
                                                                                                     |                                                                                                                                    
tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                                                               <  1:1045971/109‑1 (MQ=255)
tCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGc                                                                                                                               >  2:1045971/1‑109 (MQ=255)
                      gccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                             >  3:770350/1‑89 (MQ=255)
                      gccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAATCCCCGCtt                                                                                                                             <  4:770350/89‑1 (MQ=255)
                       ccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                      <  3:654102/79‑1 (MQ=255)
                       ccagccagccGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGAt                                                                                                                                      >  4:654102/1‑79 (MQ=255)
                                                                                             aGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt  >  1:948366/1‑141 (MQ=255)
                                                                                             aGCGGGGAATCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCgtcgtggt  <  2:948366/141‑1 (MQ=255)
                                                                                                 gggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCt                                                            >  2:956173/1‑79 (MQ=255)
                                                                                                 gggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGCTGCCTTTACCGCCACAGCACCCt                                                            <  1:956173/79‑1 (MQ=255)
                                                                                                  ggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCaa            <  1:1097255/126‑1 (MQ=255)
                                                                                                  ggATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCaa            >  2:1097255/1‑126 (MQ=255)
                                                                                                    aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCa                           >  1:123878/1‑109 (MQ=255)
                                                                                                    aTTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCa                           <  2:123878/109‑1 (MQ=255)
                                                                                                     |                                                                                                                                    
TCTGCTGACCGAGAAGTTAAAAGCCAGCCAGCCGTCGAACATCGCGTCGCAGGCTGAAGAAACGCCACCGCCACATTATTGAGGCGTAAAAAAAGCGGGGATTCCCCGCTTTTTTGTCACTTTTTCGGTATTAGTGGCAACCGCAACCGCCGTTGCCTTTACCGCCACAGCAGCCTTCGCCACCGTGTTCATGACCGTGATCGTGGCCATGACCGCCGCAGCAACCGTCGTGGT  >  NC_000913/3477777‑3478010

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