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

ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAA  >  NC_000913/357725‑358011
                                                                                                                                                         |                                                                                                                                           
aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                             <  1:933201/151‑1 (MQ=255)
 cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                            >  4:116608/1‑151 (MQ=255)
   gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                          <  2:667798/151‑1 (MQ=255)
    aaCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                                         <  2:790144/151‑1 (MQ=255)
       tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                     >  1:1122626/1‑151 (MQ=255)
       tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGGGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                      >  1:1049657/1‑150 (MQ=255)
             gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                               <  1:221564/151‑1 (MQ=255)
                   tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                          >  1:1145391/1‑150 (MQ=255)
                    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                        >  1:747380/1‑151 (MQ=255)
                    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                        <  2:898811/151‑1 (MQ=255)
                         ttgcccATCACGCAACTCCAGCCGCAGACCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                       <  2:262774/142‑1 (MQ=18)
                             gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                               >  4:432392/1‑151 (MQ=255)
                                 cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                           >  2:362887/1‑151 (MQ=255)
                                 cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                           >  3:335910/1‑151 (MQ=255)
                                      atttCCAGCCGAAGCCTTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCAAACACCCGCGTTTTTATCCAACCCCTTGTGAAAAAATCCCTCGGCGGCGGCTGTCTGCCACTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                      <  1:1091573/149‑1 (MQ=255)
                                        tcccACCCGCAGCCCTGTGCCAGCGGCTACCGTGATTCATTCCCGCCAACACCCGCGAATTTACCCAACGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                       <  4:395173/146‑1 (MQ=18)
                                         cccACCCGCAGCCCTTTGCAGGGGGCTCCGGTGATTCTTTCCCCCCAACACCCGCGCTTTTCCCCAACCCCTTGTCAAAAAATCCCTCCGGGCCGGCTGCCTGCCATTTTTTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                   <  1:1038691/150‑1 (MQ=255)
                                              ccGCAGCCTTGTGTCAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTAAAAAAATGCCTCCGGGGCGGTTGTCTCCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                <  1:619855/149‑1 (MQ=255)
                                              ccGCAGCCTTGTGCCGGCGCCTCCCGTATTTCTTTCCCCCCACCACCCCCGCTTTTACCCACCCCCATTAAAAAAAAGCCTGCCGCGGCGGTTTCCCGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                               <  2:1271924/150‑1 (MQ=255)
                                              ccGCAGCCCTGTGTCGGCGGCTCCCGTTTTTCTTTCCCCCAAAAACCCGCGCATTTACCCAACGCCTTGTGAAAAAATCCCTCCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                <  4:203725/149‑1 (MQ=255)
                                              ccGCAGACCTGTGTCGGCGCCTACCGTTATTCATTCCCCCAAACACCCGCGCTTTTACCCACCCCCATGTCAAAAATTGCTTCCGCGGCGGCTGTCGGCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                <  1:300429/149‑1 (MQ=255)
                                               cGCAGCCGTGTGTCGGCGTTTACCGTGTTTCTTTCCCCCCAAAACCCGCGCATTTACCCAACCCAATGTGAAAAAATGCCTCCGGGGGGGCTGTCTGCCCGTTTTTTTCCCCCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                 <  1:665663/147‑1 (MQ=255)
                                                      cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggc                                                                                            >  1:874963/1‑143 (MQ=255)
                                                        gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaa                                                                                     >  1:116974/1‑150 (MQ=255)
                                                         tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGaaa                                                                                     >  1:398091/1‑149 (MQ=255)
                                                               cGGCTCCCGTGTTTCATTCCCGCCAACACCCGCGCATTTACCCACCCCCATGTGCAAAAATGCTTTCGGGGCGGCTTTCGCCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                               <  4:303011/133‑1 (MQ=255)
                                                                       gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattttaacccc                                                                       >  1:1234626/1‑135 (MQ=25)
                                                                         gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGgaaagaaaaaaattt                                                                              >  1:570421/1‑133 (MQ=25)
                                                                         gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGaaaaaaattttatc                                                                          >  4:179866/1‑133 (MQ=25)
                                                                           ttCTTTCCCCCCAACACCCGCGCTTTTCCCCACCCCATGGGGAAAAAAGCCTTTGGGGCCGTCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                          <  1:62727/142‑1 (MQ=255)
                                                                               ttccccccAACAACCGCGCTTTTACCCACCCCCATGTAAAAAAAGCCCTCCGGGGGGGCTGTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                               <  1:206885/143‑1 (MQ=255)
                                                                                  ccccaaacaacGGCGCATTTATCCACCGCCATGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                                                             <  1:457607/143‑1 (MQ=255)
                                                                                          cacccGCGTGTTTCTCCAACGCCATCCGCAAAAAGTCCTAGGCGGCGGCGGTCCGCCGTCTTTTTTTCCACCACAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                   <  1:859567/147‑1 (MQ=255)
                                                                                             ccccctatatccccACCGCTATGTACAAAATTGCTTCCGGGGCGGCTGTCCGCAGGTTATTTGCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                   <  1:1195709/135‑1 (MQ=255)
                                                                                                    tttATCAACCGCCATGTAAAAAAATGCCTGCGGGGCGCCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                            <  1:95084/146‑1 (MQ=255)
                                                                                                          cAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAATTTTTAACCCTCATGTTGTGTGGGGATCGGTGCCGGGGAAAAtgt                                  >  4:90960/1‑150 (MQ=25)
                                                                                                             ccccttgtgAAAAATTCCCTCCGCGGCGGTGTTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                  <  2:130410/142‑1 (MQ=255)
                                                                                                                      aaaaaaTGCCTGCGGGGCGGCTGTCGCCAATTTATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                     <  1:764244/150‑1 (MQ=255)
                                                                                                                      aaaaaaTCCCTCCGCGGCGGTTTTCTCCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                        <  2:901297/147‑1 (MQ=255)
                                                                                                                       aaaaaTCCCTCCGCGGCGCCTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                      <  2:1043787/149‑1 (MQ=255)
                                                                                                                                  tCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGGAAAAACTTTTAACCCGTCCTGGTGTGTGGGCATCGGCGGCGGTGAAAAAAGTGGAAAgggggggggggg                 >  2:461636/1‑136 (MQ=255)
                                                                                                                                  tCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAGATTTTAATCCCTCCTGCTGTTTGGGGAACGCTGCTGGGTAAAAACTTGGGAGGGGAAgggggg                  >  2:465247/1‑139 (MQ=255)
                                                                                                                                    gcggcgGTTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa           <  1:34922/148‑1 (MQ=255)
                                                                                                                                    gcggcgGCTCTCGCCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAg          <  2:398024/149‑1 (MQ=255)
                                                                                                                                      ggcggcTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGGGGCGAATCAGCTc       <  1:417086/150‑1 (MQ=255)
                                                                                                                                       gggggcTTTCTGCCGCTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc            <  3:93381/142‑1 (MQ=255)
                                                                                                                                          ggctgccgcccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg    <  1:476685/141‑1 (MQ=255)
                                                                                                                                          ggcTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa     <  1:73988/149‑1 (MQ=255)
                                                                                                                                          gcctgccgcccAGTTTTTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt             <  3:66087/132‑1 (MQ=255)
                                                                                                                                                cggCCAGCTTTTTTCCCGCACAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa  <  2:1056695/144‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                           
ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAA  >  NC_000913/357725‑358011

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