Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,872,074 T→C 16.0% intergenic (+122/‑25) yeaH → / → yeaI conserved hypothetical protein/predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,872,0740TC16.0% 56.3 / 5.7 25intergenic (+122/‑25)yeaH/yeaIconserved hypothetical protein/predicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (12/9);  new base C (2/2);  total (14/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 = 9.85e-01

CATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTT  >  W3110S.gb/1872009‑1872136
                                                                 |                                                              
cATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                                              >  1:1486417/1‑68 (MQ=255)
      ttaCATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAAc                                                      <  1:891682/70‑1 (MQ=255)
           tacatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                                >  1:478952/1‑71 (MQ=255)
             catTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTc                                              >  1:27549/1‑71 (MQ=255)
                  gCTTACAATCGCTTTAAATATGACAGGATAACCTTTACATAATTTAGTTCCAGaaaa                                                       >  1:998209/1‑57 (MQ=255)
                      aCAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                                       <  1:189439/69‑1 (MQ=255)
                       cAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGa                                    >  1:1430679/1‑71 (MQ=255)
                         aTCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAtt                                  <  1:711043/71‑1 (MQ=255)
                          tCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAtt                                               <  1:371462/57‑1 (MQ=255)
                              tttAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCa                                <  1:678395/68‑1 (MQ=255)
                                tAAATATGACAGCATAACCTTTACATAATTTAGTTc                                                              <  1:1608988/36‑1 (MQ=255)
                                 aaaTATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCaa                           >  1:1658637/1‑70 (MQ=255)
                                  aaTATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATc                                                  >  1:582417/1‑46 (MQ=255)
                                  aaTATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATc                                                  >  1:91733/1‑46 (MQ=255)
                                   atatGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                             <  1:435271/66‑1 (MQ=255)
                                     atGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                                   <  1:1626930/42‑1 (MQ=255)
                                      tGACAGCATAACCTTTACATAATTTAGCTCCAGAAAACAAt                                                   <  1:1608264/41‑1 (MQ=255)
                                        aCAGCATAACCTTTACATAATTTAGCTCCAGAAAACAATCa                                                 >  1:1750610/1‑41 (MQ=255)
                                              tAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCaa                           <  1:959006/57‑1 (MQ=255)
                                              tAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATg             >  1:608581/1‑71 (MQ=255)
                                              tAACCTTTACATAATTTAGCTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATg             <  1:285116/71‑1 (MQ=255)
                                                aCCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCa               >  1:374978/1‑67 (MQ=255)
                                                      aCATAATTTAGCTCCAGAAAACAATCATTCGGAAAAATGa                                    >  1:573715/1‑40 (MQ=255)
                                                         tAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCtt  >  1:1694647/1‑71 (MQ=255)
                                                         tAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCtt  >  1:1186469/1‑71 (MQ=255)
                                                                 |                                                              
CATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATTCTT  >  W3110S.gb/1872009‑1872136

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