Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 168,421 G→C 6.7% intergenic (‑90/‑278) 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,4210GC 6.7% 289.3 / 8.4 90intergenic (‑90/‑278)map/rpsBmethionine aminopeptidase/30S ribosomal subunit protein S2
Reads supporting (aligned to +/- strand):  ref base G (41/43);  new base C (4/2);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTACACAGCG  >  minE/168352‑168490
                                                                     |                                                                     
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:904968/1‑70 (MQ=255)
         gACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGcc                                                              >  1:4478081/1‑70 (MQ=255)
         gACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGcc                                                              >  1:2783615/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:2871758/1‑71 (MQ=255)
              tatatatGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                        >  1:720679/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:4110454/44‑1 (MQ=25)
                            ccGTCAATGGTATCACACCCGGGCAAATTCAGAATGATTCTCaa                                                                     <  1:4227387/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:2821829/1‑41 (MQ=255)
                                 aaTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                                   >  1:508737/1‑41 (MQ=255)
                                  ttGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                                   <  1:1798214/39‑1 (MQ=255)
                                  aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                                   <  1:2970492/40‑1 (MQ=255)
                                  aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                                   <  1:2599008/40‑1 (MQ=255)
                                  aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                                   <  1:4263996/40‑1 (MQ=255)
                                  aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                                   <  1:1911466/40‑1 (MQ=255)
                                  aTGGTATCACACCCGGGCAAATTGAGAATCATTCTGAAtt                                                                   <  1:1064877/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:1529923/45‑1 (MQ=255)
                                        tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                        <  1:4003258/45‑1 (MQ=255)
                                        tCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGt                                                        <  1:4520036/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:663302/1‑50 (MQ=255)
                                             cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAg                                              >  1:1106982/1‑50 (MQ=255)
                                             cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAg                                              >  1:4088802/1‑50 (MQ=255)
                                             cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         >  1:151807/1‑71 (MQ=255)
                                             cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         >  1:4088360/1‑71 (MQ=255)
                                             cccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGACACTGAAGGTTTTCTATAATAGAAAATTc                         >  1:1265642/1‑71 (MQ=255)
                                              ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAgg                                             >  1:1463696/1‑50 (MQ=255)
                                              ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         <  1:2746858/70‑1 (MQ=255)
                                              ccGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         >  1:2634485/1‑70 (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:717955/1‑71 (MQ=255)
                                                        tGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATg              >  1:877752/1‑71 (MQ=255)
                                                         gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATg              <  1:1090251/70‑1 (MQ=255)
                                                         gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc             >  1:1765332/1‑71 (MQ=255)
                                                         gagaATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGc             >  1:1770857/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)
                                                                aTTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTataata                                 >  1:2442417/1‑44 (MQ=255)
                                                                 ttCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTacac      <  1:2843873/70‑1 (MQ=255)
                                                                  tCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         <  1:562172/50‑1 (MQ=255)
                                                                  tCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         <  1:1617571/50‑1 (MQ=255)
                                                                  tCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         <  1:648567/50‑1 (MQ=255)
                                                                  tCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTc                         <  1:3768892/50‑1 (MQ=255)
                                                                    tGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTACACAgcg  >  1:1151832/1‑71 (MQ=255)
                                                                     |                                                                     
GGTGATACCGACAATATATATGTAAGTGCCGTCAATGGTATCACACCCGGGCAAATTGAGAATCATTCTGAATTTCGCCAAACGTGCCACTGAAGGTTTTCTATAATAGAAAATTCGACGTCTGATGCTGTACACAGCG  >  minE/168352‑168490

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: