Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,753,9860AT21.1% 77.9 / 15.3 38intergenic (+40/+22)bamE/ratBlipoprotein component of BamABCDE OM biogenesis complex/UPF0125 family protein
Reads supporting (aligned to +/- strand):  ref base A (13/17);  new base T (3/5);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGC  >  NC_000913/2753839‑2754121
                                                                                                                                                   |                                                                                                                                       
cAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAAt                                                                                                                                        <  8:72634/149‑1 (MQ=255)
  gCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGa                                                                                                                                      <  7:87802/149‑1 (MQ=255)
                aaGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttttaac                                                                                                                        <  6:394/149‑4 (MQ=255)
                aaGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCttttttaac                                                                                                                        <  7:117127/149‑4 (MQ=255)
                  ggTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCATTGAGCACCttttttaacgt                                                                                                                      <  1:36289/149‑6 (MQ=255)
                   gtgtAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtctgtct                                                                                                                     <  2:100226/149‑1 (MQ=255)
                              cAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAACCCTTCGCTGAGTGGTACCTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTc                                                                                                          >  3:125498/1‑149 (MQ=255)
                                       acgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGAtt                                                                                                   <  5:110498/147‑1 (MQ=255)
                                        cgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtct                                                                                                                         >  2:124321/1‑124 (MQ=255)
                                        cgctgacCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTtct                                                                                                                         <  1:124321/124‑1 (MQ=255)
                                              ccTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTc                                                                                          <  3:114388/149‑1 (MQ=255)
                                               cTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCg                                                                                         >  6:112770/1‑149 (MQ=255)
                                               cTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCg                                                                                         >  6:110629/1‑149 (MQ=255)
                                                             gtgtgtTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGct                                                                           >  2:72728/1‑149 (MQ=255)
                                                                  ttGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTc                                                                      >  5:67996/1‑149 (MQ=255)
                                                                     aCCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGa                                                                   >  1:37897/1‑149 (MQ=255)
                                                                     aCCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGa                                                                   <  5:112770/149‑1 (MQ=255)
                                                                     aCCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGa                                                                   <  5:110629/149‑1 (MQ=255)
                                                                       cAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATc                                                                 >  2:164694/1‑149 (MQ=255)
                                                                       cAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATc                                                                 >  4:38631/1‑149 (MQ=255)
                                                                       cAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCCTTCGGATc                                                                 >  1:32297/1‑149 (MQ=255)
                                                                         atatTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCgg                                                               <  1:164694/149‑1 (MQ=255)
                                                                                   cAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                     <  3:38631/149‑1 (MQ=255)
                                                                                                    gTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCAccc                                    >  3:93912/1‑149 (MQ=255)
                                                                                                       aCTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                     >  7:55074/1‑129 (MQ=255)
                                                                                                       aCTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAgg                                                     <  8:55074/129‑1 (MQ=255)
                                                                                                              ataaAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCcatc                          <  1:72728/149‑1 (MQ=255)
                                                                                                                aaaGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                           >  2:128186/1‑50 (MQ=255)
                                                                                                                aaaGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCtttttt                                                                                                                           <  1:128186/50‑1 (MQ=255)
                                                                                                                 aaGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATg                       <  6:67996/149‑1 (MQ=255)
                                                                                                                   gTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCCTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCTCCATCATGca                     <  2:32297/149‑1 (MQ=255)
                                                                                                                                 acagaaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTcc                                        >  4:32708/5‑116 (MQ=255)
                                                                                                                                 acagaaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTcc                                        <  3:32708/112‑1 (MQ=255)
                                                                                                                                   agaaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAAtctc                                         <  2:59142/111‑1 (MQ=255)
                                                                                                                                   agaaaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGAGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAAtctc                                         >  1:59142/3‑113 (MQ=255)
                                                                                                                                      aaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCaa                                                            <  8:131481/91‑1 (MQ=255)
                                                                                                                                      aaaaaGGTGCTCATTGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTATGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCaa                                                            >  7:131481/1‑91 (MQ=255)
                                                                                                                                      aaaaaGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGc  <  4:125498/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                       
CAGCAACCAGGTCATGAAGGTGTAACTCAGCAAACGCTGACGCTGACCTTTAACAGTAGCGGTGTGTTGACCAATATTGATAACAAACCTGCGCTGAGTGGTAACTAATAATAAAGTTGCTCTCAAAGACGTTAAAAAAGGTGCTCAATGAGCACCTTTTTTCTGTCTGTTATTTATTCGCTGATTTTTCTGCTCGTTGCCTGCGAAGCTCTTTCGGATCGGCAATGAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGC  >  NC_000913/2753839‑2754121

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑