Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC45.3% 74.2 / 22.1 64K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (27/8);  new base C (0/29);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 4.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357731‑358026
                                                                                                                                                          |                                                                                                                                                      
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                        >  1:951325/1‑150 (MQ=255)
     aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                 <  2:228031/151‑1 (MQ=255)
      tgCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                                >  3:229641/2‑151 (MQ=255)
       gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                                >  1:514694/1‑150 (MQ=255)
          tgctAAGCGTAGCGCATCCGGCAATTCCAGACGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                                                                                              >  1:89881/4‑149 (MQ=255)
          tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                             >  4:296096/1‑150 (MQ=255)
            aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                           >  1:782945/1‑150 (MQ=255)
            aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                          <  3:358069/151‑1 (MQ=255)
              aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                         <  1:1443762/150‑1 (MQ=255)
              aaGCGTAGCGCATAAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACCACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGCCGGCAGTCCTCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                         >  1:431767/1‑150 (MQ=18)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                       >  1:287807/1‑150 (MQ=255)
                 cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                     <  1:622152/151‑1 (MQ=255)
                    aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                   >  2:1414879/1‑150 (MQ=255)
                       gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                               >  2:1487317/1‑151 (MQ=255)
                       gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                >  1:764947/1‑150 (MQ=255)
                           cagccattcccAGCCGGTGACTTGTGCCAGCGGCTACCGTGATTCATTCCCGCCACCACCGCCGCATTTACCCACCCCCATGGGCAAAAAGGCCTTCGCGGCGGTTTTCTGCCAGCTTTTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                <  1:1625269/137‑1 (MQ=12)
                               cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                        >  2:51741/1‑150 (MQ=255)
                               acaTTCCACCCGCACACCTGTTCCAGTGTCTCCAGTTTTTAATCTCCGCAAACACCGGCGCATTTATCCACCCCATGTGCCAAAAGTGCCTTGCGGGCGCCTCTCTGCCAATATTTTTTCCCCCAAAAGAAAACCGCCGTTCTCTCCagt                                                                                                                        <  2:1740576/148‑1 (MQ=255)
                                attcccaccggCAGTCGTGTGTGGGCGGCCCCCTTTATTCTTTCCCCCAAACACCCGCCATTTCACCCACCGCAATGTGAAAAATGCCCTCCGGGGGGGTTGTCGCCCAGTTTTTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCag                                                                                                                         <  1:1591554/138‑1 (MQ=255)
                                aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                        >  1:1354919/1‑149 (MQ=255)
                                    ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                  >  4:391851/1‑151 (MQ=255)
                                    ccACCCGCAGCCTTGTGTCGGGGGCTCCGGTGTTTCTTTCCCCCCAACACCCGCGATTTTACCCACCCCCTTGTGAAAAAATGCCTCCGGGGCGGCTTTCTCCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                  <  2:964332/151‑1 (MQ=255)
                                      aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                >  2:1468973/1‑151 (MQ=255)
                                           cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggcg                                                                                                              >  2:1322673/1‑146 (MQ=255)
                                            agagttgggtCAGGGCCTCCCGTGATTTATCTCCGCAAACACCCGGCCATTTACCCACCCCTATGTGCAAAAATGCTTTCGGGGGGGTTGTCCGCCAGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              <  4:76039/139‑1 (MQ=255)
                                            accccTGTGTCAGCGGCTCCCGTGTTTCTTCCCCGCAAACCCCCGCGCATTTACCCACCGCCATGTAAAAAAACGCCTCCGCGGCGCCTGTCGCCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                            <  1:1313644/146‑1 (MQ=255)
                                            aGACCTGTGCCGGCGGCTCCGGTGATTCATTCCCCCCAACACCCGCCTTTTTACCCAACCCCAGGTGAAAAAAGCCTTTGGGGGGGGTTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                           <  2:726091/150‑1 (MQ=255)
                                             gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                              >  3:141947/1‑146 (MQ=255)
                                                  gtgtCAGCGGTTACCGTGTTCCATTCCCCCCAACCGCCGCGCTTTTACCCAACCCAATGGCCAAAATTGCCTCCGCGGCGCCTGCCTGCCAGTTTTTTTTCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                     <  2:1229653/149‑1 (MQ=255)
                                                   tgtCGGCGCCTACCGTATTTCTTTCCCCCCAACACCCCCCCTTTTACCCACCCCCTTGTGAAAAATGCCTTCGGGGCCGGTTTTCTCCCGTTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                   <  4:592157/150‑1 (MQ=255)
                                                     tCAGCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTTCGGGGCGGTTTTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                 <  3:442843/150‑1 (MQ=255)
                                                      cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcggggaaaa                                                                                                   >  1:719173/1‑137 (MQ=255)
                                                       aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggtagtgagacaac                                                                                                 >  1:1375636/1‑136 (MQ=255)
                                                        gCGGCTCCGGTTATTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTTCGGGGCGGTTTTTTCCAGCTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                 <  2:518525/147‑1 (MQ=255)
                                                          ggCTCCGGTGATTCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCAATGTGCAAAAATGCTTCCGGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                             <  2:1729899/148‑1 (MQ=255)
                                                          ggCTCCCGTGATTCATTCCCCCAAACACCCGCGCATTTACCCAACGCCATGTAAAAAAATCCCTTCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                   <  2:1597961/143‑1 (MQ=255)
                                                           gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGACGTTCTCTCCAGTAGTGGCggtg                                                                                                              >  1:952332/1‑130 (MQ=255)
                                                             accggTGATTCTTTCCCGCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAATGGCCTTCGGGGGGGTTGTCGCCCAGCTTTTTTCCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                        <  2:1388773/146‑1 (MQ=255)
                                                               ccgtatttCATCCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCTGCCGGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                     <  4:50371/145‑1 (MQ=255)
                                                                cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAATAACTTTTAAccccc                                                                                    >  1:1201560/1‑149 (MQ=25)
                                                                     ttCTTTTCCGCAAACACCCGCGCTTTTACCCAACGCCTTGTGAAAAAATCCCTTGGCGGCGGTTGTTTGCCTTCTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                            <  3:506470/138‑1 (MQ=255)
                                                                     ttCTTTCCCCCCAACACCCGCCTATTTATCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                                 <  3:226223/149‑1 (MQ=255)
                                                                       catccccccaaacaccCGCACTTTTACCCACCGCCTGTTAAAAAAATCCTTGCGCGGCGGCTGTCGGCCATTTATTTTCCCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                                <  1:137982/138‑1 (MQ=255)
                                                                           tCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAgaatttaaaccgccc                                                                                   >  1:202671/1‑128 (MQ=25)
                                                                            ccccaaacaccCGCTCTTTTACCCACCCCAATGTGAAAAAAGGCCTTCGCGGGGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                        <  2:644180/146‑1 (MQ=255)
                                                                             ccccAACACCCGCGATTTTACCCACCGCCATGTGAAAAAATGCTTTCGCGGCGGTTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                       <  3:174882/149‑1 (MQ=255)
                                                                                      ccgctattttACCCACCCCCAGGTAAAAAAGTCCCTCCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                              <  2:1539586/144‑1 (MQ=255)
                                                                                          gcttttACCAAACCCCTTGTAAAAAATTCCCTGCGCGGGGGTGGTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                             <  2:580097/145‑1 (MQ=255)
                                                                                             tttATCCACCGCCATGTGCAAAAATCCCTCCGGGGCGCCTTTCTCCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                       <  2:456331/150‑1 (MQ=255)
                                                                                             tttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCTCCTCTTGTGTgg                                                                        >  2:1309664/1‑134 (MQ=25)
                                                                                             tttACCCAACGCCATGTGAAAAATTGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                       <  2:1041295/150‑1 (MQ=255)
                                                                                                 aTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAATTTAAAACCCGCCTTTTGTGTGGGGAACGCTGCTggggaaaa                                                    >  2:711177/1‑144 (MQ=21)
                                                                                                       cGCCTTGTAAAAAAGTGCCTCCGGGGCGGCTGTCTGCCGTCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                              <  4:382090/150‑1 (MQ=255)
                                                                                                         ccATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTAAAACCCCCCTGTTTTGTGGGGATCGCCGGTGGGGACaa                                                    >  2:1225633/1‑142 (MQ=25)
                                                                                                             gtgAAAAAATGCCTCGGGGGCGGCTGTCTGCCAGTTATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                        <  2:772815/150‑1 (MQ=255)
                                                                                                                  aaaaaTGCTTTCGCGGGGGGTCGCTCCCGGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                           <  1:1079752/143‑1 (MQ=255)
                                                                                                                  aaaaaGGCTTTCGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                    <  3:72326/150‑1 (MQ=255)
                                                                                                                      aTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAACATTTAAACCCTCCTGGTGTGTGGGGAGCGCTCCTGGGTACAAATGGGGAAAGGGagggggg                                  >  3:413335/1‑143 (MQ=255)
                                                                                                                                  ggcggcTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                     <  3:576118/151‑1 (MQ=255)
                                                                                                                                      ggcTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  2:1495079/149‑1 (MQ=255)
                                                                                                                                           tCTGCGAGTTTTTTCCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                 <  1:1642425/147‑1 (MQ=255)
                                                                                                                                             cccctgtttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc           <  1:1426048/142‑1 (MQ=255)
                                                                                                                                              gCCAGGTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  4:413561/146‑1 (MQ=255)
                                                                                                                                               ccattttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt          <  2:460283/143‑1 (MQ=255)
                                                                                                                                               ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAACTTTTAAACCCCCCTGTTGTGTGGGAATCGCCACTGGGAAAAAAGTGGGAAGGGGAGGGGGGGGGAACAACCCCAGaa                 >  2:1593069/1‑143 (MQ=255)
                                                                                                                                                cggtttttttcCCGCAGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc           <  1:616437/142‑1 (MQ=255)
                                                                                                                                                 aGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAACTTTAACCCCTCCTTTTGTTTTGGGATCCCCGCTGGGTACAAATTGGGAAGGGGAGAGCGGGGGAACAGcccccgca                 >  2:1691015/1‑137 (MQ=255)
                                                                                                                                                  gttATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt      <  2:362924/149‑1 (MQ=255)
                                                                                                                                                      tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTTGCAGAGAAAAAAAATTAAAACCCTCCTTTTTTTTTGGAATAGCTGCTGGTAAAAATATGGAAAAGGGAGGGGGGGGGAGCAGCCCCCCAACCGCGCCAAATGAGt  >  4:396701/1‑151 (MQ=255)
                                                                                                                                                      ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg   <  2:202702/150‑1 (MQ=255)
                                                                                                                                                       tttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  2:1538440/150‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357731‑358026

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