Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913272,3570GT69.2% 28.3 / 24.7 26intergenic (+167/‑223)insN/eyeACP4‑6 prophage; IS911A regulator fragment/CP4‑6 prophage; small RNA EyeA
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base T (8/10);  total (12/14)
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.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACAT  >  NC_000913/272279‑272433
                                                                              |                                                                            
ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCctct                                                                            >  2:57003/1‑81 (MQ=255)
ccAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCt                                                                            <  2:56964/81‑1 (MQ=255)
     aaTCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACaa                                                                       >  1:57011/1‑81 (MQ=255)
       tCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACg                                                                     >  1:57001/1‑81 (MQ=255)
         ccACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACgg                                                                    <  2:56965/80‑1 (MQ=255)
            ccGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGAt                                                                <  1:56966/81‑1 (MQ=255)
            ccGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGAt                                                                <  2:56967/81‑1 (MQ=255)
             cGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGAtt                                                               <  2:56968/81‑1 (MQ=255)
              gggTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTg                                                              <  2:56969/81‑1 (MQ=255)
                  gtgACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGc                                                          >  1:57006/1‑81 (MQ=255)
                    gACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGc                                                          >  2:56994/1‑79 (MQ=255)
                       ggTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGc                                                          >  1:57012/1‑76 (MQ=255)
                          ggAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCAc                                                  >  1:57025/1‑81 (MQ=255)
                           gagaACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACt                                                 <  1:56970/81‑1 (MQ=255)
                              aaCCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAg                                              <  2:56971/81‑1 (MQ=255)
                                ccTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGcc                                            <  1:56972/81‑1 (MQ=255)
                                       ggCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatga                                     <  2:56973/81‑1 (MQ=255)
                                       ggCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatga                                     <  2:56974/81‑1 (MQ=255)
                                        gCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAgatgat                                    <  2:56975/81‑1 (MQ=255)
                                          aGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATcc                                  >  1:57008/1‑81 (MQ=255)
                                          aGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATcc                                  >  2:57014/1‑81 (MQ=255)
                                                   cGGCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGtt                          <  1:56976/80‑1 (MQ=255)
                                                    ggCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGtt                          >  1:57019/1‑79 (MQ=255)
                                                    ggCATGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTaa                        <  1:56977/81‑1 (MQ=255)
                                                       aTGCCCGCAGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTa                     >  2:57013/1‑81 (MQ=255)
                                                               aGGGTGATGTAACCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAgg             >  2:57022/1‑81 (MQ=255)
                                                                          aCCCTCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACAt  >  1:57010/1‑81 (MQ=255)
                                                                              |                                                                            
CCAACAATCCCACCGGGTGTGACGGTGGAGAACCTGAGCGGCAGTGACCTGCGGCATGCCCGCAGGGTGATGTAACCCGCTGACAACGGGGATTGAGGCGAGATCACTAAGCCGAGATGATCCTCAAGGTTAAGTACTGAAAGGTTGAAGAACAT  >  NC_000913/272279‑272433

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