Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA42.0% 74.5 / 27.2 50L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (23/6);  new base A (0/21);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 5.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357741‑358034
                                                                                                                                                        |                                                                                                                                                   
tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                     >  2:318637/1‑150 (MQ=255)
  aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                  <  1:968790/151‑1 (MQ=255)
   tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                                  >  1:127642/1‑150 (MQ=255)
    aagcgcttcgcATCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                <  1:724464/143‑1 (MQ=255)
    aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                >  2:1690216/1‑151 (MQ=255)
      gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                               <  1:1409167/150‑1 (MQ=255)
      gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                <  2:1250385/149‑1 (MQ=255)
       cGTAGCGCATCAGGCCATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                             >  2:1053502/1‑151 (MQ=255)
                    gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                >  1:268393/1‑151 (MQ=255)
                       aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                             >  2:1482015/1‑151 (MQ=255)
                        ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                            >  1:502851/1‑151 (MQ=255)
                           cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                         >  1:342861/1‑151 (MQ=255)
                              cggaagtccTGTGCCGGGGGCTCCCTTGATTCTTTCCCGCAAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAAGCCCTCGGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                     <  2:1005859/144‑1 (MQ=255)
                              ccGCAGCCTTGTGCCGGCGCCTCCCGTGTTTCTTTCCCCCCAACACCCGCCCTTTTACCCACCGCAATTTGAAAAATTCCTTTCGCGGCGGCTGTCTCCCATTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                         <  2:340143/148‑1 (MQ=255)
                                 cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                     >  2:234779/1‑148 (MQ=255)
                                  aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                    >  1:828424/1‑147 (MQ=255)
                                    cctggtgtCGGCGCCTCCCGTGTTTCATTCCCCCAAACACCCGCGCTTTTCCCCACCCCCATGTCAAAAAATCCCTCCGCGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                                  <  2:58109/144‑1 (MQ=255)
                                       tgtgTCAGCGGCTCCGGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCAACCCCTTGTAAAAAAATCCCTTCGCGGCGGCTGTCTGCAAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                              <  2:858750/149‑1 (MQ=255)
                                       tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggggggg                                                                                                               >  1:48924/1‑142 (MQ=255)
                                        gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAGGAGaaa                                                                                                            >  2:1353738/1‑150 (MQ=255)
                                         tgtCGGCGGCTCCGTTAATTCTTTCCCCCAAACACCCGCGTATTTATCAACCGCAATGTGAAAAAATGCCTGGGCGGCGGTTTTCTGCAAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag                                                                                                               <  2:1546326/146‑1 (MQ=255)
                                         tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcagaaaa                                                                                                             >  2:1329351/1‑140 (MQ=255)
                                            cgggggCTCCCGTTTTTCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGTTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                                 <  2:636280/137‑1 (MQ=255)
                                            cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCACgagagaaa                                                                                                          >  2:1504475/1‑144 (MQ=255)
                                                ggctccggtttttCATTCCCCCCAACACCCGCGTATTTACCCACCGCCATGTAAAAAAATGCCTTCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                                   <  2:467460/140‑1 (MQ=255)
                                                  cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaatttt                                                                                                    >  1:1000947/1‑140 (MQ=255)
                                                     ccGTTATTCTTTCCCCCCACCACCCGCGCATTTATCCACCGCCTGTTAAAAAACTCCCTCCGCGGCGGCTGTCTCCCGGTTATTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                               <  2:1408215/150‑1 (MQ=255)
                                                       gTGATTCATTTCCGCCAAAAACCGCGCATTTATCCAACCCCATGTGCAAAAATCCCTACGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc                                                                                            <  2:1176847/151‑1 (MQ=255)
                                                               ttccccccAACACCCGCCCTTTTACCCACCGCCATGTAAAAAATTGCCTCCGGGGCGGTTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                        <  2:98028/141‑1 (MQ=255)
                                                                  cccccAACACCCGCATTTTTACCCACCGCTATGGGAAAAATTGCCTGCGGGGCGGCTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                                 <  2:751113/148‑1 (MQ=255)
                                                                      caacaaCCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACttttaaaccccc                                                                                            >  2:994518/1‑127 (MQ=25)
                                                                            ccccgtttttACCCACCCCTATGTGAAAAAATCCTTCCGCGGCGGCTTTCCGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt                                                                             <  2:935348/138‑1 (MQ=255)
                                                                            ccGCGATTTTACCCACCGCCATGTGAAAAAATGCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                        <  2:1001024/150‑1 (MQ=255)
                                                                                 cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGGAAAAAATTTTAACCCCCCCTGTTGTGTggggatagccggtgg                                                                    >  1:42786/1‑136 (MQ=21)
                                                                                   tttcccaaccgcattgtgAAAAAATGCCTCGGGGGCGGCTTTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                    <  1:378935/132‑1 (MQ=255)
                                                                                   tttACCCAACCCATTGTAAAAAAGTGCTTTGGGGGCGCTTGTCGGCCAGTTATTTTCCCGCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                               <  1:220809/150‑1 (MQ=255)
                                                                                               cgaggagaaaaaaGCCTTCCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                                   <  2:895920/144‑1 (MQ=255)
                                                                                                     cAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAACACCCCCTGTTTTGTGGGAAATGCTGCTGGGAaaaa                                                            >  2:1533516/1‑135 (MQ=255)
                                                                                                     aaaaaaTGCCCCCGGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                                <  1:1401790/146‑1 (MQ=255)
                                                                                                        aaaTGCCTCGGGGGCGGCTTTTGCCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                         <  2:1532354/151‑1 (MQ=255)
                                                                                                          aTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGACGGAAAAAAACTTAATGCCGCCCTTTTCTGTGGGGAACGCTGCGGGTCGAAAATGTGGAAAg                                                   >  2:252001/1‑140 (MQ=255)
                                                                                                                ttCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACCTTAAAACCGCCCTGTTTTGTGGGGATCGCTGCTGGTGACAAAGTGGAAAGGGGAgggggggggaac                                   >  2:20743/1‑139 (MQ=255)
                                                                                                                      ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAACCGCCCTGTTTTGTGGGCATTGCTGCTGGGGAAAAAGTTGGAAGGGGAGGGGGGGGGAACAGCTc                              >  2:1398953/1‑150 (MQ=255)
                                                                                                                        cggcTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                            <  2:214765/150‑1 (MQ=255)
                                                                                                                        cggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAACCCCTCCTGTTGTGTGGGGATCACCGCTGGGGAAAAAGAGGGAAGGGGGGGGGGGGGGATCAGCTcc                             >  1:634014/1‑149 (MQ=255)
                                                                                                                              tctgcagttttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                            <  1:120386/133‑1 (MQ=255)
                                                                                                                                      cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGGAGGGAAAAAAGATTTAACCCCGCCAGTTTTTGGGGCGTTGCTGCTTGGGAAAAATGTGGAAGGGGAGGGGGGGGGAACGACTCCACAACCgcgcacaat              >  1:1214422/1‑145 (MQ=255)
                                                                                                                                        ttttcttccgcgcaaacaaATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt               <  2:261810/134‑1 (MQ=255)
                                                                                                                                         tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAGACCGCCCCTTTGTGTGGGGATCGCCGCTGGGGACAATGTGGGAAGGGAAGGGGGGGGGATCCACTCCCGaaaccccccaatt              >  1:601692/1‑136 (MQ=255)
                                                                                                                                         ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa            <  2:1639110/149‑1 (MQ=255)
                                                                                                                                         atttcccgccgataacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGTCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg             <  1:214386/136‑1 (MQ=255)
                                                                                                                                                  ccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc        <  2:783486/145‑1 (MQ=255)
                                                                                                                                                  ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa  <  2:867853/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357741‑358034

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