Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,218,227 G→T 6.5% intergenic (+15/‑98) 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,2270GT 6.5% 282.7 / 2.6 89intergenic (+15/‑98)ymgB/ymgChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/44);  new base T (3/3);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTA  >  W3110S.gb/1218157‑1218297
                                                                      |                                                                      
aTGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAg                                                                        <  1:3309601/71‑1 (MQ=255)
  gtgATTAGAAAAACACTGGAAATCGTGTTGAGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                      >  1:6170359/1‑71 (MQ=255)
        aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                                >  1:6312619/1‑71 (MQ=255)
        aGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGtt                                                                >  1:4260427/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:1436108/71‑1 (MQ=255)
             aaCACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAg                                                           <  1:72822/71‑1 (MQ=255)
              acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                      >  1:975691/1‑59 (MQ=255)
              acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                      >  1:7841438/1‑59 (MQ=255)
              acacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGta                                                                      >  1:699095/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:5316384/1‑64 (MQ=21)
               cacTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTATACTTAATAATTAGAAGttac                                                              >  1:4536046/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:5647634/59‑1 (MQ=255)
                      aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                              <  1:2918116/59‑1 (MQ=255)
                      aaTCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTta                                                              <  1:2153894/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:6688921/52‑1 (MQ=255)
                            gtTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGttt                                                               <  1:4651752/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)
                                 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:5071444/42‑1 (MQ=255)
                                    acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGt                                                                 <  1:3529732/42‑1 (MQ=255)
                                    acacAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTtg                                    >  1:2030201/1‑71 (MQ=255)
                                     cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                  >  1:6351234/1‑40 (MQ=255)
                                     cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                  >  1:2177517/1‑40 (MQ=255)
                                     cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                  >  1:3927517/1‑40 (MQ=255)
                                     cacaGCTGATGATATGTAACTTCTAATTATTAAGTATAAg                                                                  >  1:1589774/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:4561354/71‑1 (MQ=255)
                                             atgatATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCt                           >  1:5702830/1‑71 (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:3027862/1‑36 (MQ=255)
                                                    gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                       >  1:1478641/1‑36 (MQ=255)
                                                    gTAACTTCTAATTATTAAGTATAAGTTTATAGaaaa                                                       >  1:3037038/1‑36 (MQ=255)
                                                    gTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCtt                                           <  1:4420718/48‑1 (MQ=255)
                                                       aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                          >  1:5557294/1‑46 (MQ=255)
                                                       aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                          >  1:1259737/1‑46 (MQ=255)
                                                       aCTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTa                                          >  1:5295045/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:3629268/39‑1 (MQ=255)
                                                                tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                        <  1:3408556/39‑1 (MQ=255)
                                                                tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                        <  1:7076874/39‑1 (MQ=255)
                                                                tattaAGTATAAGTTTATAGAAAACTCATTCATCTTAtt                                        <  1:4180937/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)
                                                                     aGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAtttt   <  1:4161463/71‑1 (MQ=255)
                                                                     aGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAAtttt   <  1:6100145/71‑1 (MQ=255)
                                                                      gTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTa  <  1:4528763/71‑1 (MQ=255)
                                                                      gTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTa  <  1:5760332/71‑1 (MQ=255)
                                                                      gTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTa  <  1:4595333/71‑1 (MQ=255)
                                                                      |                                                                      
ATGTGATTAGAAAAACACTGGAAATCGTGTTGCGATACACAGCTGATGATATGTAACTTCTAATTATTAAGTATAAGTTTATAGAAAACTCATTCATCTTATTTTTGTCTGTCGCTTTAGACTATACACAGCATAATTTTA  >  W3110S.gb/1218157‑1218297

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: