Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT40.5% 43.9 / 13.8 37N272K (AAT→AAAcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (17/5);  new base T (0/15);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 3.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357727‑358019
                                                                                                                                                       |                                                                                                                                                 
tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtt                                                                                                                                                 >  1:118940/1‑151 (MQ=255)
          ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                        <  1:243324/150‑1 (MQ=255)
            ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                      >  1:920564/1‑150 (MQ=255)
             ctaaaaaGCGTAGCGCATCAGGCTATTCCACCCGCAGTCCTGGCTCCGCGCCTTGCTTGATCCATTTCCGCCAACACCCGCGCATTGTTCCACCGCCATGTGCAAAAATGCCTTCGCGGCTGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                      <  1:1172270/144‑1 (MQ=17)
               gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                   <  2:659290/150‑1 (MQ=255)
                aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGAGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                                                                                     >  1:477322/1‑147 (MQ=255)
                  aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                               >  1:1100999/1‑151 (MQ=255)
                       tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                          >  2:31148/1‑151 (MQ=255)
                        aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                         >  1:919355/1‑151 (MQ=255)
                            cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                     >  1:725895/1‑151 (MQ=255)
                            cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                     >  2:783481/1‑151 (MQ=255)
                             aTCAGGCAATGCCAGCCGCAGACCTGAGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGACGTtct                                                                                                                      >  1:355597/1‑149 (MQ=255)
                                          aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                      >  1:190301/1‑151 (MQ=255)
                                            ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                    >  2:814014/1‑151 (MQ=255)
                                                  ccgtgtgCCGGCGTCTCCGGTGATTTTTTCCCCCAAACACCCGCCCTTTTACCCAACCCCATGGCAAAAAGTCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                <  2:1219728/146‑1 (MQ=255)
                                                   ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGCGGCggggg                                                                                                   >  1:5610/1‑144 (MQ=255)
                                                     tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggacagg                                                                                             >  1:693371/1‑142 (MQ=255)
                                                             ctgCTACGGTATTTCTTTCCCGCCAACACCCCCGCTTTTACCCAACCCCTGGTGAAAAAAGCCTTCCGCGGGGGCTTTCTCCAGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                          <  2:731707/142‑1 (MQ=255)
                                                                ctccggtgtTTCTTTCCCGCCAACACCCGCCCTTTTACCCACCGCCTTGTGAAAAAATCCCTCCGGGGCGCTTGTCTCCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                 <  2:836702/145‑1 (MQ=255)
                                                                  aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaaattaaaagc                                                                              >  2:585065/1‑138 (MQ=25)
                                                                      attttCATCCCCCCCACCACCCGCCCTTTCACCCACCGCATTGTCAAAAACGGCCTCCGGGGCGCCTGTCTCCCAGTTTTTTTCCCGCCACAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                <  1:247189/142‑1 (MQ=255)
                                                                         taaatcccccccAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATCCTTCGGGGGGGGCTGTCTGCCAGTTTTTTTTCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                         <  2:431267/139‑1 (MQ=255)
                                                                           cATTCCCCCCACCACCCGCGCATTTCTCCAACGCAATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                      <  1:1141385/150‑1 (MQ=255)
                                                                                       caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAATAAATTTTAACCCCGCCTGTTTTGTgg                                                                >  2:502570/1‑144 (MQ=25)
                                                                                                aTTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGAAAAAAGCTTTAACACCGCCTGTTTTGTgggggagtct                                                        >  2:39833/1‑135 (MQ=25)
                                                                                                      caacccccATGTAAAAAAATCCCTCCGGGGGGCCTGCCTCCCATTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                               <  1:918972/140‑1 (MQ=255)
                                                                                                            ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAATCCCCATGTTGTGTGGGAATAGCTGCTGGGaagcaagagggaa                                     >  2:491663/1‑139 (MQ=25)
                                                                                                            ccATGTAAAAAAATCCTTCCGGGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                    <  2:625501/151‑1 (MQ=255)
                                                                                                                        attcccTCCGGGGCGGCGTTCGGCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                         <  1:308755/147‑1 (MQ=255)
                                                                                                                         acgCCTTCGCGGCGCGTGTCGGCCATCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                        <  2:318772/149‑1 (MQ=255)
                                                                                                                                   ggggcggcTTTCCGCCGGCTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                        <  2:775044/141‑1 (MQ=255)
                                                                                                                                     ggcgtcTGTCCGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                      <  1:1070772/143‑1 (MQ=255)
                                                                                                                                     ggcggtTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc               <  2:875240/150‑1 (MQ=255)
                                                                                                                                     ggcggtTGTCGCCCATTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                <  1:977023/149‑1 (MQ=255)
                                                                                                                                     ggcggcTGTCTGCCATTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc              <  2:919422/151‑1 (MQ=255)
                                                                                                                                     ggcggcTGTCTCCCTGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                <  1:473264/149‑1 (MQ=255)
                                                                                                                                                  cAGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa  <  1:498929/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
TGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTA  >  NC_000913/357727‑358019

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