Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,710,2790AC30.8% 61.2 / 27.2 39intergenic (‑115/‑141)rseD/nadBrpoE leader peptide/L‑aspartate oxidase
Reads supporting (aligned to +/- strand):  ref base A (27/0);  new base C (0/12);  total (27/12)
Fisher's exact test for biased strand distribution p-value = 2.56e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTA  >  NC_000913/2710214‑2710346
                                                                 |                                                                   
gAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccacc                                                                   >  1:2412413/1‑68 (MQ=255)
  aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTccaccct                                                                 >  1:1281049/1‑66 (MQ=255)
  aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                 >  1:1595095/1‑68 (MQ=255)
  aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                 >  1:1849174/1‑68 (MQ=255)
  aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                 >  1:3934739/1‑68 (MQ=255)
  aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                 >  1:2120301/1‑68 (MQ=255)
  aaTGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCAccag                                                                 >  1:2058262/1‑68 (MQ=255)
    tGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCTTTTTTATTTCACCACCTAATCCTCCAccagcc                                                               >  1:25715/1‑68 (MQ=255)
           ttAACTGCCTGAAGAGTACCCCAGCGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACt                                                        <  1:1033878/68‑1 (MQ=255)
              aCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACNTAATCCTCCACTAGCCAGTAACTtct                                                     >  1:3257127/1‑68 (MQ=255)
              aCTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtct                                                     >  1:2886310/1‑68 (MQ=255)
              aCTGCCTGAAGAGAAACCCACCGGCTTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTtct                                                     <  1:3520121/68‑1 (MQ=255)
               cTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtctc                                                    >  1:1228036/1‑68 (MQ=255)
               cTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTtctc                                                    >  1:249811/1‑68 (MQ=255)
                        gagTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCAAGTAACTTCTCTTTTTCTcg                                           >  1:3324935/1‑68 (MQ=255)
                        gagCAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcg                                           >  1:3607732/1‑68 (MQ=255)
                         agTAACCCAGCGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                          >  1:2576823/1‑68 (MQ=255)
                         agTAACCCACCGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTtctccc                                          >  1:474534/1‑66 (MQ=255)
                         agTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                          >  1:3061393/1‑68 (MQ=255)
                         agTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                          >  1:2671031/1‑68 (MQ=255)
                         agTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                          >  1:2529833/1‑68 (MQ=255)
                         agTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                          >  1:4061691/1‑68 (MQ=255)
                         agTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgc                                          >  1:3952287/1‑68 (MQ=255)
                         agTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCACCCAGTAACTTCTCTTTTTCTcgc                                          >  1:3032542/1‑68 (MQ=255)
                          gTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgcc                                         >  1:2232818/1‑68 (MQ=255)
                            accccATCGGCCTTTTTTTTCCCCCCCCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTcgccgc                                       <  1:389164/66‑1 (MQ=255)
                            accccAACGGCTTTTTTTTTCCCCCCCCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTcgccgc                                       <  1:3433976/66‑1 (MQ=255)
                            aaCCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTcgccgc                                       >  1:2027520/1‑68 (MQ=255)
                             ccccACCGGCTTTTTTATTCCCCCACCAAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCgccgcc                                      <  1:2383723/67‑1 (MQ=255)
                             ccccACCGCCTTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCgccgcc                                      <  1:3500074/67‑1 (MQ=255)
                             ccccACCGCCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCgccgcc                                      <  1:1140437/67‑1 (MQ=255)
                             ccccAACGGCTTTTTTTTTCCCCCCCATAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCgccgcc                                      <  1:3967633/67‑1 (MQ=255)
                               ccACCGGCCTTTTTATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                    <  1:2075809/68‑1 (MQ=255)
                               ccAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                    >  1:561817/1‑68 (MQ=255)
                               ccAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCt                                    >  1:3644615/1‑68 (MQ=255)
                                 aaCGACTTTTTTATTACCCCCCCTACCCCTCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGc                                  <  1:2186798/68‑1 (MQ=255)
                                 aaCCGCCTTTTTACTTCACCACCTACCCCCCCCCCAGCCAGTAACTTCTNTTTTTCTCGCCGCCCTGc                                  <  1:3099577/68‑1 (MQ=255)
                                   cGGCTTTTTTTTTCCCCCACCTAACCCCCCCCCAGCCAGTAAGTTCTCTTTTTCTCGCCGCCCTGCGt                                <  1:2743552/68‑1 (MQ=255)
                                         ttttATTCCCCCCCCTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt                          <  1:1959891/68‑1 (MQ=255)
                                         ttttATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt                          <  1:493954/68‑1 (MQ=255)
                                         ttttATTCCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt                          <  1:2243208/68‑1 (MQ=255)
                                         tttaATTTCCCCCCCAAATCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgt                          <  1:3688861/68‑1 (MQ=255)
                                          tttATTTCCCCCCCTAACCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCgtg                         <  1:295385/68‑1 (MQ=255)
                                                    caccTACCCCCCCCCCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACt               <  1:1157508/68‑1 (MQ=255)
                                                            cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCCGCGTGTTTAGCAACTGTAACaaa       >  1:3219024/1‑68 (MQ=255)
                                                            cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACaaa       >  1:74083/1‑68 (MQ=255)
                                                            cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACaaa       >  1:3337907/1‑68 (MQ=255)
                                                            cctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACaaa       >  1:3120602/1‑68 (MQ=255)
                                                             ctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAat      >  1:4194666/1‑68 (MQ=255)
                                                             ctccACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAat      >  1:538426/1‑68 (MQ=255)
                                                                 accaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTa  >  1:2550286/1‑68 (MQ=255)
                                                                 accaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTa  >  1:3368519/1‑68 (MQ=255)
                                                                 accaGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTa  >  1:3345103/1‑68 (MQ=255)
                                                                 |                                                                   
GAAATGCCCATTTAACTGCCTGAAGAGTAACCCAACGGCCTTTTTATTTCACCACCTAATCCTCCACCAGCCAGTAACTTCTCTTTTTCTCGCCGCCCTGCGTCAGCGTGTTTAGCAACTGTAACAAATATTA  >  NC_000913/2710214‑2710346

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑