Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,427,8750CT25.0% 110.9 / 29.2 52R284R (CGG→CGA) insH‑5Rac prophage; IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base C (21/18);  new base T (5/8);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.12e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTACTGAGATCTCTCCCACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAG  >  NC_000913/1427746‑1428013
                                                                                                                                 |                                                                                                                                          
ttAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTtgatga                                                                                                                                  <  2:328901/140‑1 (MQ=9)
ttACTGAGATCTCTCCCACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTtgatga                                                                                                                                  <  2:684333/140‑1 (MQ=255)
       gaTCTCTCCCACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGaaa                                                                                                                           >  2:93593/1‑140 (MQ=21)
       gaTCTCTCCCACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGaaa                                                                                                                           >  1:330560/1‑140 (MQ=25)
       gaTCTCTACCACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGaaa                                                                                                                           >  2:731053/1‑140 (MQ=255)
               ccACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCt                                                                                                                   <  1:1625/140‑1 (MQ=2)
               ccACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCt                                                                                                                   <  2:330560/140‑1 (MQ=2)
               ccACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCt                                                                                                                   <  1:88126/140‑1 (MQ=2)
               ccACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCt                                                                                                                   <  1:41734/140‑1 (MQ=17)
                cACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCt                                                                                                                   <  2:97321/139‑1 (MQ=255)
                cACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCt                                                                                                                   >  1:97321/1‑139 (MQ=255)
                cACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTc                                                                                                                  <  1:731053/140‑1 (MQ=255)
                  cTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCa                                                                                                                >  1:432241/1‑140 (MQ=18)
                  cTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCa                                                                                                                >  2:231544/1‑140 (MQ=255)
                  cTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCa                                                                                                                >  2:618321/1‑140 (MQ=255)
                    gACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCAcc                                                                                                              <  1:750312/140‑1 (MQ=2)
                     aCGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCt                                                                                                                   <  2:337022/134‑1 (MQ=255)
                     aCGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCt                                                                                                                   >  1:337022/1‑134 (MQ=255)
                     aCGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCAccc                                                                                                             <  1:231544/140‑1 (MQ=255)
                      cGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTc                                                                                                                  >  1:260967/1‑134 (MQ=255)
                      cGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTc                                                                                                                  <  2:260967/134‑1 (MQ=255)
                       gTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCACCCTg                                                                                                           >  1:328272/1‑140 (MQ=18)
                       gTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTg                                                                                                           >  2:103655/1‑140 (MQ=255)
                         aTCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACCGTCGCTTGATGATGCGAAATGGGTGCt                                                                                                                   >  1:562447/1‑130 (MQ=12)
                         aTCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACCGTCGCTTGATGATGCGAAATGGGTGCt                                                                                                                   <  2:562447/130‑1 (MQ=12)
                          tCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCAcc                                                                                                              <  2:403229/134‑1 (MQ=2)
                          tCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGTCGCTTGATGATGCGAAATGGGTGCTCCAcc                                                                                                              >  1:403229/1‑134 (MQ=18)
                               tGGTCCACCCGAAACAGGTTGGCCAGGGCGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCACCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGTAACGGGTGCTCCACCCTGGCACGGAt                                                                                                   >  2:491482/1‑140 (MQ=255)
                                    cACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTgg                                                                                              >  1:573480/1‑140 (MQ=255)
                                      cccGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCt                                                                                            >  2:22967/1‑140 (MQ=255)
                                             cAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGt                                                                                     <  1:22967/140‑1 (MQ=255)
                                                 ttGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTc                                                                                 >  1:764796/1‑140 (MQ=255)
                                                 ttGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTc                                                                                 >  1:419863/1‑140 (MQ=255)
                                                    gCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGAt                                                                              >  2:45800/1‑140 (MQ=255)
                                                             aaTAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGc                                                                     >  2:312197/1‑140 (MQ=40)
                                                              aTAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGcc                                                                    >  2:682328/1‑140 (MQ=40)
                                                                  cATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGttt                                                                <  1:103655/140‑1 (MQ=40)
                                                                   aTCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCTGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGtttt                                                               >  2:320441/1‑140 (MQ=255)
                                                                   aTCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGtttt                                                               >  1:325225/1‑140 (MQ=40)
                                                                    tCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTg                                                              >  1:532841/1‑140 (MQ=40)
                                                                            tGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTgcg                                                      <  2:325225/140‑1 (MQ=40)
                                                                            tGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTgcg                                                      <  1:312197/140‑1 (MQ=40)
                                                                                  tCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGAtt                                                >  1:263543/1‑140 (MQ=255)
                                                                                  tCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGAtt                                                <  1:45800/140‑1 (MQ=255)
                                                                                          cagcagCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCa                                        <  1:766959/140‑1 (MQ=255)
                                                                                          cagcagCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCa                                        <  1:618321/140‑1 (MQ=255)
                                                                                                  ccTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGttttt                                >  2:298505/1‑140 (MQ=255)
                                                                                                            ggCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCgg                      >  2:72834/1‑140 (MQ=255)
                                                                                                              cTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGttt                                  <  2:446032/126‑1 (MQ=255)
                                                                                                              cTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGttt                                  >  1:446032/1‑126 (MQ=255)
                                                                                                              cTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGa                    <  1:72834/140‑1 (MQ=255)
                                                                                                               tttCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGAc                   <  2:532841/140‑1 (MQ=255)
                                                                                                                              cTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGcc    <  2:573480/140‑1 (MQ=255)
                                                                                                                                gCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAg  <  1:320441/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                          
TTACTGAGATCTCTCCCACTGACGTATCATTTGGTCCACCCGAAACAGGTTGGCCAGGGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAGCCCCTTGTATCTGGCTTTCACGAAGCCGAACTGCCGCTTGATGATGCGAAACGGGTGCTCCACCCTGGCACGGATGCTGGCTTTCATGTATTCGATGTTGATGGCCGTTTTGTTCTTGCGCGGATTCTGCTTCAAGGTTTTTACCTTGCCGGGACGCTCGGCGATCAGCCAG  >  NC_000913/1427746‑1428013

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