Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 168,409 G→C 7.4% intergenic (‑78/‑290) 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,4090GC 7.4% 258.1 / 8.8 81intergenic (‑78/‑290)map/rpsBmethionine aminopeptidase/30S ribosomal subunit protein S2
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base C (4/2);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

TCCATCAGCGTCGGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGC  >  minE/168340‑168479
                                                                     |                                                                      
tCCATCAGCGTCGGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTg                                                                        <  1:816934/70‑1 (MQ=255)
       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:2769455/70‑1 (MQ=255)
                   ccGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCg                                                     <  1:156460/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)
                     gACAATATATATGTAAGTGCAGTCAATGGTATCACACCCGGGCAAATTGAGaa                                                                    >  1:251794/1‑53 (MQ=255)
                      aCAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCaa                                                 <  1:3465452/71‑1 (MQ=255)
                      aCAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCa                                                  >  1:779564/1‑70 (MQ=255)
                        aaTATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAAc                                               <  1:894310/71‑1 (MQ=255)
                         atatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAAt                                                                   <  1:3227966/50‑1 (MQ=255)
                          tatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACg                                              <  1:896601/70‑1 (MQ=255)
                          tatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                             >  1:2871758/1‑71 (MQ=255)
                          tatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                             >  1:720679/1‑71 (MQ=255)
                            tatatGTAAGTGCCGTCAATGGTATCACAGCCGGGCAAATTGAGAATCATTCTGAATTTCGCCaa                                                 <  1:1739592/65‑1 (MQ=255)
                                tGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTga                                                                       >  1:4444075/1‑39 (MQ=255)
                                tGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTga                                                                       >  1:326538/1‑39 (MQ=255)
                                tGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTga                                                                       >  1:4218977/1‑39 (MQ=255)
                                tGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGaa                                                                    >  1:2894666/1‑42 (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:2095226/71‑1 (MQ=255)
                                        ccGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGttt                               <  1:4193212/71‑1 (MQ=255)
                                        ccGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGttt                               <  1:1051786/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:3001849/1‑44 (MQ=21)
                                          gTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCAAtt                                                        >  1:2406428/1‑44 (MQ=21)
                                          gTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCAAtt                                                        >  1:1898710/1‑44 (MQ=21)
                                            cAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCta                           >  1:1510324/1‑71 (MQ=255)
                                             aaTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        >  1:508737/1‑41 (MQ=255)
                                             aaTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        >  1:2821829/1‑41 (MQ=255)
                                             aaTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        >  1:4674065/1‑41 (MQ=255)
                                              ttGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        <  1:1798214/39‑1 (MQ=255)
                                              aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        <  1:4263996/40‑1 (MQ=255)
                                              aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        <  1:1064877/40‑1 (MQ=255)
                                              aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        <  1:2970492/40‑1 (MQ=255)
                                              aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        <  1:1911466/40‑1 (MQ=255)
                                              aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                        <  1:2599008/40‑1 (MQ=255)
                                                ggTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataa                        >  1:180078/1‑70 (MQ=255)
                                                 gTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataata                      <  1:3720988/71‑1 (MQ=255)
                                                 gTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataata                      <  1:1367233/71‑1 (MQ=255)
                                                    tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                             <  1:4003258/45‑1 (MQ=255)
                                                    tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                             <  1:4520036/45‑1 (MQ=255)
                                                    tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                             <  1:1529923/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:2319643/1‑58 (MQ=255)
                                                      acacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGtttt                              >  1:2899849/1‑58 (MQ=255)
                                                      acacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGaaaa                 >  1:4226054/1‑71 (MQ=255)
                                                       cacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGa                                     <  1:3354600/50‑1 (MQ=255)
                                                       cacCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGa                                     <  1:998523/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:1765057/70‑1 (MQ=255)
                                                          ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc              <  1:2746858/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:717955/1‑71 (MQ=255)
                                                                    tGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATg   >  1:877752/1‑71 (MQ=255)
                                                                     gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATg   <  1:1090251/70‑1 (MQ=255)
                                                                     gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc  >  1:3903866/1‑71 (MQ=255)
                                                                     gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc  >  1:1770857/1‑71 (MQ=255)
                                                                     gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc  >  1:1765332/1‑71 (MQ=255)
                                                                     |                                                                      
TCCATCAGCGTCGGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGC  >  minE/168340‑168479

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: