Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	51.7%	38.8 / 104.1	60	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (23/6);  new base A (0/31);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 2.03e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-08

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAT  >  NC_000913/357741‑358036                                                                                                                                                          |                                                                                                                                                     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)                                                                                                                                                    ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa    <  2:96224/151‑1 (MQ=255)                                                                                                                                                     cgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAGTTTTAAACCC‑TCCTGTTTTGTGGGGAACGCTGCTGGGGACAATGTGGGAAGGGGAGGGGGGCGGATCAGCTCCCCAACCGCGCCAATTTAGGTTGCCTc     >  2:344253/1‑149 (MQ=255)                                                                                                                                                       caacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAAAAACTTTAAACCCCCCATCTTTTGTGGCGAACGCTTCTGGGTAAAATTTGGGAAGGGGAGGGGGGGGCATCAACTCCCCAACCCCCCCAATTAAGTTTGCCTc     >  2:954600/1‑147 (MQ=255)                                                                                                                                                       caaaaaATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAt  <  2:996180/150‑1 (MQ=255)                                                                                                                                                         acaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAt  <  1:934589/148‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                     TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTATCCAACGCCATGTGC‑AAAAATG‑CCTTCGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑G‑AAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAAT  >  NC_000913/357741‑358036

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