Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT49.2% 37.0 / 37.5 64N272K (AAT→AAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (24/8);  new base T (0/31);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 8.39e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357731‑358023
                                                                                                                                                       |                                                                                                                                                       
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑c                                                                                                                                                  <  1:413633/151‑1 (MQ=255)
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑c                                                                                                                                                  <  2:1037507/151‑1 (MQ=255)
 tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑c                                                                                                                                                  >  1:956012/1‑150 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑ccgc‑cg                                                                                                                                            >  1:1324960/1‑150 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑ccgc‑cg                                                                                                                                            >  2:735856/1‑150 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑cccc‑cg                                                                                                                                            >  1:1065075/1‑150 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGcaacaa                                                                                                                                      >  2:1553464/1‑151 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGcaacaa                                                                                                                                      >  1:569531/1‑151 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGAGATTCATTTCCGCCAACAACCGCGAATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCCGCTATTTTT‑CCGC‑CGcaacaa                                                                                                                                      >  2:578935/1‑151 (MQ=255)
               aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAAc                                                                                                                                   >  1:415924/1‑150 (MQ=255)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAAccg                                                                                                                                 >  1:775338/1‑151 (MQ=255)
                             ggCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCag                                                                                                                    >  1:505602/1‑151 (MQ=255)
                             ggCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCag                                                                                                                    >  1:195687/1‑151 (MQ=255)
                                attcccggCCGAAGCCCGGTGC‑AAGCGCCTCCCTTAATTCTTTCCCCCCAACACCCGCTCTTTTCCCC‑AA‑CCCCATTTG‑AAAAATG‑CCCT‑CGGCGGGGG‑CTGTCGGCCATTTTTTTTC‑CCGC‑CCCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  1:1327964/144‑1 (MQ=255)
                                attccccgCCGAAGACCTGTGT‑CGGGGCCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCC‑AA‑CCCATTGTGAAAAAAGG‑CCTC‑C‑GGGGGGG‑TTTTCTGCCAGCATTTTTC‑CCGC‑CCAAACAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  1:958434/144‑1 (MQ=255)
                                atccccACCCGCAGCCGTGTGC‑CGGCGCCTACGGTGTTTCATTCCCCCCAACACCCGCGATTTTACCC‑AC‑CCCATTGTCCAAAAATC‑CCCT‑C‑GGGGCGG‑CTGTCTGCAGGTTTTTTTC‑CCGC‑CGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  2:371554/147‑1 (MQ=255)
                                actcccggCCGCAGCCTTGTGC‑CGGCGGCTACGGTGATTCATCCCCCCCAACCCCCGCGCTTTCCCCC‑AA‑CGCCATGTGAAAAAAGG‑CCTC‑C‑GCGGCGG‑CTGTCGGCAGGTTTTTTTC‑CCGC‑CGAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  2:44042/144‑1 (MQ=255)
                                aaccccAGCCGAACCCTTGTGA‑CGGCGGCCACTGTGATACTTCCCCGCAAACAACCGCTCTTTTACCA‑AA‑CGCATTGTAAAAAAATC‑CCTC‑C‑GCGGCGG‑CTGCCTGCCGGTTTTTTTC‑CCGC‑CGAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  2:1090162/147‑1 (MQ=255)
                                aaTCCCAGCCGCAGA‑CTGGGTCCAGCGCCACCCGGAATTCATTCCCGCAAACACCCGCGCTTTTCCCA‑AC‑CGCCTGGTGCAAAAATC‑CCTC‑C‑GCGGGGG‑CTTTCGGCAGGTTTTTTTC‑CCGC‑CGAAAAAAACCCGCCGTTCTCTCCagtag                                                                                                                 <  2:368585/151‑1 (MQ=255)
                                 aTTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                >  1:740216/1‑151 (MQ=255)
                                 aTTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                >  1:1623062/1‑151 (MQ=255)
                                        ccGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                         >  2:281465/1‑151 (MQ=255)
                                             gACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcagg                                                                                                    >  1:1106021/1‑146 (MQ=255)
                                             gACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgaca                                                                                                    >  2:402903/1‑147 (MQ=255)
                                              aCCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggtcgag                                                                                                   >  2:802305/1‑145 (MQ=255)
                                               cttgtgT‑CGGCGGATCCCGTGATTCTTTCCCCCCAACACCCGCCCTTTTACCC‑AC‑CGCCAGGTGAAAAAATG‑CCTT‑C‑GGGGGGG‑CTTTCTGCAAGTTTTTTTC‑CCCC‑CGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a                                                                                                 <  1:509576/149‑1 (MQ=255)
                                               ccTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaggga                                                                                                  >  1:21410/1‑144 (MQ=255)
                                                  gtgt‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaggaaa                                                                                               >  2:355973/1‑141 (MQ=255)
                                                  gtgt‑CAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaggggaaa                                                                                               >  1:452466/1‑141 (MQ=18)
                                                      ccggCGTCTACGGTATTTCATTCCCCCCAACACCCGGCCATTACCCC‑AC‑CGCGATG‑GCAAAAAAG‑GCTT‑G‑GGGGGGGTGGGTCGACCATCTGTTTTC‑CCCC‑CGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg                                                                                            <  1:337003/147‑1 (MQ=255)
                                                              tcccgtttttCTTTCCCCCCAACACCCGCG‑ATTTTTCC‑AA‑CCGCAAGGTAAAAAATG‑GCTTCC‑GGGGGGG‑CTGCCTCCCAGTTTTTTTC‑CCGC‑GGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                   <  1:411282/143‑1 (MQ=255)
                                                              tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑GAAAAACTTTaaaa                                                                                    >  1:1233587/1‑148 (MQ=25)
                                                              ccctgtgtTTCATTCCCCCAAACACCCG‑GCTTTTTACC‑AACCCCCGTGTGAAAAAATC‑CTTC‑C‑GGGGGGG‑CGGCCTCCCGGTTTTTTTC‑CCGC‑CAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                   <  2:980553/147‑1 (MQ=255)
                                                              cccggtatttCATCCCCGCAAACACCCGCGTTTTTCCCC‑AC‑CCCCTGGAGAAAAATTC‑CCTC‑C‑GCGGCGG‑TTGTCGCCCGGTTTTTTTC‑CCGC‑CAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                   <  2:602424/143‑1 (MQ=255)
                                                               aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAa‑aaaaaacattttacc                                                                                   >  1:1366820/1‑138 (MQ=25)
                                                                     aTTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTAAGAG‑AAACAACTTTaacaccccgcgt                                                                            >  2:1404033/1‑141 (MQ=25)
                                                                      ttcttcccccacaacaaCCGC‑TATTTAACCAAA‑CCCCGTGGGCAAAACTG‑CGTG‑G‑GGGGGGG‑CTGTCGGCAGGTTTTTTTC‑CCGC‑CAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgt                                                                           <  2:863313/140‑1 (MQ=255)
                                                                       tctttccccccaccaCCCGCCCATTTACCC‑AC‑CCCCTTGTGCAAAAATC‑CTTC‑C‑GCGGCGG‑CTGTCCCCCAGCTTTTTTC‑CCGC‑CCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgtt                                                                          <  1:763080/142‑1 (MQ=255)
                                                                        ttttccccccAACCCCCGCGCTTTTCCCC‑AA‑CCCCATGTAAAAAAATG‑CCTT‑C‑GGGGGGG‑CTGTCTGCCAGTTTTTTTC‑CCCC‑CGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                         <  1:347240/143‑1 (MQ=255)
                                                                        ctttccccccaccaCCCGCGCTTTTACCA‑AC‑CGCCATGTGAAAAAATC‑CCTC‑C‑GCGGCGG‑CTGTCTGCGAGTTATTTTC‑CCGC‑CGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                         <  2:191470/143‑1 (MQ=255)
                                                                        ctttccccccAACACCCGCGCTTTTCCCC‑AA‑CGCCATGTAAAAAAATG‑CCTT‑C‑GCGGGGG‑CTTCCTCCCAGTTTTTTTC‑CCGC‑CGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                         <  1:199571/143‑1 (MQ=255)
                                                                        catcccccccAACACCCGCGCTTTTACCC‑AC‑CGCCAGGTAAAAAAATC‑CCTG‑C‑GGGGCGG‑CTGTCTCCCAGTTTTTTTC‑CCGC‑CGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                         <  2:78311/143‑1 (MQ=255)
                                                                             ccaccAACACCCGCGCTTTCACCC‑AA‑CCCCATGTAAAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTC‑CCGC‑CACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                    <  2:1529981/148‑1 (MQ=255)
                                                                              ccccaccaCCCGCGCTTTCCCCC‑AC‑CCCCTTGTGAAAAAATG‑CCTT‑C‑GGGGCGG‑CTGTCTGCCAGCTTTTTTC‑CCGC‑CACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                   <  2:499741/149‑1 (MQ=255)
                                                                                              tttACCC‑AC‑CGCTATGAAAAAAATGG‑CTTT‑C‑GGGGCGG‑CTGTCGGCCTGCTTTTTTC‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  1:176602/151‑1 (MQ=255)
                                                                                              tttACCC‑AA‑CCCCTTGT‑AAAAAAAG‑CCTTCC‑GGGGCGC‑CTGTCTGCCAGTTTTTTTC‑CCGC‑CACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  2:597374/151‑1 (MQ=255)
                                                                                               tttacc‑AA‑CCCCATTGGAAAAAAAG‑CCTTGC‑GCGGCGG‑TTGGATGCCAG‑TTTTTTC‑CCGCGCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  2:1336393/147‑1 (MQ=255)
                                                                                               ttATCC‑AA‑CCCCCTGGGAAAAAAGGCCCTT‑C‑GCGGGGC‑CTGTCGCCGAGTTTTTTTC‑CCCC‑CGCAAAAAAACCGCCGTTCTCTCCAGTAATGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  1:1366406/151‑1 (MQ=255)
                                                                                                  cccccaacccatggtgCAAAACTC‑CTTC‑C‑GCGGCGC‑CTTTCTGCTAGTTTTTTTC‑CCGC‑CAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                 <  2:727415/138‑1 (MQ=255)
                                                                                                    cccaccccctgggaaaaaaaGC‑CCTC‑C‑GGGGCGG‑TTTTC‑GCCCGCTTTTTTTCCCGC‑CGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                               <  1:706090/136‑1 (MQ=255)
                                                                                                             ttgtcaaaaaaCC‑CCTC‑C‑GCGGCGG‑CTTTCTGCCATTTTTTTTC‑CCGC‑CGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                      <  2:1206674/145‑1 (MQ=255)
                                                                                                              tgtcaaaaagTG‑CTTC‑C‑GGGGGGG‑TTTTCGGCCAGTTTTTTTC‑CCGC‑CCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                     <  1:1598891/146‑1 (MQ=255)
                                                                                                              tgcgaaaaagTC‑CCTC‑C‑GCGGCGG‑CTGTCTGCCAGCTTTTTTC‑CCGC‑CGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                     <  2:166928/146‑1 (MQ=255)
                                                                                                              tgaacAAAAGTC‑CCCT‑C‑GCGGCGG‑CTGTCTGCCAGTTTTTTTC‑CCGC‑CGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                     <  2:109731/147‑1 (MQ=255)
                                                                                                               gtgCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAACTTTAACCCCG‑CCTGTTTTGTGGGCGACGCCGCCGGTTAAAAAGTGGGAAGGGag                                    >  2:1399586/1‑150 (MQ=255)
                                                                                                                   aaaaaTC‑CCTT‑C‑GGGGGGG‑TTTTCGCCCAGTTTTTTTC‑CCGC‑CCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                <  1:677849/151‑1 (MQ=255)
                                                                                                                                          ttgtctccaatttATTTTC‑CCGC‑CGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  1:357415/139‑1 (MQ=255)
                                                                                                                                          ttgcctgcagtttATTTTC‑CCCC‑CAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  2:441985/139‑1 (MQ=255)
                                                                                                                                          ttGTCTCCAAGTTATTTTC‑CCGC‑CGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  1:1270052/150‑1 (MQ=255)
                                                                                                                                          ctgcctgcGAGTTTTTTTC‑CCGC‑CGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  2:1403422/151‑1 (MQ=255)
                                                                                                                                          cTGTCTGCAGGTTTTTTTC‑CCGC‑CGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  2:1105498/151‑1 (MQ=255)
                                                                                                                                          cTGTCTCCCATTTTTTTTC‑CCCC‑CGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCGAGAAc             <  1:1053335/151‑1 (MQ=255)
                                                                                                                                          cTGTCGGCCAGTTTTTTTC‑CCGC‑CGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  1:844665/151‑1 (MQ=255)
                                                                                                                                                     cTATTTTC‑CCGC‑CGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg  <  2:1077036/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357731‑358023

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