Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,434,459 C→T 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,4590CT93.9% 58.0 / ‑2.7 34intergenic (‑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 (9/22);  total (10/24)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATAC  >  NC_000913/1434362‑1434506
                                                                                                 |                                               
gCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTg                                               <  1:2965399/100‑1 (MQ=18)
 cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                              <  1:1153129/100‑1 (MQ=17)
 cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                              <  1:2563524/100‑1 (MQ=25)
 cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                              >  2:3242889/1‑100 (MQ=18)
 cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                              <  2:3240672/100‑1 (MQ=18)
 cTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTgg                                              <  1:2874433/100‑1 (MQ=18)
  tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                             <  2:1417018/100‑1 (MQ=25)
  tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                             <  2:1792243/100‑1 (MQ=17)
  tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTggg                                             <  2:2439805/100‑1 (MQ=25)
  tGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATATGTTATATATGATATAGTTGACAATTTTTATCTTggg                                             <  1:1866303/100‑1 (MQ=255)
   gAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTCTTATATAGGATATAGTTGACAATTTTTATCTTTGgt                                            >  2:3367249/1‑100 (MQ=11)
    aaTATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGgtg                                           <  2:2586077/100‑1 (MQ=17)
       atGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTc                                        <  1:2962089/100‑1 (MQ=17)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                       <  2:1674142/100‑1 (MQ=14)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                       <  2:446497/100‑1 (MQ=18)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGACTCGAATTAGGGATGTGTTATATATGAGATAGTTGACAATTTTTATCTTGGGTGTTCt                                       >  1:1322578/1‑100 (MQ=11)
        tGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATAGAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCt                                       >  1:2136389/1‑100 (MQ=255)
           tAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaa                                    >  2:636033/1‑100 (MQ=18)
            aaGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTaaa                                   <  1:2792085/100‑1 (MQ=17)
             aGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAg                                  <  2:2923938/100‑1 (MQ=17)
              ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGt                                 >  2:2196280/1‑100 (MQ=17)
              ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGt                                 >  2:957209/1‑100 (MQ=17)
              ggAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGt                                 <  1:35679/100‑1 (MQ=255)
               gAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGtt                                <  2:742810/100‑1 (MQ=17)
                aGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTc                               >  1:2686182/1‑100 (MQ=14)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                              <  2:1268513/100‑1 (MQ=14)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                              >  2:1221801/1‑100 (MQ=14)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                              <  1:265429/100‑1 (MQ=11)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                              <  1:1636006/100‑1 (MQ=11)
                 gCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCg                              <  1:1635685/100‑1 (MQ=11)
                  cGAAGTTATTATCAGTATGCTGTTCATATAGCCTCGAATTATTAATGTGTTATATATGATATAGTTGGCAATTTTTATCTTGGGTGTTCTTata                                   <  1:898815/94‑3 (MQ=11)
                  cGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCTTGGGTGTTCTTAAAGTTCGt                             <  1:1626958/100‑1 (MQ=1)
                                             ataGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATAc  >  1:3371491/1‑100 (MQ=32)
                                                           gttaTGTGTTATATATGAGACAGTTGACGATTTTTATTGTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGt      <  1:1430592/79‑1 (MQ=1)
                                                                                                 |                                               
GCTGAATATGATAAGGAGCGAAGTGATTATCAGTATGCTGTTCATATAGCCTCGAATTAGTAATGTGTTATATATGATATAGTTGACAATTTTTATCCTGGGTGTTCTTAAAGTTCGTAGATAAACATTGTCGTTTCAGGTATAC  >  NC_000913/1434362‑1434506

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: