Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,065,979 A→T 6.1% intergenic (‑150/+17) ygfI ← / ← yggE predicted DNA‑binding transcriptional regulator/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,065,9790AT 6.1% 181.1 / 3.7 65intergenic (‑150/+17)ygfI/yggEpredicted DNA‑binding transcriptional regulator/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/26);  new base T (2/2);  total (37/28)
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.17e-01

TCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGG  >  W3110S.gb/3065913‑3066047
                                                                  |                                                                    
tCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                  <  1:1677033/71‑1 (MQ=255)
    ggTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                  >  1:825624/1‑67 (MQ=255)
      ttttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGt                                                                   <  1:450925/64‑1 (MQ=255)
      ttttACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTa                                                            <  1:7843969/71‑1 (MQ=255)
          acagaaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACaa                                                         <  1:6990938/70‑1 (MQ=255)
              aaaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAatat                                                     >  1:7698886/1‑70 (MQ=255)
               aaattacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAAtatta                                                   >  1:156363/1‑71 (MQ=255)
                  ttacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtg                                                <  1:7009463/71‑1 (MQ=255)
                   tacagaaaaAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATtgt                                               >  1:6126835/1‑71 (MQ=255)
                         aaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgctg                                           >  1:2637/1‑69 (MQ=255)
                         aaaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCa                                         <  1:3675236/71‑1 (MQ=255)
                          aaaGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAg                                        >  1:5978788/1‑71 (MQ=255)
                               aggCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtt                                   <  1:1879033/71‑1 (MQ=255)
                                 gCAATATCGGGTAAAGGCATTAGCCCGACGTATTCGTCggg                                                               >  1:1560559/1‑41 (MQ=25)
                                 gCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTc                                                                  <  1:787196/38‑1 (MQ=255)
                                   aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                >  1:2887482/1‑38 (MQ=255)
                                   aaTATCGGGTAAAGGCATTAGCCCGACGAATACGTCgg                                                                >  1:1052574/1‑38 (MQ=255)
                                    atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtttt                                 <  1:7057984/68‑1 (MQ=255)
                                    atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGc                              <  1:512349/71‑1 (MQ=255)
                                    atatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGc                              >  1:4442299/1‑71 (MQ=255)
                                     tatCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                             >  1:6797741/1‑71 (MQ=255)
                                      atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGtgc                                             <  1:1671475/54‑1 (MQ=255)
                                      atCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCgg                            >  1:5195344/1‑71 (MQ=255)
                                        cGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCggg                           >  1:7243138/1‑70 (MQ=255)
                                         gggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgtt                         <  1:3890399/71‑1 (MQ=255)
                                          ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                             <  1:5532412/66‑1 (MQ=255)
                                          ggTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgtt                         <  1:7095456/70‑1 (MQ=255)
                                               aGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGa                    >  1:3736102/1‑70 (MQ=255)
                                               aGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGAt                   >  1:6921970/1‑71 (MQ=255)
                                                 gCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgt                                      <  1:2932787/50‑1 (MQ=255)
                                                 gCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATcc                 <  1:2162935/71‑1 (MQ=255)
                                                 gCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATc                  >  1:7622053/1‑70 (MQ=255)
                                                   aTTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtgt                                    <  1:1914890/50‑1 (MQ=255)
                                                    ttAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGgtg                                     <  1:1040417/48‑1 (MQ=255)
                                                    ttAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                               <  1:4912364/54‑1 (MQ=255)
                                                     tAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAg             >  1:5083827/1‑71 (MQ=255)
                                                      aGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCgg                            >  1:6378654/1‑55 (MQ=255)
                                                       gCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAgg            <  1:3881366/70‑1 (MQ=255)
                                                        cccGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac               <  1:6333582/66‑1 (MQ=255)
                                                         ccGACGTATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGa                    <  1:2794715/60‑1 (MQ=39)
                                                          cgacgtATTCGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCg                             >  1:306314/1‑50 (MQ=38)
                                                            acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGgt                                      >  1:630953/1‑39 (MQ=255)
                                                            acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGgt                                      >  1:3934027/1‑39 (MQ=255)
                                                            acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGtttt                                 >  1:4081943/1‑44 (MQ=255)
                                                            acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTaa      >  1:1195218/1‑71 (MQ=255)
                                                            acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTaa      >  1:796768/1‑71 (MQ=255)
                                                            acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTaa      >  1:4834639/1‑71 (MQ=255)
                                                            acgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTaa      >  1:3664949/1‑71 (MQ=255)
                                                             cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                               >  1:1970609/1‑45 (MQ=255)
                                                             cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                               >  1:4711786/1‑45 (MQ=255)
                                                             cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAg                               >  1:1961341/1‑45 (MQ=255)
                                                             cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac               <  1:98553/61‑1 (MQ=255)
                                                             cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac               <  1:1448706/61‑1 (MQ=255)
                                                             cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCcac               <  1:6135580/61‑1 (MQ=255)
                                                             cgaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAAc     >  1:7191213/1‑71 (MQ=255)
                                                              gaATACGTCGGTCTACAAATATTATTGTGCTGCAGGTGtt                                   <  1:6935673/40‑1 (MQ=255)
                                                              gaATACGTCGGGCTACAAATATTATTGTGCTGCAGgtg                                     >  1:4330746/1‑38 (MQ=255)
                                                              gaATACGTCGGGCTACAAATATTATTGTGCTGCAGgtg                                     >  1:2589931/1‑38 (MQ=255)
                                                              gaATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCgg                            >  1:7717649/1‑47 (MQ=255)
                                                                aTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTa                                >  1:3887123/1‑41 (MQ=255)
                                                                aTACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTa                                >  1:5072666/1‑41 (MQ=255)
                                                                 tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGgttgt                       <  1:1919862/49‑1 (MQ=255)
                                                                 tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGtt          <  1:1418581/62‑1 (MQ=255)
                                                                 tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTgg  >  1:3790581/1‑70 (MQ=255)
                                                                 tACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGAGGGTTGTTGATCCACAGGTTCTAACTgg  >  1:2313190/1‑70 (MQ=255)
                                                                  |                                                                    
TCATGGTTTTACAGAAAATTACAGAAAAAGGAGGCAATATCGGGTAAAGGCATTAGCCCGACGAATACGTCGGGCTACAAATATTATTGTGCTGCAGGTGTTTTAGCGGGTTGTTGATCCACAGGTTCTAACTGG  >  W3110S.gb/3065913‑3066047

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