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,4590CT86.4% 29.0 / 1.9 22intergenic (‑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/2);  new base T (10/9);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTGGCGGCAAAAATCAGGC  >  NC_000913/1434363‑1434539
                                                                                                |                                                                                
cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                                                               <  1:2400742/100‑1 (MQ=25)
 tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                                                              >  1:258340/1‑100 (MQ=25)
 tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTAATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                                                              <  2:2117630/100‑1 (MQ=12)
  gAATATGATAAGTAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                                                             >  1:1686127/1‑100 (MQ=255)
  gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                                                             >  1:1223937/1‑100 (MQ=18)
  gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                                                             >  2:1274855/1‑100 (MQ=18)
   aaTATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgtg                                                                            <  2:356903/100‑1 (MQ=14)
   aaTATGATAAGGAGCGAAGTGATTATCAGTATGATGTTCATATAGCCTCGAATTAGTAATGTGTTATATATTATATAGTTGACAATTGTTATCTTGGgtg                                                                            >  2:1912985/1‑100 (MQ=11)
     tatGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGtt                                                                          <  2:1605683/100‑1 (MQ=14)
      ntGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGtt                                                                          >  1:1045563/2‑99 (MQ=25)
      atGATAAGGAGCGAAGTTATTATCAGTATGCTGTTCATATAGCCTGGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                                                             <  2:2039572/96‑1 (MQ=255)
      atGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTc                                                                         <  2:334255/100‑1 (MQ=18)
       tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                                                        >  2:1692066/1‑100 (MQ=18)
       tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                                                        <  1:492981/100‑1 (MQ=14)
          tAAGGAGGCAAGTGATTATCAGTCTGCTGCATATATAGCCACGAATTAGTAAAGTGTTATATATGACATCGTTGACAATTTTTTTTTTGGGTGTCCTTaa                                                                     >  2:1941669/1‑100 (MQ=255)
          tAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                                                                     >  2:1061996/1‑100 (MQ=18)
            aGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAg                                                                   <  1:551500/100‑1 (MQ=18)
             ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGt                                                                  <  2:930395/100‑1 (MQ=14)
                 cGAAGTGATTATCAGTATGCTGTTCATATAGCATCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCGt                                                              >  1:523113/1‑100 (MQ=2)
                                tATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCg                                               >  1:1284561/1‑100 (MQ=33)
                                                                  ttATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTggcggc             <  2:2023330/100‑1 (MQ=35)
                                                                             ataGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTGGCGGCAAAAATCAGGc  <  2:1284553/100‑1 (MQ=38)
                                                                                                |                                                                                
CTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTGGCGGCAAAAATCAGGC  >  NC_000913/1434363‑1434539

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