Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,763,5560TC51.3% ‑1.4 / 174.8 117noncoding (288/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (20/37);  new base C (11/49);  total (31/86)
Fisher's exact test for biased strand distribution p-value = 5.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 122 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGG  >  minE/1763502‑1763623
                                                      |                                                                   
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                      <  1:3989165/70‑1 (MQ=14)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                      <  1:3050573/70‑1 (MQ=14)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                      <  1:2086299/70‑1 (MQ=14)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                      <  1:2027441/70‑1 (MQ=14)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCtt                                                      <  1:1703884/70‑1 (MQ=14)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                     <  1:4260170/71‑1 (MQ=17)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTg                                                     <  1:4750009/71‑1 (MQ=17)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:2700721/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:4665251/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:3983000/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:3774416/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:366732/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:3525553/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:693289/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:1024923/71‑1 (MQ=34)
  aGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGt                                                   <  1:928629/71‑1 (MQ=34)
        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                <  1:2770062/68‑1 (MQ=25)
        gACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                <  1:4612295/68‑1 (MQ=25)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:2661925/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:1252008/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:1192494/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:1169012/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:4786864/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:3111503/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:3724679/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:2038277/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:1116522/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:3446000/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:136882/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:4843616/71‑1 (MQ=21)
            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                         <  1:684751/71‑1 (MQ=34)
            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                         <  1:3319267/71‑1 (MQ=34)
            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                         <  1:342974/71‑1 (MQ=34)
            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                         <  1:922475/71‑1 (MQ=34)
            gtgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTAc                                         <  1:3004397/71‑1 (MQ=34)
             tgtCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAg                                                <  1:2973196/63‑1 (MQ=25)
                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTg                                                  >  1:1893443/1‑57 (MQ=11)
                 tcAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTGGTg                                                  >  1:3580532/1‑57 (MQ=11)
                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:2570723/1‑70 (MQ=35)
                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:1588333/1‑70 (MQ=35)
                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:4139508/1‑70 (MQ=35)
                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:4220991/1‑70 (MQ=35)
                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:2070195/1‑70 (MQ=35)
                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:2675008/1‑70 (MQ=35)
                  cAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccac                                    >  1:697721/1‑70 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATGCTCTCAGACCAGCCAGGGATCGTCGCCTTGGGGAGCCGt                                            <  1:3877680/61‑1 (MQ=11)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGATCCGTTACCccacc                                   <  1:3542744/70‑1 (MQ=25)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                   <  1:3940821/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                   <  1:1699637/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                   <  1:2373558/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                   <  1:486897/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCAGTTACCccacc                                   <  1:3249255/70‑1 (MQ=25)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                   >  1:3933628/1‑70 (MQ=25)
                       ccAGTGTGGCTGGTCCTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                >  1:4818729/1‑69 (MQ=255)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                >  1:4674264/1‑69 (MQ=40)
                       ccAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                                >  1:2064708/1‑69 (MQ=40)
                         aGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                            <  1:446415/71‑1 (MQ=255)
                         aGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGc                            <  1:1636062/71‑1 (MQ=255)
                               gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGGCGTTAcc                                        <  1:958844/53‑1 (MQ=11)
                               gCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTAcc                                        <  1:125961/53‑1 (MQ=255)
                               gCTGGTCATCCTCTCAGACCAGCCAAGGATCGTCGCCTTGGTGAGCCGTTAcc                                        <  1:2321073/53‑1 (MQ=11)
                                  ggTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                            <  1:357544/46‑1 (MQ=255)
                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                >  1:84103/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                >  1:1779672/1‑71 (MQ=255)
                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                <  1:527622/70‑1 (MQ=255)
                                      aTCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg                <  1:446688/70‑1 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca            >  1:1814930/1‑71 (MQ=255)
                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac           <  1:96548/69‑1 (MQ=255)
                                            tcAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcac           <  1:2147622/69‑1 (MQ=255)
                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                   >  1:1543065/1‑44 (MQ=255)
                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc        >  1:632780/1‑71 (MQ=255)
                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc        >  1:4697370/1‑71 (MQ=255)
                                             cAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc        >  1:1752411/1‑71 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc        <  1:2527389/70‑1 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATc        <  1:1375165/70‑1 (MQ=255)
                                                aCCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGa     >  1:2634072/1‑71 (MQ=2)
                                                aCCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGa     >  1:528837/1‑71 (MQ=2)
                                                aCCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGa     >  1:4577701/1‑71 (MQ=2)
                                                aCCAGTTACGGATCGTCGCCTTGGTAGGCCTTTACCCCACCAACTAGCTAATCCGACCTAGGCTCATCTGa     >  1:4895139/1‑71 (MQ=2)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:1096161/71‑1 (MQ=255)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:4131741/71‑1 (MQ=255)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:1668571/71‑1 (MQ=255)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:2286864/71‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGTCACATCCGa     <  1:3394989/71‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:3766304/71‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:3635507/71‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:1472413/71‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:2910186/71‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:404205/71‑1 (MQ=255)
                                                aCCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:1902077/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:3144629/70‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:678764/70‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:1827271/70‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:152033/70‑1 (MQ=255)
                                                 ccagccagGGATAGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGa     <  1:1935913/70‑1 (MQ=255)
                                                   aGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg  <  1:1274318/71‑1 (MQ=255)
                                                   aGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg  <  1:1574070/71‑1 (MQ=255)
                                                   aGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg  <  1:4272066/71‑1 (MQ=255)
                                                   aGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg  <  1:3493018/71‑1 (MQ=255)
                                                   aGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATgg  <  1:3703968/71‑1 (MQ=255)
                                                      |                                                                   
GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGATGG  >  minE/1763502‑1763623

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: