Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,218,223 T→A 7.4% intergenic (+11/‑102) ymgB → / → ymgC hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,218,2230TA 7.4% 227.0 / 8.9 81intergenic (+11/‑102)ymgB/ymgChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/38);  new base A (3/3);  total (40/41)
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.52e-01

CTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAT  >  W3110S.gb/1218154‑1218293
                                                                     |                                                                      
cTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattatt                                                                       <  1:4829954/71‑1 (MQ=255)
cTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattatt                                                                       <  1:1637612/71‑1 (MQ=255)
   aTGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAg                                                                    <  1:3309601/71‑1 (MQ=255)
     gtgATTAGAAAAACACTGGAAATCGTGTTGAGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                  >  1:6170359/1‑71 (MQ=255)
       gATTAGAAAAACACTTGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattat                                                                        <  1:5499315/63‑1 (MQ=255)
       gATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattat                                                                        <  1:2605406/63‑1 (MQ=255)
       gATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattat                                                                        <  1:6259390/63‑1 (MQ=255)
        aTTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAattatt                                                                       <  1:6936826/63‑1 (MQ=255)
           aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                            >  1:4260427/1‑71 (MQ=255)
           aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                            >  1:6312619/1‑71 (MQ=255)
            gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                            <  1:3456177/70‑1 (MQ=255)
            gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                            <  1:6869511/70‑1 (MQ=255)
            gAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTATACTTAATAATTAGAAGt                                                             <  1:5904878/69‑1 (MQ=25)
             aaaaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAagtttaagttta                                                          <  1:5905985/71‑1 (MQ=255)
               aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                          >  1:7762676/1‑69 (MQ=255)
               aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                          >  1:4043939/1‑69 (MQ=255)
               aaaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                          >  1:5623980/1‑69 (MQ=255)
                aaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAg                                                       <  1:72822/71‑1 (MQ=255)
                aaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAg                                                       <  1:1436108/71‑1 (MQ=255)
                 acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                  >  1:7841438/1‑59 (MQ=255)
                 acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                  >  1:699095/1‑59 (MQ=255)
                 acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                  >  1:975691/1‑59 (MQ=255)
                 acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGa                                                      <  1:5537798/71‑1 (MQ=255)
                  cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGttt                                                           >  1:7819930/1‑65 (MQ=255)
                  cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGaa                                                     <  1:227899/71‑1 (MQ=255)
                  cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTATACTTAATAATTAGAAGttac                                                          >  1:4536046/1‑64 (MQ=21)
                  cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTATACTTAATAATTAGAAGttac                                                          >  1:5316384/1‑64 (MQ=21)
                     tGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                   >  1:3762897/1‑70 (MQ=255)
                     tGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAAc                                                  >  1:339713/1‑71 (MQ=255)
                     tGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAAc                                                  >  1:3684896/1‑71 (MQ=255)
                       gAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                            <  1:6151840/59‑1 (MQ=255)
                       gAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtc                                                <  1:2550131/71‑1 (MQ=255)
                         aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                          <  1:2153894/59‑1 (MQ=255)
                         aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                          <  1:5647634/59‑1 (MQ=255)
                         aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                          <  1:2918116/59‑1 (MQ=255)
                         aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcat                                              <  1:4396429/71‑1 (MQ=255)
                         aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcat                                              <  1:7699811/71‑1 (MQ=255)
                             gtgtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcattca                                           >  1:2820014/1‑70 (MQ=255)
                              tgtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcattcat                                          >  1:4402613/1‑70 (MQ=255)
                               gtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGttt                                                           <  1:4651752/52‑1 (MQ=255)
                               gtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGttt                                                           <  1:6688921/52‑1 (MQ=255)
                                ttGCGATACACAGCTGATTATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCt                                        <  1:7543406/70‑1 (MQ=255)
                                ttGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                             <  1:2728834/49‑1 (MQ=255)
                                ttGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCt                                        >  1:615799/1‑70 (MQ=255)
                                  gCGATACACAGCTGATGATATGTAACTTCTAAttatta                                                                      >  1:2598859/1‑38 (MQ=255)
                                    gATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttt                                   <  1:5623838/71‑1 (MQ=255)
                                     aTACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAtttt                                  <  1:4666071/71‑1 (MQ=255)
                                     aTACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAtttt                                  <  1:1369947/71‑1 (MQ=255)
                                      tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtca                                               <  1:4278364/57‑1 (MQ=255)
                                      tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                                 <  1:2331540/71‑1 (MQ=255)
                                      tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                                 <  1:6856924/71‑1 (MQ=255)
                                      tACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTAttttt                                 <  1:7722735/71‑1 (MQ=255)
                                       acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                             <  1:3529732/42‑1 (MQ=255)
                                       acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                             <  1:5071444/42‑1 (MQ=255)
                                       acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtg                                >  1:2030201/1‑71 (MQ=255)
                                        cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                              >  1:2177517/1‑40 (MQ=255)
                                        cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                              >  1:1589774/1‑40 (MQ=255)
                                        cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                              >  1:6351234/1‑40 (MQ=255)
                                        cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                              >  1:3927517/1‑40 (MQ=255)
                                        cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgt                               >  1:5354050/1‑71 (MQ=255)
                                         acaGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTtata                                                        >  1:1951783/1‑45 (MQ=255)
                                         acaGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgtc                              <  1:2697953/71‑1 (MQ=255)
                                           aGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACtcattc                                            >  1:1633118/1‑55 (MQ=255)
                                            gCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACt                                                 <  1:2991533/49‑1 (MQ=255)
                                            gCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgtctg                            >  1:6986766/1‑70 (MQ=255)
                                                atgatATGTAACTTCTAATTATTAAGTATAAGTTta                                                          <  1:6489890/36‑1 (MQ=255)
                                                atgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt                       >  1:5702830/1‑71 (MQ=255)
                                                atgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt                       <  1:4561354/71‑1 (MQ=255)
                                                 tgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt                       <  1:3016612/70‑1 (MQ=255)
                                                    tatGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgtctgt                           >  1:4541322/1‑63 (MQ=255)
                                                       taaaCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGAc                 >  1:2404236/3‑70 (MQ=25)
                                                       gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                   >  1:1478641/1‑36 (MQ=255)
                                                       gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                   >  1:3037038/1‑36 (MQ=255)
                                                       gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                   >  1:3027862/1‑36 (MQ=255)
                                                       gTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCtt                                       <  1:4420718/48‑1 (MQ=255)
                                                          aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                      >  1:1259737/1‑46 (MQ=255)
                                                          aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                      >  1:5295045/1‑46 (MQ=255)
                                                          aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                      >  1:5557294/1‑46 (MQ=255)
                                                          aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata             <  1:5555411/71‑1 (MQ=255)
                                                          aCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata             <  1:3605096/71‑1 (MQ=37)
                                                          aCTTATACTTAATAATTAGAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata             <  1:3614554/71‑1 (MQ=37)
                                                           cTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCttt                     >  1:4967469/1‑62 (MQ=255)
                                                           cTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACtata             <  1:6298013/70‑1 (MQ=255)
                                                           cTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATac            <  1:434514/71‑1 (MQ=255)
                                                              cTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATaca           <  1:1429108/69‑1 (MQ=255)
                                                              cTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATaca           >  1:1770393/1‑69 (MQ=255)
                                                               tAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGc                        >  1:570792/1‑55 (MQ=255)
                                                               tAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGc                        >  1:5041284/1‑55 (MQ=255)
                                                                 attattAACTATAAGTTTATAGAAAACTCATTCATCtt                                       <  1:2362052/38‑1 (MQ=255)
                                                                   tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                    <  1:7076874/39‑1 (MQ=255)
                                                                   tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                    <  1:4180937/39‑1 (MQ=255)
                                                                   tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                    <  1:3629268/39‑1 (MQ=255)
                                                                   tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                    <  1:3408556/39‑1 (MQ=255)
                                                                   tattaAGTATAAGCTTATAGAAAACTCATTCATCTTAtt                                    <  1:705084/39‑1 (MQ=255)
                                                                     ttaAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtgt                               >  1:6396115/1‑42 (MQ=255)
                                                                     ttaAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAt  <  1:1220540/71‑1 (MQ=255)
                                                                     |                                                                      
CTGATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAT  >  W3110S.gb/1218154‑1218293

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: