Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,631,7720AC35.5% 47.0 / 13.8 31intergenic (‑267/‑141)rpoE/nadBRNA polymerase, sigma 24 (sigma E) factor/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base A (17/3);  new base C (0/11);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 4.30e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG  >  REL606/2631697‑2631845
                                                                            |                                                                         
tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca                                                                           >  2:320741/1‑76 (MQ=255)
                cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct                                                           >  2:1220652/1‑76 (MQ=255)
                  cATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTctct                                                         >  2:53719/1‑76 (MQ=255)
                     ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtttt                                                      >  1:795012/1‑76 (MQ=255)
                     ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTATTTCACCACCTGATCCTCCACCAGCTAGTAACTTCTCttttt                                                     <  1:199075/76‑1 (MQ=255)
                      tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCCCCTAATCCCCCCCCAGCCAGTAACTTCTCttttt                                                     <  1:199059/76‑1 (MQ=255)
                      tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCttttt                                                     <  2:929104/76‑1 (MQ=255)
                      tAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCttttt                                                     >  1:401265/1‑76 (MQ=255)
                      tAACGGCCTGAAGAGTACCCCACCGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCttttt                                                     <  1:199056/76‑1 (MQ=255)
                       aaCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtc                                                    >  2:611233/1‑76 (MQ=255)
                               gAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCgccgcc                                            >  2:622005/1‑76 (MQ=255)
                                 agagTACCCCAACGCCCTTTTTATTTCCCCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                          <  2:1137350/76‑1 (MQ=255)
                                 agagTAACCCACCGGCTTTTTTATTTCCCCCCCTATTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                          <  1:143879/76‑1 (MQ=255)
                                 agagTAACCCACCGGCTTTTTTATTTCCCCACCTAACCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                          <  2:1137344/76‑1 (MQ=255)
                                  gagTACCCCACCGGCCTTTTTATTCCCCCCCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg                                         <  2:250505/76‑1 (MQ=255)
                                     tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt                                      >  1:424181/1‑76 (MQ=255)
                                       ccccAACGCCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCa                                    <  2:951063/75‑1 (MQ=255)
                                              ggCCTTTTTATTTCACCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt                             <  2:597861/76‑1 (MQ=255)
                                               gCCTTTTTATTCCCCCCCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt                            <  1:1121564/76‑1 (MQ=255)
                                               gCCTTTTTATTCCCCCCCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGttt                            <  2:241203/76‑1 (MQ=255)
                                                cttttttATTTCCCCCCCTAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa                           <  2:724622/74‑1 (MQ=255)
                                                ccTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa                           >  2:147167/1‑76 (MQ=255)
                                                 tttttttttTCCCCAACCAAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTa                           <  2:724634/75‑1 (MQ=255)
                                                 ttttttATCTCCCCCCCAACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg                          <  2:812286/75‑1 (MQ=255)
                                                 cTTTTTATTCCCCCCCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg                          <  2:816126/76‑1 (MQ=255)
                                                 cTTTTTATTCCCCCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAg                          <  2:812297/76‑1 (MQ=255)
                                                     ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc                      >  2:519991/1‑76 (MQ=255)
                                                     ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc                      >  2:668921/1‑76 (MQ=255)
                                                       aTTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTg                    >  1:1109096/1‑76 (MQ=255)
                                                          tCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTaa                 >  2:902577/1‑76 (MQ=255)
                                                            accaccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCCCTGTAAcg               >  1:1012117/1‑75 (MQ=255)
                                                            accaccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACa               >  1:1012118/1‑76 (MQ=255)
                                                              caccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACaaa             >  1:696142/1‑76 (MQ=255)
                                                                cccTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAat            <  2:1146207/75‑1 (MQ=255)
                                                                 ccTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAata           >  2:357576/1‑76 (MQ=255)
                                                                   tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt         >  1:858796/1‑76 (MQ=255)
                                                                   tAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATAtt         >  1:621140/1‑76 (MQ=255)
                                                                     aTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaa       >  1:453527/1‑76 (MQ=255)
                                                                       cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaaa     >  1:1178552/1‑76 (MQ=255)
                                                                          ccaccaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAg  >  2:959254/1‑76 (MQ=255)
                                                                            |                                                                         
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAAATAG  >  REL606/2631697‑2631845

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: