Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,067,300 C→T G8R (GGA→AGA)  insH‑6 ← CP4‑44 prophage; IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,067,3000CT97.0% 108.1 / ‑4.1 33G8R (GGA→AGA) insH‑6CP4‑44 prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/14);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

TTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAGCGC  >  NC_000913/2067224‑2067434
                                                                            |                                                                                                                                      
tttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAAcacc                                                                         >  2:756515/1‑140 (MQ=21)
         tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCg                                                                <  2:226163/140‑1 (MQ=255)
           aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGct                                                              >  2:375911/1‑140 (MQ=255)
           aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGct                                                              <  2:790144/140‑1 (MQ=255)
             ggCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctct                                                            <  2:46129/140‑1 (MQ=255)
               cgcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGctctct                                                          <  2:116965/140‑1 (MQ=255)
                gcTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCtctctc                                                         >  2:1093577/1‑140 (MQ=255)
                   tACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGt                                                      <  1:1453542/140‑1 (MQ=255)
                    aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTa                                                     <  2:914022/140‑1 (MQ=255)
                        ctgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGc                                                 <  2:58478/140‑1 (MQ=255)
                         tgAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCa                                                >  1:1281847/1‑140 (MQ=255)
                            aTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCa                                                <  2:753566/137‑1 (MQ=255)
                            aTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCa                                                >  1:753566/1‑137 (MQ=255)
                             ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGt                                            >  2:298713/1‑140 (MQ=255)
                                aCTGTCGGCGAAGGTAAGTTGATTACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTcac                                         >  1:356219/1‑140 (MQ=255)
                                aCTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTcac                                         >  1:1392689/1‑140 (MQ=255)
                                  tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTca                                          >  1:278129/1‑137 (MQ=255)
                                  tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTca                                          <  2:278129/137‑1 (MQ=255)
                                            ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAg                             >  1:1215754/1‑140 (MQ=255)
                                                     atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACt                    >  1:607249/1‑140 (MQ=255)
                                                      tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACtt                   >  1:691878/1‑140 (MQ=255)
                                                       gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCGCTGTCACAGCATTCAACAGAATGTGACGTg                  >  2:930986/1‑140 (MQ=255)
                                                       gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTg                  <  2:1392689/140‑1 (MQ=255)
                                                         cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGcc                <  2:682653/140‑1 (MQ=255)
                                                         cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGcc                <  2:737141/140‑1 (MQ=255)
                                                         cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGcc                <  1:1507313/140‑1 (MQ=255)
                                                         cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGcc                <  1:1034990/140‑1 (MQ=255)
                                                            atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGa             >  2:609283/1‑140 (MQ=255)
                                                               atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGAttt          >  2:687837/1‑140 (MQ=255)
                                                                    cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAg     >  1:318645/1‑140 (MQ=255)
                                                                      cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAgcg   >  1:310079/1‑140 (MQ=255)
                                                                      cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAgcg   >  1:1236726/1‑140 (MQ=255)
                                                                       tGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAgcgc  >  2:375699/1‑140 (MQ=255)
                                                                            |                                                                                                                                      
TTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCCTAGAACACCACAATTTCGCTCTCTCGGTAAGGCACTGTCACAGCATTCAACAGAATGTGACTTGCCAGATTTATTAGCGC  >  NC_000913/2067224‑2067434

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