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

TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAA  >  REL606/2631697‑2631841
                                                                           |                                                                     
tATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTcca                                                                       <  1:1417544/76‑1 (MQ=255)
 aTCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccac                                                                      <  1:1040751/76‑1 (MQ=255)
    tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                   >  2:836986/1‑76 (MQ=255)
    tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                   <  2:1944362/76‑1 (MQ=255)
    tgtgTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                   <  1:1537594/76‑1 (MQ=255)
      tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc                                                                 >  2:1045965/1‑76 (MQ=255)
      tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc                                                                 <  1:257267/76‑1 (MQ=255)
      tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc                                                                 <  1:1466561/76‑1 (MQ=255)
      tgtAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagcc                                                                 >  1:1128393/1‑76 (MQ=255)
        tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag                                                               <  2:46719/76‑1 (MQ=255)
        tAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccagccag                                                               >  2:675820/1‑76 (MQ=255)
           aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa                                                            >  2:322389/1‑76 (MQ=255)
           aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa                                                            >  1:2237022/1‑76 (MQ=255)
           aaaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTaa                                                            >  1:2054098/1‑76 (MQ=255)
           aaaTACCTATTTACCTGCTGGAAGACAACCCCAACGCCCTTTTTATTTCTCCCCCTAATCCCCCCCCAGCCAGTaa                                                            <  2:1625009/76‑1 (MQ=255)
              tGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCCCCTAATCCTCCACCAGCCAGTAACtt                                                         <  2:532275/76‑1 (MQ=255)
               gCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtc                                                        >  1:1695750/1‑76 (MQ=255)
                cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct                                                       >  2:441655/1‑76 (MQ=255)
                cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct                                                       >  1:784661/1‑76 (MQ=255)
                cccATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct                                                       >  2:1986730/1‑76 (MQ=255)
                  cATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTctct                                                     >  1:1921813/1‑76 (MQ=255)
                   aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtt                                                    >  2:1230386/1‑76 (MQ=255)
                   aTTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCtt                                                    >  1:2096419/1‑76 (MQ=255)
                    tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTTTTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCttt                                                   <  2:2132573/76‑1 (MQ=255)
                    tttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCCCCACCAGCCAGTAACTTCTCttt                                                   <  2:1878111/76‑1 (MQ=255)
                     ttAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTCCACCACCAAATCCCCCCCCAGCCAGTAACTTCTCtttt                                                  <  2:2048529/76‑1 (MQ=38)
                     taaaCTGCCAGAAGAGTACCCCACCGGCCTTTTTATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCtttt                                                  <  2:776654/74‑1 (MQ=255)
                       aaCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtc                                                >  2:697746/1‑76 (MQ=255)
                        aCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtct                                               >  1:110044/1‑76 (MQ=255)
                        aCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtct                                               >  1:394430/1‑76 (MQ=255)
                          cgCCTGAAGAGTACCCCAACGGCCTTTTTATTTCCCCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTcg                                             <  1:811305/75‑1 (MQ=255)
                           gCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                            >  1:1898896/1‑76 (MQ=255)
                           gCCTGAAGAGTAACCCAACGGCCTTTTTATGTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                            >  2:1192994/1‑76 (MQ=255)
                            ccTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgcc                                           >  1:1442049/1‑76 (MQ=255)
                              tGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccgc                                         >  2:2178898/1‑76 (MQ=255)
                              tGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccgc                                         >  2:241875/1‑76 (MQ=255)
                                 agagTACCCCAACGGCCTTTTTATTTCCCCCCCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                      <  1:2138503/76‑1 (MQ=255)
                                 agagTACCCCAACGGCCTTTTTATTCCACCCCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                      <  1:1597660/76‑1 (MQ=255)
                                 agagTAACCCAACGGCTTTTTTATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                      <  1:714937/76‑1 (MQ=255)
                                 agagTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                      >  1:1411036/1‑76 (MQ=255)
                                  gagTACCCCACCGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg                                     <  1:383527/76‑1 (MQ=255)
                                  gagTACCCCAACGGCCTTTTTATTTCACCCCCTAGCCCCCCTCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTg                                     <  2:943991/76‑1 (MQ=255)
                                   agTACCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGc                                    <  2:1827859/76‑1 (MQ=255)
                                    ccaaCCCAGCGGCCTTTTTTTTCCACCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCg                                   <  1:292334/74‑1 (MQ=255)
                                     tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCgc                                  >  1:1412391/1‑75 (MQ=255)
                                     tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt                                  >  1:1260669/1‑76 (MQ=255)
                                     tAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGt                                  >  1:698294/1‑76 (MQ=255)
                                      aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTc                                 >  2:936904/1‑76 (MQ=255)
                                       ccccACCGGCCTTTTTATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCa                                <  1:1355163/75‑1 (MQ=255)
                                        cccACCGGCCTTTTTATTTCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg                               <  2:1433123/76‑1 (MQ=255)
                                        cccAACGGCCTTTTTATTTCCCCACCTAATCCCCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg                               <  2:721075/76‑1 (MQ=255)
                                        cccAACGCCTTTTTTATTCCCCCACCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAg                               <  2:1476060/76‑1 (MQ=255)
                                          cAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCg                             >  2:498487/1‑76 (MQ=255)
                                           aaCGGCCTTTTTATTCCACCCCCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt                            <  1:1102085/76‑1 (MQ=255)
                                              ggCTTTTTTATTTCCCCACCTAACCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt                         <  2:1747652/76‑1 (MQ=255)
                                              ggCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt                         >  1:302680/1‑76 (MQ=255)
                                              ggCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGtt                         >  2:1248198/1‑76 (MQ=255)
                                               gCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTCCTCGCCGCCCTGCGTCAGCGTGttt                        <  2:104735/76‑1 (MQ=255)
                                                   ttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc                    >  1:423347/1‑76 (MQ=255)
                                                   ttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc                    >  1:2063002/1‑76 (MQ=255)
                                                   ttttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc                    >  1:122849/1‑76 (MQ=255)
                                                   ttttATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGcc                    <  1:1096810/76‑1 (MQ=255)
                                                    tttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCa                   >  1:1026387/1‑76 (MQ=255)
                                                     ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCAc                  >  1:1928840/1‑76 (MQ=255)
                                                     ttATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTCCGTCAGCGCGTTTAGCCAc                  >  2:793728/1‑76 (MQ=255)
                                                        tttCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGt               <  1:361063/76‑1 (MQ=255)
                                                         ttCACCACCTAATCCTCCACCAGCCAGTAACCTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTa              >  1:1588445/1‑76 (MQ=255)
                                                                 cTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat      >  2:535426/1‑76 (MQ=255)
                                                                 cTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAatat      <  1:1481256/76‑1 (MQ=255)
                                                                   aaTCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTa    >  2:222910/1‑76 (MQ=255)
                                                                     tcctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTaaa  >  2:742951/1‑76 (MQ=255)
                                                                           |                                                                     
TATCTGTGTAGAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCCACTGTAACAAATATTAAA  >  REL606/2631697‑2631841

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