Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	57.8%	30.8 / 87.0	64	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (24/3);  new base C (0/37);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 1.17e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-11

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357731‑358025                                                                                                                                                             |                                                                                                                                                      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                                                                                                                                               >  2:735856/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‑cccc‑cg                                                                                                                                               >  1:1065075/1‑150 (MQ=255)            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGc‑aacaa                                                                                                                                        >  2:1553464/1‑151 (MQ=255)            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGc‑aacaa                                                                                                                                        >  1:569531/1‑151 (MQ=255)            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGAGATTCATTTCCGCCAACAACCGCGAATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCCGCTATTTTT‑CCGC‑CGc‑aacaa                                                                                                                                        >  2:578935/1‑151 (MQ=255)                aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAAc                                                                                                                                     >  1:415924/1‑150 (MQ=255)                 gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAAccg                                                                                                                                   >  1:775338/1‑151 (MQ=255)                              ggCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCag                                                                                                                      >  1:505602/1‑151 (MQ=255)                              ggCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCag                                                                                                                      >  1:195687/1‑151 (MQ=255)                                 attcccggCCGAAGCCCGGTGC‑AAGCGCCTCCCTTAATTCTTTCCCCCCAACACCCGCTCTTTTCCCC‑AA‑CCCCATTTG‑AAAAATG‑CCCT‑CGGCGGGGG‑CTGTCGGCCATTTTTTTTC‑CCGC‑CCC‑AACAAAACCGCCGTTCTCTCCagtag                                                                                                                   <  1:1327964/144‑1 (MQ=255)                                 attccccgCCGAAGACCTGTGT‑CGGGGCCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCC‑AA‑CCCATTGTGAAAAAAGG‑CCTC‑C‑GGGGGGG‑TTTTCTGCCAGCATTTTTC‑CCGC‑CCA‑AACAAAACCGCCGTTCTCTCCagtag                                                                                                                   <  1:958434/144‑1 (MQ=255)                                 atccccACCCGCAGCCGTGTGC‑CGGCGCCTACGGTGTTTCATTCCCCCCAACACCCGCGATTTTACCC‑AC‑CCCATTGTCCAAAAATC‑CCCT‑C‑GGGGCGG‑CTGTCTGCAGGTTTTTTTC‑CCGC‑CGC‑AAAAAAACCGCCGTTCTCTCCagtag                                                                                                                   <  2:371554/147‑1 (MQ=255)                                 actcccggCCGCAGCCTTGTGC‑CGGCGGCTACGGTGATTCATCCCCCCCAACCCCCGCGCTTTCCCCC‑AA‑CGCCATGTGAAAAAAGG‑CCTC‑C‑GCGGCGG‑CTGTCGGCAGGTTTTTTTC‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCagtag                                                                                                                   <  2:44042/144‑1 (MQ=255)                                 aaccccAGCCGAACCCTTGTGA‑CGGCGGCCACTGTGATACTTCCCCGCAAACAACCGCTCTTTTACCA‑AA‑CGCATTGTAAAAAAATC‑CCTC‑C‑GCGGCGG‑CTGCCTGCCGGTTTTTTTC‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCagtag                                                                                                                   <  2:1090162/147‑1 (MQ=255)                                 aaTCCCAGCCGCAGA‑CTGGGTCCAGCGCCACCCGGAATTCATTCCCGCAAACACCCGCGCTTTTCCCA‑AC‑CGCCTGGTGCAAAAATC‑CCTC‑C‑GCGGGGG‑CTTTCGGCAGGTTTTTTTC‑CCGC‑CGA‑AAAAAACCCGCCGTTCTCTCCagtag                                                                                                                   <  2:368585/151‑1 (MQ=255)                                  aTTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCagtagt                                                                                                                  >  1:1623062/1‑151 (MQ=255)                                  aTTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCagtagt                                                                                                                  >  1:740216/1‑151 (MQ=255)                                         ccGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:281465/1‑151 (MQ=255)                                              gACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcagg                                                                                                      >  1:1106021/1‑146 (MQ=255)                                              gACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgaca                                                                                                      >  2:402903/1‑147 (MQ=255)                                               aCCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggtcgag                                                                                                     >  2:802305/1‑145 (MQ=255)                                                cttgtgT‑CGGCGGATCCCGTGATTCTTTCCCCCCAACACCCGCCCTTTTACCC‑AC‑CGCCAGGTGAAAAAATG‑CCTT‑C‑GGGGGGG‑CTTTCTGCAAGTTTTTTTC‑CCCC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgag‑a                                                                                                   <  1:509576/149‑1 (MQ=255)                                                ccTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaggga                                                                                                    >  1:21410/1‑144 (MQ=255)                                                   gtgt‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaggaaa                                                                                                 >  2:355973/1‑141 (MQ=255)                                                   gtgt‑CAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaggggaaa                                                                                                 >  1:452466/1‑141 (MQ=18)                                                       ccggCGTCTACGGTATTTCATTCCCCCCAACACCCGGCCATTACCCC‑AC‑CGCGATG‑GCAAAAAAG‑GCTT‑G‑GGGGGGGTGGGTCGACCATCTGTTTTC‑CCCC‑CGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAg                                                                                              <  1:337003/147‑1 (MQ=255)                                                               tcccgtttttCTTTCCCCCCAACACCCGCG‑ATTTTTCC‑AA‑CCGCAAGGTAAAAAATG‑GCTTCC‑GGGGGGG‑CTGCCTCCCAGTTTTTTTC‑CCGC‑GGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                     <  1:411282/143‑1 (MQ=255)                                                               tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAG‑GAAAAACTTTaaaa                                                                                      >  1:1233587/1‑148 (MQ=25)                                                               ccctgtgtTTCATTCCCCCAAACACCCG‑GCTTTTTACC‑AACCCCCGTGTGAAAAAATC‑CTTC‑C‑GGGGGGG‑CGGCCTCCCGGTTTTTTTC‑CCGC‑CAA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                     <  2:980553/147‑1 (MQ=255)                                                               cccggtatttCATCCCCGCAAACACCCGCGTTTTTCCCC‑AC‑CCCCTGGAGAAAAATTC‑CCTC‑C‑GCGGCGG‑TTGTCGCCCGGTTTTTTTC‑CCGC‑CAA‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGcc                                                                                     <  2:602424/143‑1 (MQ=255)                                                                aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAa‑aaaaaacattttacc                                                                                     >  1:1366820/1‑138 (MQ=25)                                                                      aTTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTAAGAG‑AAACAACTTTaacaccccgcgt                                                                              >  2:1404033/1‑141 (MQ=25)                                                                       ttcttcccccacaacaaCCGC‑TATTTAACCAAA‑CCCCGTGGGCAAAACTG‑CGTG‑G‑GGGGGGG‑CTGTCGGCAGGTTTTTTTC‑CCGC‑CAA‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgt                                                                             <  2:863313/140‑1 (MQ=255)                                                                        tctttccccccaccaCCCGCCCATTTACCC‑AC‑CCCCTTGTGCAAAAATC‑CTTC‑C‑GCGGCGG‑CTGTCCCCCAGCTTTTTTC‑CCGC‑CCC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgtt                                                                            <  1:763080/142‑1 (MQ=255)                                                                         ttttccccccAACCCCCGCGCTTTTCCCC‑AA‑CCCCATGTAAAAAAATG‑CCTT‑C‑GGGGGGG‑CTGTCTGCCAGTTTTTTTC‑CCCC‑CGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                           <  1:347240/143‑1 (MQ=255)                                                                         ctttccccccaccaCCCGCGCTTTTACCA‑AC‑CGCCATGTGAAAAAATC‑CCTC‑C‑GCGGCGG‑CTGTCTGCGAGTTATTTTC‑CCGC‑CGC‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                           <  2:191470/143‑1 (MQ=255)                                                                         ctttccccccAACACCCGCGCTTTTCCCC‑AA‑CGCCATGTAAAAAAATG‑CCTT‑C‑GCGGGGG‑CTTCCTCCCAGTTTTTTTC‑CCGC‑CGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                           <  1:199571/143‑1 (MQ=255)                                                                         catcccccccAACACCCGCGCTTTTACCC‑AC‑CGCCAGGTAAAAAAATC‑CCTG‑C‑GGGGCGG‑CTGTCTCCCAGTTTTTTTC‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCAtgttg                                                                           <  2:78311/143‑1 (MQ=255)                                                                              ccaccAACACCCGCGCTTTCACCC‑AA‑CCCCATGTAAAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTC‑CCGC‑CAC‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                      <  2:1529981/148‑1 (MQ=255)                                                                               ccccaccaCCCGCGCTTTCCCCC‑AC‑CCCCTTGTGAAAAAATG‑CCTT‑C‑GGGGCGG‑CTGTCTGCCAGCTTTTTTC‑CCGC‑CAC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                     <  2:499741/149‑1 (MQ=255)                                                                                               tttACCC‑AC‑CGCTATGAAAAAAATGG‑CTTT‑C‑GGGGCGG‑CTGTCGGCCTGCTTTTTTC‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                     <  1:176602/151‑1 (MQ=255)                                                                                               tttACCC‑AA‑CCCCTTGT‑AAAAAAAG‑CCTTCC‑GGGGCGC‑CTGTCTGCCAGTTTTTTTC‑CCGC‑CAC‑AACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                     <  2:597374/151‑1 (MQ=255)                                                                                                tttacc‑AA‑CCCCATTGGAAAAAAAG‑CCTTGC‑GCGGCGG‑TTGGATGCCAG‑TTTTTTC‑CCGCGCCA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                     <  2:1336393/147‑1 (MQ=255)                                                                                                ttATCC‑AA‑CCCCCTGGGAAAAAAGGCCCTT‑C‑GCGGGGC‑CTGTCGCCGAGTTTTTTTC‑CCCC‑CGC‑AAAAAAACCGCCGTTCTCTCCAGTAATGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                     <  1:1366406/151‑1 (MQ=255)                                                                                                   cccccaacccatggtgCAAAACTC‑CTTC‑C‑GCGGCGC‑CTTTCTGCTAGTTTTTTTC‑CCGC‑CAA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                   <  2:727415/138‑1 (MQ=255)                                                                                                     cccaccccctgggaaaaaaaGC‑CCTC‑C‑GGGGCGG‑TTTTC‑GCCCGCTTTTTTTCCCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                 <  1:706090/136‑1 (MQ=255)                                                                                                              ttgtcaaaaaaCC‑CCTC‑C‑GCGGCGG‑CTTTCTGCCATTTTTTTTC‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                        <  2:1206674/145‑1 (MQ=255)                                                                                                               tgtcaaaaagTG‑CTTC‑C‑GGGGGGG‑TTTTCGGCCAGTTTTTTTC‑CCGC‑CCA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                       <  1:1598891/146‑1 (MQ=255)                                                                                                               tgcgaaaaagTC‑CCTC‑C‑GCGGCGG‑CTGTCTGCCAGCTTTTTTC‑CCGC‑CGA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                       <  2:166928/146‑1 (MQ=255)                                                                                                               tgaacAAAAGTC‑CCCT‑C‑GCGGCGG‑CTGTCTGCCAGTTTTTTTC‑CCGC‑CGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                       <  2:109731/147‑1 (MQ=255)                                                                                                                gtgCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAACTTTAACCCCG‑CCTGTTTTGTGGGCGACGCCGCCGGTTAAAAAGTGGGAAGGGag                                      >  2:1399586/1‑150 (MQ=255)                                                                                                                    aaaaaTC‑CCTT‑C‑GGGGGGG‑TTTTCGCCCAGTTTTTTTC‑CCGC‑CCA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                  <  1:677849/151‑1 (MQ=255)                                                                                                                                           ttgtctccaatttATTTTC‑CCGC‑CGA‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  1:357415/139‑1 (MQ=255)                                                                                                                                           ttgcctgcagtttATTTTC‑CCCC‑CAA‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  2:441985/139‑1 (MQ=255)                                                                                                                                           ttGTCTCCAAGTTATTTTC‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  1:1270052/150‑1 (MQ=255)                                                                                                                                           ctttctgccgttttttttc‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  2:769005/137‑1 (MQ=255)                                                                                                                                           ctttctgccgttttttttc‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  1:340963/137‑1 (MQ=255)                                                                                                                                           ctgcctgcGAGTTTTTTTC‑CCGC‑CGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  2:1403422/151‑1 (MQ=255)                                                                                                                                           cTGTCTGCAGGTTTTTTTC‑CCGC‑CGA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  2:1105498/151‑1 (MQ=255)                                                                                                                                           cTGTCTCCCATTTTTTTTC‑CCCC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCGAGAAc               <  1:1053335/151‑1 (MQ=255)                                                                                                                                           cTGTCGGCCAGTTTTTTTC‑CCGC‑CGA‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc               <  1:844665/151‑1 (MQ=255)                                                                                                                                                 cccgttttttttc‑CCGC‑CGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc          <  1:1097336/142‑1 (MQ=255)                                                                                                                                                      tttttttc‑CCCC‑CCC‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg    <  1:1379978/148‑1 (MQ=255)                                                                                                                                                      cTATTTTC‑CCGC‑CGA‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg    <  2:1077036/151‑1 (MQ=255)                                                                                                                                                        tttttc‑CCCC‑CCC‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg  <  1:403014/150‑1 (MQ=255)                                                                                                                                                         ttttc‑CCGC‑CCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg  <  2:676442/151‑1 (MQ=255)                                                                                                                                                             |                                                                                                                                                      CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGT‑CAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCC‑AA‑CGCCATGTGCAAAAATG‑CCTT‑C‑GCGGCGG‑CTGTCTGCCAGCTATTTTT‑CCGC‑CGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAG‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAG  >  NC_000913/357731‑358025

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