Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,359,8670GA35.9% 99.5 / 79.1 78noncoding (1228/2905 nt)rrlA23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (18/10);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/2359815‑2359934
                                                    |                                                                   
cccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   >  1:1544043/1‑71 (MQ=9)
 ccGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTca                                                   <  1:240262/70‑1 (MQ=31)
   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTGCAGCATGCCTCACAGCACACCTTCACAg                                                >  1:190410/1‑71 (MQ=21)
   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCATCATGCCTCACAGCACACCTTCACAg                                                >  1:1130353/1‑71 (MQ=21)
   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAg                                                >  1:327724/1‑71 (MQ=34)
   gCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAg                                                <  1:168356/71‑1 (MQ=255)
    cTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                               >  1:1211336/1‑71 (MQ=34)
     tttATCGTTACTTATGACAGCATTCTCACTTATGATACCTCCAGCATGCCTCACAGCACACCTTCACAgg                                               >  1:1192818/1‑70 (MQ=11)
      ttATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                             >  1:881155/1‑71 (MQ=255)
        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTcaca                                                 >  1:2320836/1‑65 (MQ=35)
        aTCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCtt                                            <  1:509702/70‑1 (MQ=255)
         tCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTAc                                          <  1:765026/71‑1 (MQ=255)
           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                         <  1:1691612/70‑1 (MQ=35)
           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:2076606/71‑1 (MQ=34)
           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAg                                        <  1:184195/71‑1 (MQ=34)
           gTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTa                                           <  1:2046372/68‑1 (MQ=255)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                              >  1:1846034/1‑62 (MQ=35)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                              >  1:1751694/1‑62 (MQ=35)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                              >  1:1580867/1‑62 (MQ=35)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGc                                              >  1:641685/1‑62 (MQ=35)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGaa                                      <  1:1016382/70‑1 (MQ=255)
              aCTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAAc                                     >  1:2258860/1‑71 (MQ=255)
               cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACg                                    <  1:271315/71‑1 (MQ=34)
               cTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   >  1:2048281/1‑56 (MQ=11)
                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                 <  1:534699/70‑1 (MQ=35)
                   tGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTc                                 <  1:1128502/70‑1 (MQ=35)
                    gTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                               >  1:394107/1‑71 (MQ=34)
                    gTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                               >  1:2325664/1‑71 (MQ=34)
                    gTCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                               >  1:1670198/1‑71 (MQ=255)
                     tCAGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTca                                                   <  1:478644/50‑1 (MQ=11)
                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTa                            >  1:259972/1‑71 (MQ=34)
                       aGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACCCTCCCCt                             <  1:1267579/70‑1 (MQ=25)
                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                  >  1:168796/1‑65 (MQ=11)
                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                  >  1:165807/1‑65 (MQ=11)
                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCt                                  >  1:612125/1‑65 (MQ=11)
                       aGCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTa                            >  1:1028352/1‑71 (MQ=255)
                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACa                                         >  1:1259725/1‑57 (MQ=35)
                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                           <  1:491382/71‑1 (MQ=34)
                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                           <  1:471056/71‑1 (MQ=34)
                        gCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAc                           <  1:851242/71‑1 (MQ=34)
                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCt                                             <  1:1916265/53‑1 (MQ=11)
                        gCATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAc                           >  1:434990/1‑71 (MQ=255)
                         cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:45435/1‑71 (MQ=255)
                         cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:2083635/1‑71 (MQ=255)
                         cATTCGCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:1330736/1‑71 (MQ=255)
                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                          <  1:1897716/52‑1 (MQ=35)
                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                          <  1:353913/52‑1 (MQ=35)
                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                          <  1:458668/52‑1 (MQ=35)
                            tCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTAc                                          <  1:1464835/52‑1 (MQ=35)
                            tCGCACATCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                               >  1:1800/1‑63 (MQ=11)
                             cGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCt                                  >  1:1182083/1‑59 (MQ=35)
                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      <  1:1481976/70‑1 (MQ=35)
                              gCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      <  1:631445/70‑1 (MQ=35)
                              gCACTTCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaac                      <  1:1566882/70‑1 (MQ=255)
                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                               >  1:1127678/1‑59 (MQ=35)
                                aCTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTccc                               >  1:1006831/1‑59 (MQ=35)
                                aCTTCTCATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTc                                 >  1:3060/1‑57 (MQ=11)
                                  ttCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTCCCCTACCCaacaacg                  >  1:1560134/1‑69 (MQ=11)
                                   tCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                >  1:705443/1‑71 (MQ=39)
                                   tCTGATACCTCCAGCATACCTCACAGTACACCTTCACAGGCTTACa                                         >  1:1941963/1‑46 (MQ=11)
                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                >  1:867568/1‑70 (MQ=40)
                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                <  1:1014025/70‑1 (MQ=40)
                                    cTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAAcaca                <  1:1871677/70‑1 (MQ=40)
                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATa              <  1:2430768/69‑1 (MQ=40)
                                       aTACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgt            >  1:742668/1‑71 (MQ=255)
                                        tACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtg           <  1:1413208/71‑1 (MQ=255)
                                         aCCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAgtgt          >  1:979347/1‑71 (MQ=255)
                                          ccTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                       <  1:98624/41‑1 (MQ=14)
                                          ccTCCAGCATACCTCACAGTACACCTTCACAGGCTTACAGa                                       <  1:416633/41‑1 (MQ=14)
                                           cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCcaacaa                    >  1:1684020/1‑59 (MQ=35)
                                           cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg        <  1:1782875/71‑1 (MQ=255)
                                           cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg        <  1:1189538/71‑1 (MQ=255)
                                           cTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCg        <  1:2114114/71‑1 (MQ=255)
                                             ccAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt      >  1:1578569/1‑71 (MQ=255)
                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCt      <  1:1479757/70‑1 (MQ=255)
                                              cAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTg     >  1:2131936/1‑71 (MQ=255)
                                               aGCATACCTCACAGTACACCTTCACAGGCTTACAGAACGCTccc                               >  1:2074426/1‑44 (MQ=11)
                                                 cATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCg  <  1:1723320/71‑1 (MQ=255)
                                                    gCCTCACAGCACCCCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:1127842/1‑44 (MQ=255)
                                                    gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:2112582/1‑44 (MQ=38)
                                                    gCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTAcc                          >  1:946393/1‑44 (MQ=38)
                                                    |                                                                   
CCCGCTTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCACAGGCTTACAGAACGCTCCCCTACCCAACAACACATAGTGTCGCTGCCG  >  minE/2359815‑2359934

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: