Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC41.8% 153.3 / 30.8 103K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (43/17);  new base C (0/43);  total (43/60)
Fisher's exact test for biased strand distribution p-value = 3.54e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357733‑358027
                                                                                                                                                      |                                                                                                                                                       
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                      >  4:873020/1‑151 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                       <  2:1100113/150‑1 (MQ=255)
 gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                     <  4:721136/151‑1 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                    <  2:351848/150‑1 (MQ=255)
      ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaa                                                                                                                                                <  4:77399/151‑1 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                              >  2:723870/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                             >  4:904540/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                            <  2:1324119/151‑1 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                            <  2:368107/151‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                           <  2:707722/150‑1 (MQ=255)
              gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                        >  4:761616/1‑151 (MQ=255)
              gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                        >  4:262137/1‑151 (MQ=255)
              gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                        >  2:991939/1‑151 (MQ=255)
               cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                       >  2:360030/1‑151 (MQ=255)
                gTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTGCCGCC‑AA‑CAACCGCGCATTTCTCCAACGCCATGTGCAAAAATGCCTT‑CGCAGCGGCTGTCTGCCAGCTGTTGTTTAGCCCCAACGCGAccgcc                                                                                                                                      >  4:1226736/1‑151 (MQ=14)
                  aGCGCATCAGGCAATTCCAGCCGCAGCCTTGTTTCAGCGGCACCAGTGATTCATTCCCGCC‑AA‑CACCCGCGCATTAATCCAACTCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                     <  4:217965/150‑1 (MQ=18)
                  aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                    >  2:1447584/1‑151 (MQ=255)
                    cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                  >  2:1422943/1‑151 (MQ=255)
                          aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCa                                                                                                                            >  1:457302/1‑151 (MQ=255)
                              atttccccccGCAGCCGTGTGTCGGCGCCTCCCGTGTTTCTTTCCCCCC‑AC‑CACCCGCGCTTTTACCCAACCCCTTGTGCAAAAATCCCTC‑CGCGGCGGCTGTCTGCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                        <  4:1367532/143‑1 (MQ=255)
                              atttCCAGCCGAAGACCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCC‑AA‑CACCCGCGCTTTTATCCACCGCCATGTGCAAAAATGCTTT‑CGGGGCGGCTGTCTGCCACCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                        <  1:525304/149‑1 (MQ=18)
                              aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                          >  4:730292/1‑149 (MQ=255)
                                  ccAGCCTCAGCCCTGTGTCAGGGGCTCCCGTGATTAATTCCCGCC‑AA‑CACCCGCGCATTTACCCAACCCCTTGTGAAAAAATCCCTC‑CGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                    <  4:29075/151‑1 (MQ=255)
                                  ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                    >  4:972568/1‑151 (MQ=255)
                                   cAGCCGCAGCCCTGTGCCAGCGGCTACCGTTATTCTTTCCCGCA‑AA‑CACCCGCGCATTTACCCAACCCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTTCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                   <  2:789862/151‑1 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CACCCGCGCATTTATCCACCGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTATTTTTCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                   <  2:1337820/151‑1 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                   >  2:1423815/1‑151 (MQ=255)
                                    aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  >  4:848652/1‑151 (MQ=255)
                                        gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg                                                                                                              >  4:1106068/1‑149 (MQ=255)
                                        gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg                                                                                                              >  2:212207/1‑149 (MQ=255)
                                        ccggACGTGTGAGAGGGGCGACTGTGATCCTTTCCCGCC‑AA‑CAACCGCGCATTAATCCAGCGCGATGTAAAAAAATGCCTT‑CCGGGCGGCTGTCTACCGGGTATTTGTCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                              <  2:1393002/148‑1 (MQ=16)
                                        ccagccttgtgTCGGCGGCTACCGTGTTTCTTTCCCCCC‑AA‑CACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTC‑CGCGGCGGCTGCCGGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                              <  4:1233297/144‑1 (MQ=255)
                                           gACTTGTGTCAGCGGCTCCCGT‑ATTCCATTCCCCC‑AACCACCCGCGCTTTTACCCACCCCCATGTGAAAAAATCCCTC‑CGCGGCGCCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                           <  4:1040781/151‑1 (MQ=255)
                                            actggtgtCACCGGCTCCCGTGATTCATTCCCCCA‑AA‑CACCCCCGCTTTTATCCACCGCCATGTGCAAAAATGCCTC‑CGCGGCGGCTGTCGGCAAGCTTTTTTCCCGCGACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g                                                                                                         <  3:20288/147‑1 (MQ=18)
                                             ccTGTGTCGGCGGCTCCCGTGTTTCATTTCCGCC‑AA‑CAACCGCGCATTTACCCACCGCCATGGGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAga‑g‑a                                                                                                       <  2:1423092/151‑1 (MQ=255)
                                              cTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaga                                                                                                         >  4:122787/1‑143 (MQ=255)
                                               tgtgTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcggg                                                                                                          >  4:1051733/1‑142 (MQ=255)
                                                gtgtCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaggaa                                                                                                       >  4:1216314/1‑141 (MQ=255)
                                                gtgtCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGccaaaaaaa                                                                                                       >  2:1080993/1‑142 (MQ=255)
                                                gtgtCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGACga‑g‑ataa                                                                                                    >  3:474650/1‑148 (MQ=255)
                                                  gtCAGCGGCTCCCGTGATTCATTCCCGCC‑AA‑CAACCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAg                                                                                                  <  2:982231/151‑1 (MQ=25)
                                                  gccggcgcCTCCCGTGTTTCTTTCCCGCC‑AA‑CACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCTTC‑CGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAg                                                                                                  <  4:331378/147‑1 (MQ=255)
                                                   tCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑gaaaaa                                                                                                  >  4:801666/1‑147 (MQ=255)
                                                       cGGCTACCGTGATTCATTTCCCCC‑AA‑CACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTa                                                                                             <  2:1233179/151‑1 (MQ=255)
                                                          cTACCGTGATTCATTCCCGCC‑AA‑CACCCGCGCATTTACCCAACGCCATGTGCAAAAATCCTTT‑CGCGGCGGCTGTCTGCCACTTATTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGc                                                                                          <  1:190382/151‑1 (MQ=25)
                                                           tACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGA‑G‑AAAAAACttttaacc                                                                                         >  2:1284696/1‑146 (MQ=255)
                                                            cccGTGTTTCATTTCCGCC‑AA‑CAACCGCGCTTTTACCAAACCCCATGTGCAAAAAGGCTTC‑CGCGGCGGCTGTCTCCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGc                                                                                          <  2:1104562/148‑1 (MQ=255)
                                                            aCCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑aaaaaaacttaaacc                                                                                         >  1:414803/1‑138 (MQ=255)
                                                             ccgtatttCATTCCCGCC‑AA‑CACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCTTC‑CGGGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCg                                                                                       <  4:1364152/145‑1 (MQ=255)
                                                               gTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAACttttaaaccccc                                                                                     >  2:628715/1‑142 (MQ=25)
                                                                 gATTCATTTCCGCC‑AA‑CACCCGCGCATTTATCCAACGCCATGTGCAAAAATGCTTT‑CGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAt                                                                                   <  2:912299/151‑1 (MQ=37)
                                                                   ttCATTCCCGCC‑AC‑CACCCGCGCTTTTACCCAACCCCATGTGAAAAAATGCTTT‑CGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtg                                                                                  <  4:439113/150‑1 (MQ=255)
                                                                     cATTTCCGCC‑AA‑CACCCGCCCATTTACCCAACGCCATGTGCAAAAATGCCTC‑CGCGGCGGCTTTCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgtt                                                                                <  2:1215652/150‑1 (MQ=38)
                                                                     cATTCCCCCC‑AA‑CACCCGCGCATTTACCCACCGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttg                                                                               <  1:358644/151‑1 (MQ=37)
                                                                      ctttCCCCCAAA‑CACCCGCTCATTTACCCACCGCCATGGCAAAAAATGCCTC‑CGCGGCGGTGGCCTGCAGGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtgttg                                                                               <  3:69406/150‑1 (MQ=255)
                                                                      aTTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGG‑AAAAAGATTTAAccccccctgtttt                                                                               >  4:722925/1‑142 (MQ=25)
                                                                        ttCCCCC‑AA‑CACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgtg                                                                             <  2:498897/150‑1 (MQ=255)
                                                                         cccccc‑AA‑CACCCCCGCTTTTTCCCACCCCCATGTGAAAAATTGCCTT‑CGCGCCGGCTGCCTCCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                           <  2:1417728/147‑1 (MQ=255)
                                                                          cccccaaa‑caccCGCCCTTTTACCCACCCCCATGTG‑AAAAAAGCCTTCCGGGGGGGCTTTCGGCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                           <  4:1429867/145‑1 (MQ=255)
                                                                          ccccc‑AA‑CACCCGCGCATTTACCCAACGCAATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                          <  1:551792/148‑1 (MQ=255)
                                                                          ccGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAAATTTTAAACCGCCATGTTTTTTgg                                                                          >  2:1289536/1‑151 (MQ=25)
                                                                             ccaaa‑caccCGCGCATTTACCCAACCCCATGTGAAAAAAGCCCTT‑CGCGGCGGCTGTCTGCCGGCTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                        <  2:896112/148‑1 (MQ=255)
                                                                               caaaaacccGCGAATTTACCCACCGCCATGTGAAAAAAGGCCTT‑CGGGGCGGCTGTCGGCCAGCTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                       <  1:224229/144‑1 (MQ=255)
                                                                                aa‑caccCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTATTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                       <  2:669551/103‑1 (MQ=255)
                                                                                aa‑caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                       >  4:669548/1‑103 (MQ=255)
                                                                                   cacccGCGCTTTTACCCACCGCCATGTCCAAAAAGCCCTC‑CGGGGCGGCTGCCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgc                                                                   <  4:1336974/148‑1 (MQ=255)
                                                                                        gctcATTTCCCCACCGCCTTGTAAAAAAATGCCTC‑CGGGGCGGTTCGCCCCCCCTTTTTTTCCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGAAATAGctgg‑tg                                                             <  4:830259/148‑1 (MQ=255)
                                                                                         cgcttttcaccAACGCCATGTCAAAAAATCCCTC‑CGCGGCGGCTGTCGGCCGGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑Tgg                                                            <  3:243782/142‑1 (MQ=255)
                                                                                          gcttttaccaaccGCCATGTAAAAAAATGCTTT‑CGCGGCGGCTGTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGt                                                           <  2:1207005/139‑1 (MQ=255)
                                                                                           tatttccccAGCGCTATGTGCAAAAATGCCTT‑CGGGGCGGTTGTTTACAAGATATTTTTCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGt                                                           <  2:1427528/143‑1 (MQ=255)
                                                                                           cATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAAAAAATTTATCCCGTCATATTTTGTGGGaaaggcttccggg                                                            >  2:77141/1‑139 (MQ=25)
                                                                                            aTTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGG‑G‑AAAAAAGTTTTAACCCCCCCTGTTGTGTgggggtcgcggccgggg                                                           >  4:1460577/1‑135 (MQ=25)
                                                                                             tttaccaacccccATGTGAAAAAATGCCTG‑CGCGGGGGCTTTCGGCCAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTaa                                                         <  4:546657/139‑1 (MQ=255)
                                                                                             tttACCCACCCCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAAc                                                        <  1:91023/151‑1 (MQ=255)
                                                                                               tacccacccccATGTGCAAAAAGGCCTC‑CGCGGCGGTTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACaa                                                      <  1:514636/142‑1 (MQ=255)
                                                                                                acccAACGCCATGTAAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGc                                                                                          <  3:436155/113‑1 (MQ=255)
                                                                                                 ccaaccGCAATGTGAAAAAAGGCTTC‑CGGGGCGGCTTTCTGCCGGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAAtg                                                    <  3:419756/146‑1 (MQ=255)
                                                                                                       cccATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGaaa                                             <  3:87145/150‑1 (MQ=255)
                                                                                                        ccATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑A‑AAAAAGATTTTAGCCCGGCACGTTGTGTGGGAAACGCTGCTTGGTACCaaaaggaa                                                >  2:232322/1‑143 (MQ=37)
                                                                                                          aTGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAggg                                          <  2:768918/151‑1 (MQ=255)
                                                                                                             tgCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAAATTTTAACCCCCCCCGTTGTGTGGGGATCGCTGC‑TGGTGACAATGT‑GGGAAgggg                                         >  4:1458743/1‑148 (MQ=255)
                                                                                                               cAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAA‑G‑AAAAAAATTTTAAACCCCCATGTTGTGTGGGAATCGCTGC‑TGGGAACAATGT‑Gggaagggaag                                       >  2:412348/1‑149 (MQ=255)
                                                                                                                 aaaaTGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGG‑AAAAAACTTTAA‑CCCCCCATGGTTTGTGGGGAACGCCGC‑TTGGAAACATTT‑GGGAAGGGAAGGgcg                                   >  4:1139846/1‑151 (MQ=255)
                                                                                                                      gcttt‑CGGGGCGGCTGTCTCCCAGCTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGa                               <  2:1226729/147‑1 (MQ=255)
                                                                                                                       ccTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGAA‑GAAAAAAACTTTAA‑AACGCCATGTTGTGTGGGGATAGCCGC‑TGGGAAAAAAGT‑GGGAAGGGGAggggggggga                              >  2:132701/1‑141 (MQ=255)
                                                                                                                           ccgcgGCGGTTGTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAg                          <  4:1423818/150‑1 (MQ=255)
                                                                                                                             gcgccgGCTGTCCCCCAGTTTTTTTCCCGCCCCACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCt                        <  2:306317/151‑1 (MQ=255)
                                                                                                                              ggggggcTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTc                       <  4:1271465/147‑1 (MQ=255)
                                                                                                                                ggggcTTTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCa                     <  4:1243140/149‑1 (MQ=255)
                                                                                                                                    cggtcggccactttttttcCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  4:498793/137‑1 (MQ=255)
                                                                                                                                    cTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑G‑AAAAAAATTTTAACCCCCCCTTTTTTGTGGGAATAGCTTC‑TGGTGACAATGTGGGAAAGGGAGGGGGGGGGAA‑CAGCTCCaaaa                  >  2:780670/1‑150 (MQ=255)
                                                                                                                                      gccTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                 <  2:1297782/147‑1 (MQ=255)
                                                                                                                                      gTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCAGA‑G‑AAAAAACTTTAAACCCCGCAAGTTGTGTGGGGATTGCCGC‑TGGTGAAAAAGT‑GGGAAGGGGAGGGCGGCGGAACAGCTCCAGAAccc               >  2:574983/1‑150 (MQ=255)
                                                                                                                                          gCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg             <  4:183887/149‑1 (MQ=255)
                                                                                                                                             agtttttttcCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt        <  2:801669/146‑1 (MQ=255)
                                                                                                                                             aGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAA‑G‑GAAAAACTTTTAAACCCCCATGTTGTGGGGGGAATGCCGG‑TGGTAACAATGT‑GGGAAGGGAAGGGCGGGGGAACAACTCCCGAAACgcgcccat        >  2:467379/1‑147 (MQ=255)
                                                                                                                                                ttttttcCCGCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg      <  4:1248303/148‑1 (MQ=255)
                                                                                                                                                 aTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGTGGCAGA‑G‑AAAAAAACTTAAAACCCTCCCGTTGTGTGGGGGTCCCTGC‑TGGGGACAAAAT‑GGGAAGGGGAAGGGGGCGGAGCACCTCCCGAACCCCGCCaaatta     >  4:171259/1‑146 (MQ=255)
                                                                                                                                                   ttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑G‑AAAAAAAATTAAACCCCTCCTGTTTTGTGGGAATTGCTGG‑TGGTGAAAAAGT‑GGGAAGGGGAGGGGGGGGGAACAACTCCCCAACCGCGCCAAATGAggt  >  2:1463101/1‑149 (MQ=255)
                                                                                                                                                   tttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt  <  1:178900/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATGT‑GGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357733‑358027

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: ‑