Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 168,415 C→G 6.1% intergenic (‑84/‑284) map ← / → rpsB methionine aminopeptidase/30S ribosomal subunit protein S2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE168,4150CG 6.1% 273.0 / 5.4 86intergenic (‑84/‑284)map/rpsBmethionine aminopeptidase/30S ribosomal subunit protein S2
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base G (3/2);  total (44/42)
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.88e-01

GCGTCGGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTACA  >  minE/168347‑168485
                                                                    |                                                                      
gCGTCGGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCAt                                                                      <  1:3844444/71‑1 (MQ=255)
gCGTCGGTGATACCGACAATATATATGTAAGTGCCGGCAATGGTATCACACCCGGGCAAATTGAGAATCAt                                                                      <  1:1972023/71‑1 (MQ=255)
   tCGGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCt                                                                   >  1:2036739/1‑71 (MQ=255)
     ggTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTg                                                                  >  1:1431332/1‑70 (MQ=255)
       tGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAt                                                               <  1:1984200/71‑1 (MQ=255)
          tACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTc                                                            >  1:3177109/1‑71 (MQ=255)
           aCCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCg                                                           >  1:4222670/1‑71 (MQ=255)
            ccGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCg                                                           <  1:156460/70‑1 (MQ=255)
            ccGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCg                                                           <  1:2769455/70‑1 (MQ=255)
              gACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGcc                                                         >  1:2783615/1‑70 (MQ=255)
              gACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGcc                                                         >  1:904968/1‑70 (MQ=255)
              gACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGcc                                                         >  1:4478081/1‑70 (MQ=255)
              gACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCa                                                        <  1:844703/71‑1 (MQ=255)
               aCAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCaa                                                       <  1:3465452/71‑1 (MQ=255)
               aCAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCa                                                        >  1:779564/1‑70 (MQ=255)
                 aaTATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAAc                                                     <  1:894310/71‑1 (MQ=255)
                   tatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACg                                                    <  1:896601/70‑1 (MQ=255)
                   tatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                   >  1:720679/1‑71 (MQ=255)
                   tatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                   >  1:2871758/1‑71 (MQ=255)
                     tatatGTAAGTGCCGTCAATGGTATCACAGCCGGGCAAATTGAGAATCATTCTGAATTTCGCCaa                                                       <  1:1739592/65‑1 (MQ=255)
                         tGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:579000/54‑1 (MQ=255)
                         tGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:134023/54‑1 (MQ=255)
                          gTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:3399784/53‑1 (MQ=255)
                              gTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAgg                                        <  1:1592850/71‑1 (MQ=255)
                                gCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGtt                                      >  1:1415737/1‑71 (MQ=255)
                                 ccGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGttt                                     <  1:1051786/71‑1 (MQ=255)
                                 ccGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGttt                                     <  1:4193212/71‑1 (MQ=255)
                                 ccGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGttt                                     <  1:2095226/71‑1 (MQ=255)
                                 ccGTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCaa                                                                <  1:4227387/44‑1 (MQ=25)
                                 ccGTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCaa                                                                <  1:4110454/44‑1 (MQ=25)
                                   gTCAATGGTATCACACCCGGGCAAATTCAGAATTATTCTCAAtt                                                              >  1:4169795/1‑44 (MQ=25)
                                   gTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCAAtt                                                              >  1:1898710/1‑44 (MQ=21)
                                   gTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCAAtt                                                              >  1:3001849/1‑44 (MQ=21)
                                   gTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCAAtt                                                              >  1:2406428/1‑44 (MQ=21)
                                     cAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCta                                 >  1:1510324/1‑71 (MQ=255)
                                      aaTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              >  1:4674065/1‑41 (MQ=255)
                                      aaTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              >  1:508737/1‑41 (MQ=255)
                                      aaTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              >  1:2821829/1‑41 (MQ=255)
                                       ttGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:1798214/39‑1 (MQ=255)
                                       aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:2970492/40‑1 (MQ=255)
                                       aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:1911466/40‑1 (MQ=255)
                                       aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:2599008/40‑1 (MQ=255)
                                       aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:4263996/40‑1 (MQ=255)
                                       aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                              <  1:1064877/40‑1 (MQ=255)
                                         ggTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataa                              >  1:180078/1‑70 (MQ=255)
                                          gTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataata                            <  1:1367233/71‑1 (MQ=255)
                                          gTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataata                            <  1:3720988/71‑1 (MQ=255)
                                             tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                   <  1:4520036/45‑1 (MQ=255)
                                             tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                   <  1:1529923/45‑1 (MQ=255)
                                             tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                   <  1:4003258/45‑1 (MQ=255)
                                             tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataat                             <  1:1812842/67‑1 (MQ=255)
                                             tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGa                          >  1:4762029/1‑70 (MQ=255)
                                               acacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGtttt                                    >  1:41808/1‑58 (MQ=255)
                                               acacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGtttt                                    >  1:2899849/1‑58 (MQ=255)
                                               acacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGtttt                                    >  1:2319643/1‑58 (MQ=255)
                                               acacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGaaaa                       >  1:4226054/1‑71 (MQ=255)
                                                cacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGa                                           <  1:998523/50‑1 (MQ=255)
                                                cacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGa                                           <  1:3354600/50‑1 (MQ=255)
                                                  cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAg                                         >  1:4088802/1‑50 (MQ=255)
                                                  cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAg                                         >  1:1106982/1‑50 (MQ=255)
                                                  cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAg                                         >  1:663302/1‑50 (MQ=255)
                                                  cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                    >  1:4088360/1‑71 (MQ=255)
                                                  cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                    >  1:151807/1‑71 (MQ=255)
                                                  cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGACACTGAAGGTTTTCTATAATAGAAAATTc                    >  1:1265642/1‑71 (MQ=255)
                                                   ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAgg                                        >  1:1463696/1‑50 (MQ=255)
                                                   ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                    >  1:2634485/1‑70 (MQ=255)
                                                   ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                    <  1:2746858/70‑1 (MQ=255)
                                                   ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                    <  1:1765057/70‑1 (MQ=255)
                                                    cGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGa                  >  1:704739/1‑71 (MQ=255)
                                                     gggCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGAc                 <  1:1609067/71‑1 (MQ=255)
                                                      ggCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGa                                           <  1:168072/44‑1 (MQ=255)
                                                        cAAATTGGGAATCATTCTGAATTTCGCCAAACGTGc                                                 >  1:287417/1‑36 (MQ=255)
                                                         aaaTTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCt             >  1:3995589/1‑71 (MQ=255)
                                                          aaTTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCt             >  1:4383859/1‑70 (MQ=255)
                                                           aTTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGa                  >  1:430261/1‑64 (MQ=255)
                                                           aTTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGa           >  1:1580488/1‑71 (MQ=255)
                                                             tGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATg         >  1:877752/1‑71 (MQ=255)
                                                             tGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATg         >  1:717955/1‑71 (MQ=255)
                                                              gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATg         <  1:1090251/70‑1 (MQ=255)
                                                              gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc        >  1:1770857/1‑71 (MQ=255)
                                                              gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc        >  1:1765332/1‑71 (MQ=255)
                                                              gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc        >  1:3903866/1‑71 (MQ=255)
                                                                 aaTCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGt     <  1:2346325/71‑1 (MQ=255)
                                                                  aTCATTCTGAATTTCGCCAAACGTGCCACTGAATGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGt     <  1:3533677/70‑1 (MQ=255)
                                                                  aTCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTa    <  1:4528602/71‑1 (MQ=255)
                                                                  aTCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTa    <  1:4320886/71‑1 (MQ=255)
                                                                    cATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTaca  >  1:1050501/1‑71 (MQ=255)
                                                                    |                                                                      
GCGTCGGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTACA  >  minE/168347‑168485

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