Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	44.4%	59.7 / 32.5	63	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (23/12);  new base A (0/28);  total (23/40)
Fisher's exact test for biased strand distribution p-value = 1.65e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357741‑358033                                                                                                                                                           |                                                                                                                                                 tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGc‑aac                                                                                                                                                  >  1:1133585/1‑150 (MQ=255)  gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGc‑aaca                                                                                                                                                 >  2:595562/1‑150 (MQ=255)  gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGc‑aaca                                                                                                                                                 <  1:123213/150‑1 (MQ=255)     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAAc                                                                                                                                             >  1:96296/1‑151 (MQ=255)     aaGCGGAGCGCATCAGGCATTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATTTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACaaa                                                                                                                                               <  2:562095/149‑1 (MQ=255)      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAAcc                                                                                                                                            >  2:1250324/1‑151 (MQ=255)      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAAcc                                                                                                                                            <  2:1064981/151‑1 (MQ=255)          tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAAccgccg                                                                                                                                        >  2:87014/1‑151 (MQ=255)                   aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTcc                                                                                                                                >  2:1312274/1‑150 (MQ=255)                      cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAGAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCagt                                                                                                                             >  1:574490/1‑150 (MQ=255)                       aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCagt                                                                                                                             >  2:1395962/1‑149 (MQ=255)                           ccggCCGCAGACCTGTGTCGGCGCCTCCGGTGTTTCTTCCCCGCCAACAC‑CCGCCCTTTTCCCCACCCCCATGTGAAAAAATG‑CCTC‑CGCGGCGCCTTTCGGCCGTTTTTTTTC‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCagtagt                                                                                                                          <  1:601368/145‑1 (MQ=255)                           ccccccGCGGCCCTGTGTCGGCGGCTCCCGTGATTCTTTCCCCCAAACAC‑CCCCGCTTTTACCCACCCCAATGTAAAAAAATG‑CTTC‑CGGGGCGGCTTTCGCCCAGTTATTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                         <  2:655384/145‑1 (MQ=255)                             aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                     >  2:1020942/1‑151 (MQ=255)                              gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                    >  1:496131/1‑151 (MQ=255)                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                     >  2:802343/1‑148 (MQ=255)                                    gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                   >  2:988112/1‑146 (MQ=255)                                     aCCTGTGTCAGCGGCTACCGTGATTCATTTCTGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                                  >  1:1522006/1‑145 (MQ=255)                                       ttgtgTCGGGGACTCCCTTAATTCTTTCCCCCAAACCA‑CCGCCCATTTCTCCAACGCCAGGTGCAAAATTG‑CTTC‑CGGGGCGGTTGCCTGCCACCTTTTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                   <  2:613834/142‑1 (MQ=255)                                          tgtCGGCGGCTCCCGTGTTTCATTCCCGCCACAAC‑CCGGCCTTTTACCCAACCCGTTATGAAAAATTG‑CCTC‑CGCGGGGGTTCTCAGCCA‑TTATTTTTCCCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                    <  2:1311783/139‑1 (MQ=255)                                           gtaggcggccccggTGATTCTTTCCCGCAAACAC‑CCGCGCTTTTACCCACCGCCATGTGAAAAAATG‑CTTT‑CGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAg                                                                                                      <  2:1319759/138‑1 (MQ=255)                                                cGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑Gaaaaaaattt                                                                                                  >  2:1395417/1‑143 (MQ=255)                                                 gtctcccGTGTTTCATTCCCCCCAACAC‑CCGCGCTTTTTTCAACCCCCATGTCAAAAAATC‑CCTT‑CGGGGGGGCTTTCGCCCAGCTTTTTTC‑CCGCCGA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTaa                                                                                                <  2:178878/146‑1 (MQ=255)                                                 ggCTCCCGTGTTTCATTCCCGCCAACAC‑CCGCGCTTTTCCCCACCCCCATGTGAAAAAATC‑CCTT‑CGCGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTaa                                                                                                <  2:361296/151‑1 (MQ=255)                                                   cccccGTGATTC‑TTTCCCCCAACAACCCGCTCTTTCACCCACCGCCTTGTAAAAAAGTG‑CCGC‑CGGGGCGGCGGCCGGCCATTTTTTTTC‑CCGCCAA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAg                                                                                               <  1:197636/147‑1 (MQ=255)                                                      ccGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAA‑GAAAAAACttttaaac                                                                                             >  2:687251/1‑144 (MQ=25)                                                            ttCATTCCCCCCAACAC‑CCGCGCTTTTACCCACCCCCTTGTGAAAAAATC‑CCTC‑CGGGGGGGTTGTCTGCCACTTTTTTTC‑CCCCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtg                                                                                      <  2:1483911/150‑1 (MQ=255)                                                             tCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑Gaaaaaaat                                                                                                    >  1:332741/1‑130 (MQ=255)                                                                   ccGCCAACAC‑CCGCGATTTTACCCACCGCCTTGTAAAAAAATC‑CCTG‑CGGGGCGGCTTTCTGCCAGTTTTTTTC‑CCCCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTtgtg                                                                                 <  2:1432990/148‑1 (MQ=255)                                                                       caacaa‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAA‑GAAAAAAATTTTAACCCCTCCTGTTGTgggggga                                                                           >  1:1408250/1‑147 (MQ=25)                                                                             cccGCTCATTTACCAAACCCCATGTAAAAAAATG‑CCTT‑CGCGGCGTCTGTCTGCCATCTTTTTTC‑CCCCCAC‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                     <  1:74642/150‑1 (MQ=255)                                                                              ccGCGCTTTTACCCACCCCAATGGAAAAAAATG‑CCTC‑CGGGGGGGCTGTTGGCCAGTTTTTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                     <  2:815592/150‑1 (MQ=255)                                                                              ccGCCCTTTTACCCAACGCCATGTAAAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                     <  2:790094/150‑1 (MQ=255)                                                                               cgcctaTTTATCCAACGCCATGTACAAAAATG‑CTTT‑CGGGGCGGCTGTCGCCCAGTTATTTTG‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                     <  2:1340595/144‑1 (MQ=255)                                                                                 cgtaTTTACCCAACCCCTTGTG‑AAAAATGCCCTT‑CGGGGGGGTTGTCTGCCGGTTTTTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                  <  2:622650/147‑1 (MQ=255)                                                                                 cgcATTTCCCCACCCCTATGTGAAAAAAGG‑CTTC‑CGGGGGGGCTTTCGGCAAGCTTTTTTC‑CCCCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                                 <  1:1090480/151‑1 (MQ=255)                                                                                    aTTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAA‑GAAAAAAATTTAAAACCGTCATCTTTTGTggggaaatctgcttg                                                                 >  2:655476/1‑135 (MQ=25)                                                                                      ttccccAACGCATTGTG‑AAAAAAG‑CCTTCCGGGGCGGCTTTCTGCCATTTTTTTTC‑CCGCCCC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                              <  2:728644/145‑1 (MQ=255)                                                                                          ccACCGCCATGTCAAAAAAGC‑CCTT‑CGCGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                             <  2:546709/146‑1 (MQ=255)                                                                                               gCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑GAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGATCGCTGCTgggtaaaa                                                           >  2:182263/1‑137 (MQ=25)                                                                                                ccATGTGCAAAAATG‑CCTC‑CGCGGCGGCTGTCTGCCAGCTATTTTC‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                                  <  2:532216/151‑1 (MQ=255)                                                                                                    gtaaaaaaaTC‑CCTC‑CGGGGGGCCTGTCTCCCAGCTTTTTTC‑CCCCCCA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                              <  1:134515/147‑1 (MQ=255)                                                                                                        aaaaaTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAAAAAATTTTAACCCGGCATGTTTTGTGGCTATCAGCGGTGGTAGAAAAATGGAAAAGGGAGTgg                                          >  2:45074/1‑150 (MQ=255)                                                                                                         aaaaTG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAATTTTAACCCCCCCTGTTTTTTGGGGAAAGCTGCTGGGGAAAATGTGGGAAGGGGAggggg                                         >  2:1369875/1‑147 (MQ=255)                                                                                                          aattG‑CTTT‑CGGGGCGGCTGTCTCCCAGTTTTTTTC‑CCGCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                        <  2:1522155/148‑1 (MQ=255)                                                                                                          aattG‑CCTT‑CGCGGCGGCTGTCTCCCATCTTTTTTC‑CCCCCGC‑AAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                                        <  2:881688/148‑1 (MQ=255)                                                                                                                        ggcggcTGTCTGCCAGTTTTTTTC‑CCGCCAA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                             <  1:350071/150‑1 (MQ=255)                                                                                                                        ggcggcTGTCTGCCAGTTTTTTTC‑CCCCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGaa                                    <  2:1104037/143‑1 (MQ=255)                                                                                                                        ggcggcTGTCTGCCAGTTATTTTC‑CCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTc                             <  2:764330/150‑1 (MQ=255)                                                                                                                               gTCTGCCAGCTTTTTTC‑CCCCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                          <  2:396007/146‑1 (MQ=255)                                                                                                                                  cggccgtttttttt‑CCCCCACAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg                     <  1:582496/140‑1 (MQ=255)                                                                                                                                    ccattttttttc‑CCCCCGC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                    <  1:1513881/140‑1 (MQ=255)                                                                                                                                    ccAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAATTTTAAACCCCCCTGTTGTGTTGGGATCGGCGCTGGGGAAAATATGGGAAGGGGAGGGGGGGGGAACAACTCCCCAACCccgcc                 >  2:127003/1‑149 (MQ=255)                                                                                                                                         ttttttc‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt             <  2:609478/147‑1 (MQ=255)                                                                                                                                         ttttttc‑CCGCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg            <  1:182241/148‑1 (MQ=255)                                                                                                                                         ttttttc‑CCCCCGA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                      <  2:627371/138‑1 (MQ=255)                                                                                                                                           ttttc‑CCGCCGA‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa           <  2:1023247/149‑1 (MQ=255)                                                                                                                                            tttt‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAAAAGCTTTAAGCCCCTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGGAAGGGAAGAGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGtt        >  1:1319630/1‑151 (MQ=255)                                                                                                                                              ttcccgcccC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt         <  2:1133702/146‑1 (MQ=255)                                                                                                                                               t‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCGA‑GAAAAAACATTAAAACCCCCCACTTTTGTGGGGATTGCCGCTGGGTACAAAGTGGGAAGGGGAGGGGGGGGGCACAGGTCCAAAACCGCGCTAATTTAGTTCg      >  1:791616/1‑150 (MQ=255)                                                                                                                                                cccgccgA‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg      <  1:1031988/149‑1 (MQ=255)                                                                                                                                                 ccgccaC‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg      <  2:212194/149‑1 (MQ=255)                                                                                                                                                 cccccga‑AAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGcc    <  1:1229248/148‑1 (MQ=255)                                                                                                                                                  cgccgc‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAAAATTAA‑AACGTCCCGCTGTGCGGCGATAATTGGTGGTAACAACAGGGTAAGGGGAGGGGGGGGGATCAGGTCCAGAACCGAGCAACTTTAGttt       >  1:1013265/1‑146 (MQ=255)                                                                                                                                                  ccccga‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt   <  2:574568/149‑1 (MQ=255)                                                                                                                                                    ccgc‑AAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt   <  2:635369/149‑1 (MQ=255)                                                                                                                                                         aacaaAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc  <  1:353221/146‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                 TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATG‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGC‑AACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC  >  NC_000913/357741‑358033

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