Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	357,834	G→C	17.7%	A286G (GCG→GGG)	cynR ←	transcriptional activator of cyn operon, autorepressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,834	0	G	C	17.7%	69.5 / 3.7	44	A286G (GCG→GGG)	cynR	transcriptional activator of cyn operon, autorepressor
Reads supporting (aligned to +/- strand):  ref base G (25/10);  new base C (2/6);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 3.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTT  >  NC_000913/357599‑358062                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           attcattcCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTCTCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCACCCGCAGACCTGTGTCCGCGGCGACCGTGGTTCCTTTCCGCCCACCACCGCCCATTTTTCCAC‑CGcc                                                                                                                                                                                                                                          >  2:1328848/1‑238 (MQ=21)                           tcctccAACTCCATGTGCAGAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATCCCCGATACGGCGTTCAACCCTTTTCCCACCGACGTAGAACACCAACCTCGTAGGCGTGAGAAGCGGAGTCGC‑CCAGGCGGTTTGAGGCGCAGAACTGCGGCGGCGGCTGCCGGGATTTATTTCCTCCAACAACCCCGAATTTATACTA‑CTCACT‑GTGCAAAAATCCCCT‑CCCCGCGG‑Ct                                                                                                                                                                                                            >  2:1102823/1‑239 (MQ=255)                           tcctccAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACAACCGCGCAATTATCCAA‑CGCCAT‑GTTCAAAAATGCCTT‑Cgcggcg                                                                                                                                                                                                                >  1:729728/1‑236 (MQ=255)                                       aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGCACGCCTTATCCGGCCTACGTAGAGCACTGAAATCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCCGCCGCAGACCTGTGTCAGCGGCTAACGTGATTCATTTCCGCCAAAAACCCCGCATTTATCCCA‑CGCCAT‑TTGCAAAAATGCCCT‑CGCGGCGG‑CTGGCTCCCCGCTa                                                                                                                                                                                                >  2:952476/1‑239 (MQ=18)                                       aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCCTTTATCCAA‑CCCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTa                                                                                                                                                                                                >  1:84942/1‑239 (MQ=255)                                       aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGcca                                                                                                                                                                                                >  1:1083380/1‑237 (MQ=255)                                       aTGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTa                                                                                                                                                                                                >  1:1305601/1‑239 (MQ=255)                                        tgtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCCTGATTCATTTCCGCCAACAACCCCCCATTTATCCAA‑CGCCA‑‑‑TGCAAAAATGCCTT‑CGCGGCGG‑CGGTCTGCCACCCAttt                                                                                                                                                                                             >  1:1204456/1‑239 (MQ=18)                                        tgtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGAGGCTACCGTGATTCATTTCCGCCAACAACCCCGCATTTATCCCA‑CGCCCT‑GCGCAAAAATGCCTT‑CGCGGCGG‑CTGGCTGCCAGCttg                                                                                                                                                                                               >  2:243771/1‑237 (MQ=255)                                         gtgCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAAAGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTGCCGCCCACAACCGCCCCTTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGACGG‑GCGTCTGCCAGCTAtt                                                                                                                                                                                              >  1:1059123/1‑239 (MQ=17)                                                    ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCCCCCATTTATCCAA‑CGCCAT‑GTGCAAAAACGCCTT‑CGCGGCGG‑GCGCCTGCCCGCTCTTTTTtcgccccc                                                                                                                                                                                   >  2:1240645/1‑236 (MQ=17)                                                    ccTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACGTAGAGCACTGGACCCGTAGGCCTGATAAACGGAG‑CGCATCAGGCAATTCCAGCCGCAGACCTTGGTCAGCGGCTACCGTGATTCATTTCCGCCAACCACCGCGCCTTTAGCCAA‑CGCCAT‑GTGCAACAATGCCTT‑CGCCGCGG‑CTGTCTGCCAGCCATGTTTccgcctc                                                                                                                                                                                    >  2:501507/1‑236 (MQ=17)                                                                          gCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCCGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CCCCCT‑TGGCAAAAATGCCTT‑CGCGGCGG‑CCGTCTGCCAGCCATTTTTTCCCCCCAACAAAAACGCCCGTCCCTccc                                                                                                                                                              >  2:696177/1‑237 (MQ=14)                                                                               aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGcc                                                                                                                                                                                                                                          <  2:1300732/160‑1 (MQ=40)                                                                               aGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGcc                                                                                                                                                                                                                                          >  1:1300732/1‑160 (MQ=40)                                                                                 tttATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                                                      >  1:175765/1‑239 (MQ=255)                                                                                 tttATGCCGGATGCGGCGTGAACCCCTTATCCGGGCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGAGGAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAAAGCCTT‑CGCGGCGG‑CTGTCTCGCAGCTATTTTTCCGCCGCAAAAAAACCCCCGGtctctctcgaagttgg                                                                                                                                                      >  2:1248628/1‑229 (MQ=17)                                                                                     tGCCGGCTGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAAGGCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAAACCTGTGTCAGCGGCTCCCGGGGTTCATTTTCGCCCTCAAACGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGACGCGG‑GTGTCTGCCACCTCTTTTTTCGCCGCAACCAAAACCCCGTTGTCTCCCGTCGTGCT‑gggg                                                                                                                                                 >  2:1264451/1‑239 (MQ=14)                                                                                      gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGGACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCCCGCATTTATCCAA‑CGCCAT‑GTGCCAAAATGCCTT‑CGCGGCGG‑CCGTCTTCCCGCTCTTTTTTCCCCGCCACAAAACCGCCGCTCTTTCCCGGAAGGGC‑GGGGc                                                                                                                                                >  2:1227977/1‑239 (MQ=14)                                                                                      gCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCCGCTATTTTTTCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGC‑GGGGc                                                                                                                                                >  2:8340/1‑239 (MQ=255)                                                                                          gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGC‑gggggctag                                                                                                                                            >  1:128078/1‑234 (MQ=255)                                                                                           aTGCGGCGTGAACCCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCCT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCCGTCTCTCCAGTAGTGGC‑gggcggcggg                                                                                                                                           >  1:743392/1‑232 (MQ=255)                                                                                               ggcgagcacctcTTACCCGGCCTCCGAAGGGCCCTGAACTCGTAGCCCTGAAAGGCGTGG‑CGCTCCAGCATACCCCCGCCCCGGCCCTGTGTCAGGGGCTACGGTGATTCTTCCCCCCCACCACCCCCGCTTTAACCCAC‑CCCCAT‑GTCAAAAAATCCCTT‑CGGGGCGG‑TTGTCGGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGaaa                                                                                                                                       <  2:87675/227‑1 (MQ=255)                                                                                                             aTCCGGCCTACGTAGACCTCTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCATCCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑cccc                                                                                                                                                                                                                                          <  2:1275339/130‑1 (MQ=14)                                                                                                                  gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCCGTAGTGGC‑GGGGCAACGAGAAAA‑AG‑TTTAATccccgctttg                                                                                                                  >  1:525869/1‑232 (MQ=18)                                                                                                                    cTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAC‑GTGGAAT‑aaattttttaccctcccgttg                                                                                                                >  1:1326757/1‑219 (MQ=18)                                                                                                                       cGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCaaaa                                                                                                                                                                                                                               >  1:1275339/1‑130 (MQ=255)                                                                                                                                 tGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTg                                                                                                                                                                                                       <  1:893224/142‑1 (MQ=255)                                                                                                                                 tGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTg                                                                                                                                                                                                       >  2:893224/1‑142 (MQ=255)                                                                                                                                 tGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAA‑AGGAAAA‑AGATTTAA‑CCCCCCATGGTGTTTGGCGATCcct                                                                                                   >  1:87052/1‑239 (MQ=18)                                                                                                                                    ccTCGTAGGCCTGATAAGCGTAC‑CGCATCAGGCAATCCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgc                                                                                                 <  2:627455/237‑1 (MQ=255)                                                                                                                                                       gTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCA‑GAGGAAA‑AACTTTTAAACACTCTTGTTGTGTTGGGAAGGCTGCTGGGAACCACGTgggtatgg                                                                             >  1:1192581/1‑233 (MQ=18)                                                                                                                                                             gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCA‑GAGAAATAAGATTTAA‑CCCCTCACTTTGTTTGGGAATAGCTGGTGGGTACAAAGTGGGAAGGGGAGTg                                                                        >  1:470435/1‑239 (MQ=21)                                                                                                                                                                                                           tccatccccGCCAACACCCGCACTTTACCCCAC‑CCCCAT‑TTGAAAAAATGCCTT‑GGGGGCGG‑TTGTCGGCCGTCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGa                          <  1:286968/232‑1 (MQ=255)                                                                                                                                                                                                            tctttcccgccacccacgcctaaTTTATCCAA‑CCCCAG‑GAAAAAAAAGCCCTC‑GGCGGCGT‑CTGCCTGCCACTTTTTTCCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                         <  2:390316/217‑1 (MQ=255)                                                                                                                                                                                                            tctttccccccAACACCCGCGCTTTTACCCAA‑CGCCAT‑GTCAAAAAATGCCTT‑CGCGGCGG‑TTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                         <  1:1298912/230‑1 (MQ=255)                                                                                                                                                                                                            tctttcccGCCAACCCCCGCGCATTT‑CCCAA‑CCCCATGGTCAAAAAATCCCTT‑CGCGCGGG‑CTGTCGGCCAGTTATTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                         <  2:1305601/233‑1 (MQ=255)                                                                                                                                                                                                            tctttcccGCAAACACCCGCGCATTTACCCAA‑CCCCAT‑GTAAAAAAATGCCTG‑CGGGGCGG‑CGGTCCGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGac                         <  1:8340/233‑1 (MQ=255)                                                                                                                                                                                                              ctttccccccaacaccccccATTTAACCAACCCCCAT‑GTGCAAAATGGCCTC‑CGCGGCGG‑CTGTCTGCAATTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGacc                        <  1:221079/220‑1 (MQ=255)                                                                                                                                                                                                               ttcccGCCAACACCCCCACTTTTACCCAC‑CGCCAT‑GTAAAAAAATGCTTC‑CGCGGCGG‑CTGCCTGCCATTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGaccac                      <  1:134787/236‑1 (MQ=255)                                                                                                                                                                                                                   cccaaacaccCGCGCTTTTACCCAC‑CGCCTG‑GTGCAAAAAGGCCTC‑CG‑GGCGGTTTTTCGCCCAGTTTTTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCtgt                  <  1:508315/235‑1 (MQ=255)                                                                                                                                                                                                                            ccGCGAATATACCCAA‑‑GCCAT‑GGTCAAAAACGCCTTCCGGGGCGG‑TTTTCTGCTGGTTTTTTCCCCGCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGc         <  2:1149998/239‑1 (MQ=255)                                                                                                                                                                                                                              cccattttccccAC‑CGCCAT‑GGAAAAAAATCCCTC‑GGCGGCGG‑TTGTCGGCCATTTTTTTTCCCCCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGccc       <  1:193107/229‑1 (MQ=255)                                                                                                                                                                                                                                   tttACCCAC‑CGCATT‑GTGAAAAACGCCCTC‑CGCGGCGG‑TTGTCTGCCATCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCtt  <  1:243771/239‑1 (MQ=255)                                                                                                                                                                                                                                   tttACCCAA‑CGCAAT‑GGAAAAAAAGGCCTC‑CGCGGCGG‑CTGTCGGCCATCTTTTTCCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCtt  <  1:105761/239‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                           ATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAG‑CGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAA‑CGCCAT‑GTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAA‑GAGAAAT‑AGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACCACCTGTGGATGTAGCCCCGCTT  >  NC_000913/357599‑358062

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