Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,434,459 C→T 81.1% 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,4590CT81.1% 34.6 / 13.1 37intergenic (‑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 (3/4);  new base T (18/12);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTGGCGGCAAAAATCAGGCGGTTTATGGCG  >  NC_000913/1434362‑1434550
                                                                                                 |                                                                                           
gCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTg                                                                                           <  2:152515/100‑1 (MQ=25)
 cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                                                                          >  2:526497/1‑100 (MQ=18)
 cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                                                                          >  1:415025/1‑100 (MQ=25)
  tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTTTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTggg                                                                                         <  2:2487151/100‑1 (MQ=255)
  tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTggg                                                                                         >  2:4064437/1‑100 (MQ=255)
   gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgt                                                                                        >  2:2666036/1‑100 (MQ=17)
   gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATATTggg                                                                                         >  1:2066469/1‑99 (MQ=12)
    aaTATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgtg                                                                                       >  2:4149004/1‑100 (MQ=18)
     atatGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgtgt                                                                                      >  2:2801567/1‑100 (MQ=18)
      tatGATAAGGAGCGAAGTGATTATTAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGtt                                                                                     <  1:3427904/100‑1 (MQ=9)
       atGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCATCGAATTAGTAATGTGTTATATATGATATAGTTGAAAATTTTTATCTTGGATGTTc                                                                                    >  1:1049357/1‑100 (MQ=11)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                                                                   >  1:3063055/1‑100 (MQ=25)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                                                                   <  2:1673855/100‑1 (MQ=17)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                                                                   <  2:2990712/100‑1 (MQ=18)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                                                                   >  2:1835693/1‑100 (MQ=25)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTAATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                                                                   >  1:3001584/1‑100 (MQ=255)
         gATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCtt                                                                                  >  1:914781/1‑100 (MQ=18)
         gATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCtt                                                                                  >  2:1855468/1‑100 (MQ=17)
           tAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                                                                                >  2:2818447/1‑100 (MQ=18)
           tAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                                                                                >  1:2145668/1‑100 (MQ=18)
             aGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAg                                                                              >  2:3653035/1‑100 (MQ=18)
              ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGANTTAGNAANGNGTTANATATGATNTAGNTGACAATTTTTANCTTGGGTGTTCTTAAAGt                                                                             >  2:2845913/1‑100 (MQ=11)
              ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGt                                                                             <  1:3064758/100‑1 (MQ=18)
              ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGt                                                                             >  1:3083598/1‑100 (MQ=17)
               gAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGtt                                                                            <  1:806845/100‑1 (MQ=14)
               gAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGt                                                                             <  1:1461455/99‑1 (MQ=14)
                aGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTc                                                                           <  2:1362177/100‑1 (MQ=14)
                aGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTc                                                                           <  2:1361983/100‑1 (MQ=14)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                                                                          >  2:2852826/1‑100 (MQ=11)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                                                                          <  2:1733836/100‑1 (MQ=11)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                                                                          <  2:3029021/100‑1 (MQ=11)
                  cGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCGt                                                                         <  2:3036880/100‑1 (MQ=1)
                   gAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTa                                                                        <  2:3258717/100‑1 (MQ=38)
                                            tataGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGtata                                               >  2:765112/1‑100 (MQ=33)
                                                            tAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAgg                               <  2:4269150/100‑1 (MQ=32)
                                                              aTGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTg                             >  2:1737378/1‑100 (MQ=33)
                                                                                         tttttATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTGGCGGCAAAAATCAGGCGGTTTATGgcg  <  1:765110/100‑1 (MQ=38)
                                                                                                 |                                                                                           
GCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATACAGGAATGCTAACAGGTGGCGGCAAAAATCAGGCGGTTTATGGCG  >  NC_000913/1434362‑1434550

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