Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,591,4980TA21.4% 87.3 / 20.2 42intergenic (‑219/‑159)yjiY/tsrputative transporter/methyl‑accepting chemotaxis protein I, serine sensor receptor
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base A (4/5);  total (20/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 = 9.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGT  >  NC_000913/4591364‑4591638
                                                                                                                                      |                                                                                                                                            
tCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATAttt                                                                                                                                <  1:94751/149‑1 (MQ=255)
       gTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATtgt                                                                                                                         <  6:81943/149‑1 (MQ=255)
       gTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGACGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATtgt                                                                                                                         <  4:102028/149‑1 (MQ=255)
          aTGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGa                                                                                                                       <  2:47186/148‑1 (MQ=255)
                 cAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGTGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTa                                                                                                               >  3:74739/1‑149 (MQ=255)
                 cAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTta                                                                                                                               >  3:65142/1‑133 (MQ=255)
                 cAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTta                                                                                                                               >  3:62811/1‑133 (MQ=255)
                 cAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTta                                                                                                                               <  4:62811/133‑1 (MQ=255)
                 cAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTta                                                                                                                               <  4:65142/133‑1 (MQ=255)
                          gcgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGTGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAaataa                                                                                                      <  4:74739/149‑1 (MQ=255)
                          gcgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCGAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAaataa                                                                                                      >  1:119414/1‑149 (MQ=255)
                           cgcgcgTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAaataat                                                                                                     <  2:119414/149‑1 (MQ=255)
                                          gTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAtt                                                                                        >  2:70433/1‑147 (MQ=255)
                                            aaaTAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                      <  7:16105/147‑1 (MQ=255)
                                            aaaTAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCa                                                                                      >  8:16105/1‑147 (MQ=255)
                                                aaCTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAaataat                                                                                                     >  4:44296/1‑128 (MQ=255)
                                                aaCTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAaataat                                                                                                     <  3:44296/128‑1 (MQ=255)
                                                   tGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtat                                                                                                                              >  1:102954/1‑100 (MQ=38)
                                                   tGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATATGAATCACATATTtat                                                                                                                              <  2:102954/100‑1 (MQ=38)
                                                         gTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAAt                                                                       >  6:116060/1‑149 (MQ=255)
                                                                        tAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAAtc                                                        >  2:13111/1‑149 (MQ=255)
                                                                              gTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                    <  2:188049/115‑1 (MQ=255)
                                                                              gTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCaga                                                                                    >  1:188049/1‑115 (MQ=255)
                                                                              gTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCcgcg                                                  <  1:13111/149‑1 (MQ=255)
                                                                               tCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCgg                                                 <  5:116060/149‑1 (MQ=255)
                                                                                         ggTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCAt                                        >  3:84176/1‑148 (MQ=255)
                                                                                                   cAAAGGCGTCGTTTAAGTGACGATAAATATGTGGTTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGTGATATTCATAAAGTTtttca                             >  7:36343/1‑148 (MQ=255)
                                                                                                   cAAAGGCGTCGTTTAAGTGACGATAAATATGTGGTTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGttttt                               <  8:36343/147‑1 (MQ=255)
                                                                                                   cAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAtt                                                                                        >  5:58527/1‑90 (MQ=255)
                                                                                                   cAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAAtt                                                                                        <  6:58527/90‑1 (MQ=255)
                                                                                                      aGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATaaaa                                                                        <  4:141939/103‑1 (MQ=255)
                                                                                                      aGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATaaaa                                                                        >  3:141939/1‑103 (MQ=255)
                                                                                                                   gTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAAt             >  5:50724/1‑149 (MQ=255)
                                                                                                                      aataTAAATATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGt                                                            <  5:8894/96‑1 (MQ=37)
                                                                                                                      aataTAAATATGTGATATGAATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGt                                                            >  6:8894/4‑99 (MQ=37)
                                                                                                                      aCGATAAATATGTGATTCATATCACATATTTATATTGTGATTAATTTATGTATATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGAGGGATATTCATAAAGTTTTTCCTTTCCATGCCGAAAATCt           >  4:43404/1‑148 (MQ=255)
                                                                                                                      aCGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGTGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAAGCt           <  3:43404/148‑1 (MQ=255)
                                                                                                                      aCGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt           >  4:43402/1‑148 (MQ=255)
                                                                                                                      aCGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCt           <  3:43402/148‑1 (MQ=255)
                                                                                                                        gATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTg         >  5:142486/1‑148 (MQ=255)
                                                                                                                              tatGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGt  <  4:84176/149‑1 (MQ=255)
                                                                                                                              tatGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGt  <  6:50724/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                            
TCCTGTAGTGATGTTGCCAGAGTTACGCGCGCGTAAAGTAAGGTAAATAACTGAGTGGTTATTTTAGGGATGTAAGCGGTCAGTTTTGCGGTTGAGCGGCAAAGGCGTCGTTTAAGTGACGATAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGT  >  NC_000913/4591364‑4591638

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