Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi1,018,9860AT27.5% 83.9 / 42.6 51intergenic (+49/+27)matP/ompATer macrodomain organizer matS‑binding protein/outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base A (21/16);  new base T (7/7);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGT  >  NC_000913_3_bme_pgi/1018843‑1019126
                                                                                                                                                |                                                                                                                                            
ttCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGA‑AAAAAAAACCCCGCAGCAGCgggg                                                                                                                                        <  4:227888/150‑1 (MQ=255)
 tCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGt                                                                                                                                       <  3:292985/150‑1 (MQ=255)
 tCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGt                                                                                                                                       <  1:190397/150‑1 (MQ=255)
  cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtt                                                                                                                                      <  1:332288/150‑1 (MQ=255)
  cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGt                                                                                                                                       >  2:33957/1‑149 (MQ=255)
  cAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGt                                                                                                                                       <  1:33957/149‑1 (MQ=255)
      tGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCTGCTGCGGGGtttttt                                                                                                                                  <  4:99073/150‑2 (MQ=255)
          taaaGATGCGGAAAAAAAAGAGAAATACGGGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGttttttttgc                                                                                                                              <  2:505849/148‑6 (MQ=255)
               aTGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGAc                                                                                                                          >  3:213401/1‑150 (MQ=255)
                    gAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGAGCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAAc                                                                                                                    <  2:215354/150‑1 (MQ=255)
                           aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGcc                                                                                                             >  1:201900/1‑150 (MQ=255)
                           aaGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGA‑AAAAAAAACCCCGCAGCACCGGGGGTTTTCTACCAGATGAGAACGTAAGcc                                                                                                             >  4:338172/1‑150 (MQ=255)
                                    aCGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                    <  3:103971/150‑1 (MQ=255)
                                    aCGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGa                                                                                                    >  4:103971/1‑150 (MQ=255)
                                        gAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTa                                                                                                >  4:10087/1‑150 (MQ=255)
                                                  tcttctCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCttt                                                                                      >  3:79931/1‑150 (MQ=255)
                                                                ggATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                                                        >  1:498290/1‑150 (MQ=255)
                                                                 gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTa                                                                                                                                <  2:458154/93‑1 (MQ=255)
                                                                 gATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGATTTTCTa                                                                                                                                >  1:458154/1‑93 (MQ=255)
                                                                                   ggTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                              >  3:191184/1‑77 (MQ=255)
                                                                                   ggTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTAcc                                                                                                                              <  4:191184/77‑1 (MQ=255)
                                                                                    gTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                >  2:19414/1‑90 (MQ=255)
                                                                                    gTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTaa                                                                                                                <  1:19414/90‑1 (MQ=255)
                                                                                      aaGGAATAACTGACGAAAGTCAGTTCCATTTACTAAAGCC‑AAAACAAACCCCGCGGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCgg                                                  <  3:338172/150‑1 (MQ=255)
                                                                                      aaGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                                                        <  2:122144/128‑1 (MQ=255)
                                                                                      aaGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACtt                                                                        >  1:122144/1‑128 (MQ=255)
                                                                                        ggAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGa                                                 <  4:213401/150‑1 (MQ=255)
                                                                                         gAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                               >  3:120905/1‑150 (MQ=255)
                                                                                                tGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAg                                        >  2:406074/1‑150 (MQ=255)
                                                                                                     aaaGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCt                                                                               >  2:6313/1‑106 (MQ=255)
                                                                                                     aaaGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCt                                                                               <  1:6313/106‑1 (MQ=255)
                                                                                                             ttCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcag                           >  3:156707/1‑150 (MQ=255)
                                                                                                                 aTTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcag                           >  1:156165/1‑146 (MQ=255)
                                                                                                                 aTTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTgcag                           <  2:156165/146‑1 (MQ=255)
                                                                                                                  tttACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGc                      >  2:577666/1‑150 (MQ=255)
                                                                                                                   ttACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAgg                                           >  2:89013/1‑128 (MQ=255)
                                                                                                                   ttACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAgg                                           <  1:89013/128‑1 (MQ=255)
                                                                                                                        aaaGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTc                <  3:146903/150‑1 (MQ=255)
                                                                                                                                aaaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGt         >  1:395819/1‑150 (MQ=255)
                                                                                                                                 aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTc        <  4:156707/150‑1 (MQ=255)
                                                                                                                                 aaaaaaCCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTc        >  1:253321/1‑150 (MQ=255)
                                                                                                                                  aaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                               <  3:302655/110‑1 (MQ=255)
                                                                                                                                  aaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                               >  4:302653/1‑110 (MQ=255)
                                                                                                                                  aaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                               >  4:302655/1‑110 (MQ=255)
                                                                                                                                  aaaaaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc                                               <  3:302653/110‑1 (MQ=255)
                                                                                                                                     aaCCCCGCTGCTGCGGGGTTTTTTTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                 <  3:313824/121‑1 (MQ=255)
                                                                                                                                     aaCCCCGCTGCTGCGGGGTTTTTTTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCa                                 >  4:313824/1‑121 (MQ=255)
                                                                                                                                     aaCCCCGCTGCTGCGGGGTTTTTTTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAg    >  1:153344/1‑150 (MQ=255)
                                                                                                                                     aaCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAg    >  1:556891/1‑150 (MQ=255)
                                                                                                                                      aCCCCGCTGCTGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAgg   >  3:77539/1‑150 (MQ=255)
                                                                                                                                       ccccGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGgt  >  2:374211/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGC‑AAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGT  >  NC_000913_3_bme_pgi/1018843‑1019126

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