Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,880	0	T	C	54.1%	47.4 / 91.5	61	K271E (AAA→GAA)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (26/2);  new base C (0/33);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 3.10e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-13

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357730‑358026                                                                                                                                                          |                                                                                                                                                      aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                       <  1:265689/151‑1 (MQ=255) aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                       <  1:346958/151‑1 (MQ=255)  cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑c                                                                                                                                                     >  1:911970/1‑151 (MQ=255)    cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑ccg                                                                                                                                                   >  1:750579/1‑151 (MQ=255)    cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑ccg                                                                                                                                                   >  2:903222/1‑151 (MQ=255)            tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                            >  1:648115/1‑150 (MQ=255)                    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                   <  2:112854/151‑1 (MQ=255)                    gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgcc                                                                                                                                   >  2:1112945/1‑151 (MQ=255)                     tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgccg                                                                                                                                  >  2:332735/1‑151 (MQ=255)                         gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtct                                                                                                                              >  2:329810/1‑151 (MQ=255)                             cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                          >  2:1059497/1‑151 (MQ=255)                               ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                          >  2:194929/1‑149 (MQ=255)                                  aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                       >  2:1234288/1‑149 (MQ=255)                                  aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                       >  2:6552/1‑149 (MQ=255)                                  aaTCCCAGCCCCGGACCTGTGTCGGCGGTTCCCGTATTTCATTCCCGCCAC‑CAACCGCGCTTTACCCCACCGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGCAAAAAAAccgccg                                                                                                                                  <  1:833313/138‑1 (MQ=14)                                   aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    >  1:933705/1‑151 (MQ=255)                                   aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                     >  2:214672/1‑150 (MQ=255)                                    ttCCAGCCGCAGACCTGTGTCGGCGGCTCCCTTGTTTCTTTCCCCCCAA‑CACCCGCCCTTTTACCCACCCCCATGTGC‑AAAATTC‑CCTCCGCGGCGGCTGTCTCCCGTTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagta                                                                                                                      <  2:1194886/148‑1 (MQ=255)                                    ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                   >  1:696447/1‑151 (MQ=255)                                    tcccggCCGCGGCCTTGTGCCAGCGGCTCCCGTGTTTCTTTCCCGCCAA‑CACCCGCGCTTTTACCCACCGCCATGTCA‑AAAAATG‑CCTCCGCGGCGGTTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    <  2:1069745/145‑1 (MQ=255)                                     tCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                                  >  2:916307/1‑151 (MQ=255)                                           cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                            >  1:961136/1‑150 (MQ=255)                                              agcctgtggccggCGGCTACCGTTTTTCTTTCCCCCCAC‑CACCCGCCCATTTA‑CCAACCCCCAGGGCAAAAAATC‑CCTCCGGGGGGCTTGTCGCCCAGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  1:986350/137‑1 (MQ=255)                                              agcccttttccAGCGCCTCCGGTGTTTCTTTCCCGCCAC‑CACCCGCGCTTTTACCCACCCCCATGGGA‑AAAAATG‑CCTTGGGGGGGGTTGTCTGCCGTTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                           <  2:193202/139‑1 (MQ=255)                                              agcccTGTGTCGGCGGCTACCGTGTTTCATTCCCCCCAA‑CAACCGCGCATTTACCCAACGCCATGTGA‑AAAAATG‑CCTCCGCGGCGGCTGTCTGCCAGTTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                         <  2:402876/148‑1 (MQ=255)                                              aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                             >  1:570104/1‑147 (MQ=255)                                                 ccTGTGTAGGCGGCTACCGTGATTCATTCCCCCCAA‑CACCCGCGCTTTTACCCACCCCCATGTGC‑AAAAATG‑CTTTCGCGGCGGCTGTCTGCCAGTTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAga‑g‑a                                                                                                    <  2:175879/151‑1 (MQ=21)                                                      gtCGGCGGCTCCCGTGTTTCTTTCCCCCCAA‑CACCCGCGTATTTACCCACCGCCATGGGA‑AAAAATC‑CTTGCGGGGCGGCTGTCGGCCGGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAt                                                                                                 <  2:1160199/149‑1 (MQ=255)                                                      gccgcgggcacccGTAATTCATTCCCCCCAA‑CACCCGCCATTTTACCCAACCCCTTGTAA‑AAAAATC‑CTTCCGCGGGGCCTGCCTCCAAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑aa                                                                                                   <  1:567842/136‑1 (MQ=255)                                                            gtctcccgtttttCTTTCCCCCCAA‑AACCCGCCCTTTTCCCAAACCCCTTGTGA‑AAAATTC‑CCTGCGGGGGGGCTGTCT‑CCAGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAg                                                                                               <  1:589440/134‑1 (MQ=255)                                                            ggCTACGGTGATTCATTCCCGCCAA‑CACCCGCTCATTTACCCAACCCCATGTGA‑AAAAATG‑CCTTCGGGCCGGTTGTCTCCCAGTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCttt                                                                                           <  2:761641/149‑1 (MQ=255)                                                            ggCTACCGTGATTCTTTCCCCCCAA‑CACCCGCGCTTTTATCCAACGCCATGTAC‑AAAAATG‑CCTGCGGGGCGGCTGTCTGCCAGCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTaa                                                                                         <  2:1205044/151‑1 (MQ=255)                                                               tACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAATaaaatt                                                                                           >  2:481212/1‑141 (MQ=255)                                                                aCCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑gaaaaaaatttaaccc                                                                                     >  2:195228/1‑138 (MQ=25)                                                                  ccgtttttctttCCCCCAACCACCCGCCCTTTTACCCAACGCCATG‑GA‑AAAAATGCCCTTCGGGGGGGTTTTCTGCCGGTTATTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAga‑g‑a                                                                                                    <  1:181566/126‑1 (MQ=255)                                                                     gATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGG‑AAAAAAATTTAAGcc                                                                                      >  1:814665/1‑146 (MQ=25)                                                                       ttatttcccGCAAA‑AAACCGCGAATTTACCCACCGCCATGTAC‑AAAAATG‑CTTCCGGGGCGGTTTTTTGCCAGCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGccc                                                                                     <  2:133435/137‑1 (MQ=255)                                                                       ttCTTTCCCGCCAA‑CACCCGCCTTTTTACCCACCCCCTTTTAA‑AAAAATG‑CCTCCGGGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCAtg                                                                               <  2:872545/150‑1 (MQ=255)                                                                       tactttccccccAA‑CACCCCCCCTTTTCCCCACCGCCATTTAA‑AAAAGTC‑CCTCCGCGGGGGTTGCCTGCCATTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTaa                                                                                         <  1:158964/130‑1 (MQ=255)                                                                              cccccAA‑AACCCGCGCTTTTACCCACCGCATTGTAA‑AAAAATG‑CCTCCGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                        <  2:1006117/147‑1 (MQ=255)                                                                                  caa‑caaCCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑G‑GAAAAAATTTTAACCCGTCACGGTTTGTgggg                                                                     >  1:533778/1‑147 (MQ=21)                                                                                      cacccGCCCTTTTACCCACCCCCATGTAA‑AAAATTG‑CCTCCGGGGCGCTGGTCTCCCATTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                  <  1:573724/146‑1 (MQ=255)                                                                                         ccgcgtttttccccAACCCCTTGTGA‑AAAAATC‑CCTCCGGGGCGGCTTTCTGCCGGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                             <  1:234971/139‑1 (MQ=255)                                                                                                tttaccaaccGCCGTGTGA‑AAAAATC‑CCTGCGGGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                         <  1:632244/139‑1 (MQ=255)                                                                                                tttACCCAACGCCTTGTGA‑AAAAATG‑CTTTCGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                        <  2:809712/149‑1 (MQ=255)                                                                                                 ttATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAAAACTTTAAATCCGCCCTGTTTTGTGGCGATCg                                                                 >  1:808891/1‑139 (MQ=37)                                                                                                       acccccttgtgA‑AAAAATC‑CCTCCGGGGCGGTTGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                <  2:481975/143‑1 (MQ=255)                                                                                                       acccccATGGGA‑AAAAAGG‑CCTCCGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                <  2:163306/146‑1 (MQ=255)                                                                                                       accGCAATGTGA‑AAAAAGG‑CCTTCGCGGCGGCTGTCTGCCGTTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                   <  1:1054077/145‑1 (MQ=255)                                                                                                       aaCCCCATGTAC‑AAAAATG‑CCTTCGCGGCGGTTGTCGGCAAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                               <  1:239940/151‑1 (MQ=255)                                                                                                         ccctaTGTGA‑AAAAATG‑CCTCCGCGGGGCCTGTCTCCAAGTTTTTTTC‑CCCCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                               <  1:884427/145‑1 (MQ=255)                                                                                                           ccATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑GAAAAAAAGTTTAAACCC‑TCCTGTTTTTTGGGGATCGCCGCTGGGGACAAAGTgg                                               >  1:194294/1‑147 (MQ=19)                                                                                                           ccATGTGA‑AAAAATG‑CTTTCGCGGCGGCTGTCTGCCGGCTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                               <  2:278012/147‑1 (MQ=255)                                                                                                           ccATGTAA‑AAAAATG‑CCTCGGGGGCGGCTTTCTGCCAGTTATTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                   <  2:1037713/143‑1 (MQ=255)                                                                                                            cATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑G‑AAAAAAATTTAAACCCCCCCTGTTTTGTGGGGATTGGCGCTGGGGACAAAGTGGgaa                                            >  2:967400/1‑149 (MQ=21)                                                                                                              tgtgA‑AAAATTC‑CCTCCGCGGCGGCTTTCTCCCATTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                         <  2:199620/150‑1 (MQ=255)                                                                                                                    aaaatTG‑CCTGCGCGGCGGTCGTCGTCCAGTTTTTTTC‑CCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                      <  2:350674/148‑1 (MQ=255)                                                                                                                      aaaTG‑CCTTCGCGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                 <  2:922042/151‑1 (MQ=255)                                                                                                                      aaaTG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGG‑G‑AAAAAAATTTAAAACCCCCCTGTTTTGTGGGGAACGCTGCTGGGAACAATGTGGAAAGGGAAgggggg                                 >  2:629566/1‑146 (MQ=255)                                                                                                                            cctcgggggcggcTGTCGCCCTTTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                            <  1:55694/144‑1 (MQ=255)                                                                                                                            ccTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGG‑G‑AAAAAAATTTTAACCCCCCCTGTTTTGTTGGGATCGCTGCTGGGGAAAAAGTGGGAAGGGGAggggggggg                              >  2:956126/1‑141 (MQ=255)                                                                                                                                               ccAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  2:957230/147‑1 (MQ=255)                                                                                                                                                    tATTTTC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg     <  1:812165/150‑1 (MQ=255)                                                                                                                                                       tttc‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  2:183117/150‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                      ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357730‑358026

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