Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,359,8760CT31.0% 172.9 / 99.6 116noncoding (1219/2905 nt)rrlA23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (38/42);  new base T (30/6);  total (68/48)
Fisher's exact test for biased strand distribution p-value = 2.50e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 116 total aligned reads displayed.

CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359816‑2359945
                                                            |                                                                     
ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                              >  1:1907301/1‑70 (MQ=255)
    tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:2968886/1‑71 (MQ=255)
    tttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:2280583/1‑71 (MQ=255)
     ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                                         <  1:1967007/70‑1 (MQ=35)
     ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:3150403/1‑70 (MQ=255)
      tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        >  1:2733858/1‑70 (MQ=35)
      tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        >  1:2470773/1‑70 (MQ=35)
      tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCt                                                        >  1:1939222/1‑70 (MQ=35)
      tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                                        >  1:3037095/1‑70 (MQ=255)
      tATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCACCATACCTCACAGTACACCTTCACAGGCt                                                        >  1:1969175/1‑70 (MQ=11)
       aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                                       <  1:332117/70‑1 (MQ=35)
       aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      >  1:151163/1‑71 (MQ=34)
       aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                                       <  1:1649295/70‑1 (MQ=255)
        tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTa                                                      <  1:3181243/70‑1 (MQ=35)
          gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                                   <  1:539391/71‑1 (MQ=34)
            tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                     >  1:2370869/1‑67 (MQ=35)
            tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 >  1:335653/1‑71 (MQ=34)
            tACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 >  1:1868856/1‑71 (MQ=34)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 <  1:1399411/70‑1 (MQ=35)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 <  1:1121859/70‑1 (MQ=35)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 <  1:1204994/70‑1 (MQ=35)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 <  1:1931644/70‑1 (MQ=35)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGaa                                                 <  1:2298462/70‑1 (MQ=35)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 <  1:383993/70‑1 (MQ=255)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 <  1:428843/70‑1 (MQ=255)
             aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 <  1:1242278/70‑1 (MQ=255)
              cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               <  1:1911558/71‑1 (MQ=34)
              cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                               <  1:1684853/71‑1 (MQ=34)
                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCt                                             <  1:1304075/70‑1 (MQ=35)
                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            >  1:147692/1‑71 (MQ=34)
                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            >  1:145065/1‑71 (MQ=34)
                 aTGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                            >  1:57369/1‑71 (MQ=255)
                  tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                           <  1:691510/71‑1 (MQ=34)
                  tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGATc                                            <  1:1877005/70‑1 (MQ=25)
                     cAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                            >  1:758414/1‑67 (MQ=35)
                      aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                                       <  1:714709/55‑1 (MQ=35)
                      aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCtt                                                       <  1:3133983/55‑1 (MQ=35)
                      aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                        <  1:1618436/70‑1 (MQ=255)
                      aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCt                                        <  1:605396/70‑1 (MQ=255)
                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:3132319/1‑46 (MQ=32)
                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                                     >  1:2322363/1‑55 (MQ=35)
                        cATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:941583/1‑71 (MQ=34)
                        cATTCGCACTTCTGATACCTCCAGCATGCCTACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:1792660/1‑70 (MQ=255)
                        cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                                     >  1:1694804/1‑71 (MQ=255)
                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:200207/1‑45 (MQ=32)
                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:3088941/1‑45 (MQ=32)
                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:1437037/1‑45 (MQ=32)
                         aTTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                              >  1:1360791/1‑45 (MQ=32)
                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:2928787/1‑58 (MQ=11)
                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:3080901/1‑58 (MQ=11)
                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:2463169/1‑58 (MQ=11)
                         aTTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                                 >  1:731266/1‑58 (MQ=11)
                           tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1237256/1‑71 (MQ=255)
                           tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:23665/1‑71 (MQ=255)
                           tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:822661/1‑71 (MQ=255)
                           tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                   >  1:853895/1‑70 (MQ=255)
                           tCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCca                                   >  1:2073456/1‑70 (MQ=255)
                            cGCACTTCTGATACCTCCAGCATACCTCACAGTACCCCTTCACAGGc                                                         >  1:555927/1‑47 (MQ=11)
                            cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:1869252/1‑47 (MQ=11)
                            cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGc                                                         >  1:839819/1‑47 (MQ=11)
                            cGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:947785/1‑70 (MQ=255)
                             gCACTTCTGATACCTCGAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:1265470/70‑1 (MQ=255)
                             gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:2240640/70‑1 (MQ=35)
                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:3401/1‑70 (MQ=255)
                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:1080682/1‑70 (MQ=255)
                             gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 >  1:2137436/1‑70 (MQ=255)
                              cACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                               <  1:2844971/71‑1 (MQ=34)
                               aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                                 <  1:924353/68‑1 (MQ=35)
                                   cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAacac                            <  1:1832262/69‑1 (MQ=39)
                                        aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                     >  1:2876018/1‑71 (MQ=255)
                                             cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg                <  1:738200/71‑1 (MQ=255)
                                               gCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGcc              <  1:2290805/71‑1 (MQ=255)
                                                 aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg             >  1:2734822/1‑70 (MQ=255)
                                                 aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg             >  1:1210600/1‑70 (MQ=255)
                                                 aTGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGc            <  1:2549879/71‑1 (MQ=255)
                                                   gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                           <  1:3169475/38‑1 (MQ=40)
                                                   gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTcc                                           <  1:2171599/38‑1 (MQ=40)
                                                       cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCt                                        <  1:2088105/37‑1 (MQ=40)
                                                       cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:2632489/1‑71 (MQ=255)
                                                       cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:2198079/1‑71 (MQ=255)
                                                       cacaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTc      >  1:1953737/1‑71 (MQ=255)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1116067/1‑42 (MQ=38)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:1111019/1‑42 (MQ=38)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:270662/1‑42 (MQ=38)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:2749400/1‑42 (MQ=38)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:2906489/1‑42 (MQ=38)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaa                                  >  1:2525073/1‑42 (MQ=38)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                           >  1:2206199/1‑49 (MQ=255)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     <  1:571480/71‑1 (MQ=255)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     <  1:1727473/71‑1 (MQ=255)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     <  1:704749/71‑1 (MQ=255)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     >  1:1319868/1‑71 (MQ=255)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     <  1:2304925/71‑1 (MQ=255)
                                                        acaGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCg     <  1:2373047/71‑1 (MQ=255)
                                                          aGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGt   >  1:1626919/1‑71 (MQ=255)
                                                           gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                     <  1:1374289/52‑1 (MQ=255)
                                                           gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                     <  1:1168686/52‑1 (MQ=255)
                                                           gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                     <  1:307006/52‑1 (MQ=255)
                                                           gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt                     <  1:3066404/52‑1 (MQ=255)
                                                           gCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTg  <  1:2807555/71‑1 (MQ=255)
                                                            |                                                                     
CCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCGCAGCTTCGGTG  >  minE/2359816‑2359945

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