Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,763,5560TC36.8% 145.3 / 171.2 153noncoding (288/1542 nt)rrsG16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (52/44);  new base C (24/32);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 174 total aligned reads displayed.

GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAT  >  minE/1763502‑1763621
                                                      |                                                                 
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                   >  1:365339/1‑71 (MQ=11)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                   >  1:2489549/1‑71 (MQ=11)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCGTCCTCTCAGACCAGCTACTGATCGTCGCCTTg                                                   >  1:186385/1‑71 (MQ=11)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                   <  1:1144300/71‑1 (MQ=32)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                   <  1:3537856/71‑1 (MQ=32)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                   <  1:307566/71‑1 (MQ=32)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                   <  1:1383399/71‑1 (MQ=32)
ggAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                   <  1:3019994/71‑1 (MQ=32)
 gAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                   <  1:158302/70‑1 (MQ=33)
         aCCGTGTCTTAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                          <  1:596045/71‑1 (MQ=255)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                          <  1:2726531/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                          <  1:913546/71‑1 (MQ=21)
         aCCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGt                                          <  1:3404556/71‑1 (MQ=21)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 >  1:777098/1‑70 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 <  1:1813721/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 <  1:155075/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 <  1:899957/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 <  1:1061374/70‑1 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 >  1:522071/1‑70 (MQ=35)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 >  1:4141088/1‑70 (MQ=25)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 >  1:2310881/1‑70 (MQ=25)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 >  1:2795622/1‑70 (MQ=25)
                   aGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 >  1:1577904/1‑70 (MQ=25)
                   aGTTCAAGTGTGGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCccacc                                 >  1:2001106/1‑70 (MQ=255)
                      tCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACcaac                              >  1:1214281/1‑70 (MQ=40)
                            gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTaa                       >  1:164387/1‑71 (MQ=255)
                            gtgGCTGGTCATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTa                        >  1:3348190/1‑70 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:2660384/1‑68 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:2527480/1‑68 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:3365900/1‑68 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:3149314/1‑68 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:3221367/1‑68 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:2049019/1‑68 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:956198/1‑68 (MQ=255)
                             tgGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:1104555/1‑68 (MQ=255)
                             tgGCTGGTCATCATCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCt                         >  1:1615286/1‑68 (MQ=255)
                                    tcatcaTCTCAGACCAGCTAGGGATCGTCGCCTTg                                                   <  1:646627/35‑1 (MQ=255)
                                    tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt               >  1:477485/1‑71 (MQ=255)
                                    tCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCt               >  1:1974745/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:4336705/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:3985558/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:3483714/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:3675568/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:1913105/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:2457682/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCGCATCTg              >  1:942898/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:2917760/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:1014755/1‑71 (MQ=255)
                                     cATCCTCTCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTg              >  1:326349/1‑71 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca          <  1:2623834/71‑1 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca          <  1:3487716/71‑1 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca          <  1:409808/71‑1 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca          <  1:2258352/71‑1 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca          <  1:1505091/71‑1 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca          <  1:832450/71‑1 (MQ=255)
                                         ctctCAGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGca          <  1:485720/71‑1 (MQ=255)
                                              aGACCTGCTCGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:1540861/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCTCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:3818140/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCATCAAGCTAATCCCATCTGGGCACATcc     <  1:4053839/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                >  1:167492/1‑60 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:4221963/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:2999635/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:4374884/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:2440632/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:36014/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:3607227/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:471709/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:1896773/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:3850626/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGcaca        >  1:1500066/1‑68 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:4159663/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:4271873/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:45495/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:57270/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:2952343/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:402312/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:1228494/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:3667389/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:1729394/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:1773448/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:2361642/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:2872968/71‑1 (MQ=255)
                                              aGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:2402242/71‑1 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATc                >  1:2121248/1‑60 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:1046977/71‑1 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     >  1:1756336/1‑71 (MQ=255)
                                              aGACCAGCCAGGGATCGTCGCCTTGGTGAGCCGGTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:4148146/71‑1 (MQ=255)
                                                aCCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATcc     <  1:4306017/69‑1 (MQ=255)
                                                 ccagccagGGATCTTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:2178960/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:2612567/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:2283115/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:2886065/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:2102053/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:3045008/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:3466545/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:1283650/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:3881216/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:1207817/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:1072228/71‑1 (MQ=255)
                                                 ccagccagGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAt  <  1:3991338/71‑1 (MQ=255)
                                                      |                                                                 
GGAGTCTGGACCGTGTCTCAGTTCCAGTGTGGCTGGTCATCCTCTCAGACCAGCTAGGGATCGTCGCCTTGGTGAGCCGTTACCCCACCAACAAGCTAATCCCATCTGGGCACATCCGAT  >  minE/1763502‑1763621

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: