Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,434,459 C→T 100% intergenic (‑235/+499) ynaE ← / ← ttcC cold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,434,4590CT88.9% 24.6 / ‑0.8 18intergenic (‑235/+499)ynaE/ttcCcold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (7/9);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAA  >  NC_000913/1434363‑1434511
                                                                                                |                                                    
cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                                   >  1:2346003/1‑100 (MQ=18)
cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                                   >  2:1083842/1‑100 (MQ=18)
 tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                                  <  2:786892/100‑1 (MQ=18)
 tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                                  <  2:354508/100‑1 (MQ=17)
 tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                                  >  1:1250137/1‑100 (MQ=25)
  gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                                 <  1:1710800/100‑1 (MQ=14)
   aaTATGATGAGGAGCGAAGTGATTATCAATATGCTTTTTATATATGTTGAAATTAGTAATGTGTTATATATGATATCGTTGACAATTTTTATCTTGGgtg                                                <  2:2640827/100‑1 (MQ=11)
   aaTATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgtg                                                <  2:1781470/100‑1 (MQ=18)
       tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                            <  1:1554053/100‑1 (MQ=14)
        gATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCtt                                           >  2:1733225/1‑100 (MQ=25)
         aTAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTa                                          <  1:1010004/100‑1 (MQ=17)
            aGGAGCGAAGTGCTTATCAGTATGTTGTTCATATAGCCTTGCATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAg                                       <  1:2895885/100‑1 (MQ=11)
            aGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAg                                       >  1:1621758/1‑100 (MQ=25)
             ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGt                                      >  2:2217184/1‑100 (MQ=17)
                gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                                   >  1:347121/1‑100 (MQ=14)
                gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                                   <  1:2884886/100‑1 (MQ=11)
                                     tGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGATAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCa              <  1:1208190/100‑1 (MQ=17)
                                                 cTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGaa  >  1:1051846/1‑100 (MQ=33)
                                                                                                |                                                    
CTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAA  >  NC_000913/1434363‑1434511

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