Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913266,3550AG34.8% 32.1 / 16.4 23L67L (TTG→CTG) yafYCP4‑6 prophage; inner membrane lipoprotein YafY
Reads supporting (aligned to +/- strand):  ref base A (0/15);  new base G (3/5);  total (3/20)
Fisher's exact test for biased strand distribution p-value = 3.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCAT  >  NC_000913/266220‑266489
                                                                                                                                       |                                                                                                                                      
aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTg                                                                                                                                    <  1:603161/140‑1 (MQ=255)
 ttCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTgg                                                                                                                                   <  1:495171/140‑1 (MQ=255)
  tgcAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGAAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTggg                                                                                                                                  <  1:633015/138‑1 (MQ=255)
   ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGa                                                                                                                                 <  1:597766/140‑1 (MQ=255)
    cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGAc                                                                                                                                <  1:139760/140‑1 (MQ=255)
     aGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATGGGGCGACCGGCGGGGAAGTAAAAACTGCAAGGTgggggc                                                                                                                               >  2:272338/1‑137 (MQ=255)
       tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAACCATCGGGCGACTCACGGGTACGTAAGACCTGCAGGGt                                                                                                                                     <  2:547881/132‑1 (MQ=14)
       tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAACCATCGGGCGACTCACGGGTACGTAAGACCTGCAGGGt                                                                                                                                     >  1:547881/1‑132 (MQ=14)
        ttGATGTTATCGGCAATACACAGCGCATTCATCAATGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAACCATCGGGCGACTCACGAGTACGTAAGACCTGCAGGGTGAGACCGGa                                                                                                                            >  1:324562/1‑140 (MQ=14)
         tGATGTTATCGGCAATACACAGGGCATTCATCAGTGGCTCGACGGTTTCTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAACCATCGGGCGACTCACGGGTACGTAAGACCTGCAGGGTGAGACCGGAt                                                                                                                           <  2:255433/140‑1 (MQ=11)
         tGATGTTATCGGCAATACACAGCGCCTTCATCAGTGGTTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGAt                                                                                                                           <  1:330412/140‑1 (MQ=255)
            tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCCTCGGGCGAGTGGGGGGGACGTAAAACCAGCAAGAGGGGGGCGGATAag                                                                                                                        >  2:258737/1‑140 (MQ=255)
            tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAACCATCGGGCGACTCACGGGTACGTAAGACCTGCAGGGTGAGACCGGAtgag                                                                                                                        <  2:576283/140‑4 (MQ=11)
            tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAACCATCGGGCGACTCACGGGTACGTAAGACCTGCAGGGTGAGACCGGAt                                                                                                                           >  2:699369/1‑137 (MQ=14)
            tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAACCATCGGGCGACTCACGGGTACGTAAGACCTGCAGGGTGAGACCGGAt                                                                                                                           <  1:699369/137‑1 (MQ=11)
                aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGTCATCCTCGCAACCATCGGGCGACTCACGGGTACGTAAGACCTGCAGGGTGAGACCGGAtgagaagc                                                                                                                    <  1:391644/140‑8 (MQ=11)
                               cgcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGCCGGGTAAATAAGCAACCCGTGACCAtt                                                                                                     >  1:621007/1‑140 (MQ=255)
                                     catcaGTGGCTGGACGGTTTTTTCCCGCAAATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGc                                                                                               <  1:272338/140‑1 (MQ=255)
                                       tcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAAGGGGGGCCGGGTAAATAAGCAAAATGTGGCCAatctttcccg                                                                                             >  2:680165/1‑135 (MQ=255)
                                       tcaGTGGCTGGACGGATTTTTCGTGTATACCCACGAGTTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTg                                                                                             <  1:345956/140‑1 (MQ=21)
                                                ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGTCATGAAAACTGCCAGGGGGGGACCGAAAAGCAGTGGACCGTTGGCCATATTTCGCTGGcacagtcc                                                                                    >  2:193110/1‑134 (MQ=255)
                                                           tCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGAc                                                                         <  2:54250/140‑1 (MQ=255)
                                                             cTGCATATCCACGAGGTAGTAGGAAATGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACca                                                                       <  1:267591/140‑1 (MQ=255)
                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACAt                                                                    <  1:680165/140‑1 (MQ=255)
                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGATAAGTAGACAACCTGTGACCATTCATCGCTGACGACATCCTGGTGGGTGGCcccac                                                                    >  2:258783/1‑139 (MQ=255)
                                                                cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGCACGTCAAAACCGCAAAGAGGGGCCGGGTAAGTAGTGAAACCTCGACCATTTTTCGCCTCCATTATCCTGGTGGCTGACcaccc                                                                    >  2:604619/1‑138 (MQ=255)
                                                                             tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGGGCGGGTAAGTAATCAACCTGTGACCGTTCTTCGCAGACGTCAACCTGGGGGGTGGCGACATCACCGAAAAGCgg                                                       >  2:186718/1‑140 (MQ=255)
                                                                              agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGGAAAAACCGCAAAGGGGGGCCGGATAAAGAGTCAAACTGTGGCCAGTCTTCGCTTACAACATTCTTGTGGCTGGCCCCACCACAGAAACACggg                                                      >  2:237237/1‑139 (MQ=255)
                                                                              agtagtaACTGCCATCCTCGCAACCCTCGCGCGACTGTCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGt                                                      <  1:193110/140‑1 (MQ=255)
                                                                                tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTGCGGAAAACCCGCAAGGGGGGGCCGGGAAAAAAGGGAACCTTTGACCCTTGTTCCCTCACCCCATTCGGGGGGCTGTCCACATAACTGAAACGCGgggg                                                    >  2:231729/1‑140 (MQ=255)
                                                                                                              gACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt                      <  1:604619/140‑1 (MQ=255)
                                                                                                                      ggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcat              <  1:113184/140‑1 (MQ=255)
                                                                                                                                  cTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCAt  <  1:487351/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
ATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCAT  >  NC_000913/266220‑266489

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