Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000964	1,143,719	0	A	G	32.4%	45.0 / 12.1	37	D48G (GAT→GGT)	sbcD	DNA repair exonuclease
Reads supporting (aligned to +/- strand):  ref base A (6/19);  new base G (12/0);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGGCACATTCTGACATTATAGCGAGATCCATAAGCTCCGGAATTTCAGGCGGAGCGGGTCTCGTTTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGTGCCTTCAGCAGGAGAGCTTTTGGCGGTTGGAGCGCTGGC  >  NC_000964/1143483‑1143968                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                gAGGGCACATTCTGGCATTATAGCGAGATCCATAAGCTCCGGAATTTCAGGCGGAGCTGGTCACTTTTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTATATGATGTGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATGTGAC‑AC‑Cg                                                                                                                                                                                                                                                <  2:70721/251‑1 (MQ=255)            cTGACATTATAGCGAGATCCATAAGCTCCGGAATTTCAGGCGGAGCGGGTCTCGTTTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGGGG‑ATTTTAC‑AC‑CGTaaaacccccc                                                                                                                                                                                                                                     >  1:78342/1‑244 (MQ=255)              gACATTATAGCGAGATCCATAAGCTCCGGAATTTCAGGCGGAGCGGGTCTCGATTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACCGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGGAATGGCGGGGGGTGCCATTTTAC‑‑C‑CGTAAA‑ACCGCCaac                                                                                                                                                                                                                                  >  1:88499/1‑249 (MQ=255)                      aGCGAGATCCATAAGCTCCGGAATTTCAGGCGGAGCGGGTCTCGTTTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGGGC‑ATTTTAC‑AA‑CGTAAA‑ACCGcccaccttagccg                                                                                                                                                                                                                          >  1:308707/1‑240 (MQ=255)                                 tAAGCTCCGGAATTTCAGGCGGAGCGGGTCTCGTTTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAAATGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGTGC‑ATTTGACGAC‑CG‑AAA‑ccccacagatttaacacaggaactggtt                                                                                                                                                                                                               >  1:268001/1‑226 (MQ=255)                                                catgcaggcggggTCTGGTTTTCCATATAAAAGAGTGGTGAATGCGTTAGCGACTTTGACATGCGGTTGAGTGTCATCTTGGCAAAACGCATGAAGGGAGATTGATGCTGAGTGAACAGGCGGATGTGTGGGATGAATGGAATACAATTGTAAGGGATGACTAAACGGAGGCAATTTTATTGGCGGGGGATGC‑ATTTGGG‑CC‑CGTATA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAATGCCTGTc                                                                                                                                                                                               <  2:100440/240‑1 (MQ=255)                                                  ggcggggcggGTCTCGTTTTCTATAAAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTgc                                                                                                                                                                                            <  2:209489/251‑1 (MQ=255)                                                     ggAGCGGGTCTCGTTTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCt                                                                                                                                                                                         <  1:151139/251‑1 (MQ=255)                                                                   tttttatataAAAGAGAGGTGAAAGCATTTGCGGATTTTAAATACGGCTGCCTGGCATATTGGAAAAACGCTTGAAGGAAGAAGCAGGCGGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTGGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGa                                                                                                                                                                            <  2:210797/246‑1 (MQ=255)                                                                      cTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGTGG‑ATTTGGA‑AC‑CGTAAA‑CCCCCCCACTTTATCCG‑GGCGGCTG‑TTTTTTTAAAGACCGTCTGCGCttttcgacggcgtacac                                                                                                                                                                         >  2:347176/1‑235 (MQ=255)                                                                             aaaGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGGGG‑ATTTGAC‑AC‑CGCAAA‑AACCCCACCCTTAGCCG‑AGGCGCAG‑GTTTTTTAAAAACTGTCTGCGCTTAGAGACAGAGGAAAGCGccccat                                                                                                                                                                 >  1:240457/1‑248 (MQ=255)                                                                                 agagGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTTAC‑AC‑CGTAAA‑CCCGCCCGCTTTATCCG‑AGCAGCTG‑TTTTATTAAAGCCTGTCTGCGCTTAGCGACCAAGGGAAGCGccccctcgtc                                                                                                                                                             >  1:139643/1‑244 (MQ=255)                                                                                    ggTGAAAGCATTTGCGGATTTTACATACGGCGAACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCGGAGGGAACAGGCCGATTTGCTGGATGAACTGAAAACGATGGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑CC‑CGTAAC‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCgtcgtc                                                                                                                                                          <  2:343979/251‑1 (MQ=255)                                                                                       aaaatcATTTGCGAATTTTAAATACGGCTACCTGGCATTTTGGAAAACCGCTTAAAGGAAGAAGCAGGCTAAGTAAACAGGCCGATGTGTTGGATGAATTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTTTAATGGCGGGGAATGC‑TTTTGAC‑CC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAACGCCCGATCGTCGTCAtt                                                                                                                                                       <  2:337626/246‑1 (MQ=255)                                                                                          aGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGAGC‑ATTTGA‑‑ACACGGAAA‑CCCCCCAGCCTTAACCC‑GGAAGCTG‑GTTTATGAAAACCTGTCTTCGCTTTGCGACCAAGGGAAAAGCCCGATCGTCGTCATTGcc                                                                                                                                                    >  1:278629/1‑251 (MQ=255)                                                                                               ttGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAC‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGaa                                                                                                                                                <  1:346534/250‑1 (MQ=255)                                                                                                     aTTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGTTC‑CTTTTGC‑AC‑CGTAAA‑CCCCGCAACTTTATCCGCGGCA‑CTT‑TTTTTTGGAAACATGTCTTCGCTTTGCGAAAGCGGAAAGAGCCCCAGCGCCGTCATTGTCGGAAATCACgc                                                                                                                                         >  1:306680/1‑250 (MQ=255)                                                                                                                             cTTGGCAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGGGC‑ATTTTAC‑AC‑CCTAAA‑ACCCCCACCCTTTGCCG‑AGCAGCTGGTTTTTTGAAAGC‑TGTTTGCGCTTTGCGCCAGAGGAAAGCGCCCCCTCGTCGTCATTGGCGGAAATTACCATAATCCTTACCCGTTGTCcgccgc                                                                                                                 >  1:167270/1‑251 (MQ=255)                                                                                                                                         gtggaaggaagAAGCAGGCTGAGTGAGCAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATCCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAAGCCTGACCGTCTGTCCGCCGCTTCACCGCTGac                                                                                                     <  2:99183/248‑1 (MQ=255)                                                                                                                                            gaaggaagAAGCAGTCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATGGTAATGGCGGGGGATGC‑ATTTGAC‑CC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCg                                                                                                          <  2:231935/243‑1 (MQ=255)                                                                                                                                            gaaggaagAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGGGG‑ATTTGAC‑AC‑CGTAAA‑CCCCCCAGCTTTATCCG‑AGCAGCTG‑TTTTTTGAAAGCCTGTCTGCGCTTTGCGGCAGAGGGAAACGCCCGGTCGTCGTCATTTCCGGAAATAACGATAATCCTGACCGTTTGTCCGCCGCTTCAcccggatcggg                                                                                                  >  1:231935/1‑242 (MQ=255)                                                                                                                                                              ggTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGGGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACCCATGAAAACGGCATTCAttt                                                                                <  1:87526/250‑1 (MQ=255)                                                                                                                                                                     aGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGAAGC‑ATTTGAC‑AC‑CCTAAAACCCCCCCG‑TTTAGCCG‑AGCAGCTG‑TTTTATGAAAACCTGTCTGCGCTTTGCGACAGAGGAAAGCGCCCGGTCGTCGTCATTGACGGAAATAAAGAAAAATCTTACCGGTTGTTCGCCCCTTCACCGCTGACACAAGAAAACAGcatttcattaattcgg                                                                         >  1:252144/1‑239 (MQ=255)                                                                                                                                                                                  ggATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGacaac                                                              <  1:297756/249‑1 (MQ=255)                                                                                                                                                                                  ggATGAACTGAAGACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑CC‑CGAAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGacaac                                                              <  2:98862/249‑1 (MQ=255)                                                                                                                                                                                               aCGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTGGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCatat                                               <  2:162377/251‑1 (MQ=255)                                                                                                                                                                                                cGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGTGC‑ATTTGAC‑AC‑CGTAAA‑ACCGCCCGCTTTAGCCG‑AGCAGCTG‑TTTTTTGAAAACCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCCATCGTCGTCATTGTCGGAAATCACGATAATCCTGACCGTTTGTCCTCCGCTTCACCGCCGACACATGAAAACGGCATTCATTTTATCGGGTATTCGACAACAGAGCCGCTTCATAtt                                              >  1:178976/1‑251 (MQ=255)                                                                                                                                                                                                    tGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGAGGC‑CTTTGGC‑AC‑CGTAAA‑CCCCCCCGCCTTTTCCG‑AGGAGCTG‑TTTTTTTAAAACCCGTCTTCGCTTAGAGACCGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGGAAATCCTGACCGTTTGTCCCCCGCCTTACCGCTGACACATGAAAAAGGCATTCATTTAATCGGTTATCCCCCAAAAGAGCCGAGTCATAtttaa                                           >  1:48161/1‑247 (MQ=255)                                                                                                                                                                                                      tAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGGTGC‑ATTTGGC‑AC‑CGTAAA‑ACCGCCCGCTTTATCCG‑AGCAGCTG‑TTTTTTTAAAGGCTGTCTGCGCTTTGCGGCAGAGGAAAGCGCCCGAACGTCGTCCTTGCCGGAAATCACGATAATCCTGACCGGTTGTCCGCCGCCTCACCGCTGACACATGACAACGGGATTCATTTTATCGGGTATCCGACAACAGAGGCGACTCCTATTGAAGTg                                        >  1:75936/1‑251 (MQ=255)                                                                                                                                                                                                      gaaaGGAGGAGCAGATCGATGCCATTGTAAGGGCGGTGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGt                                         <  2:31915/249‑1 (MQ=255)                                                                                                                                                                                                       aaaGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGTg                                        <  1:36995/250‑1 (MQ=255)                                                                                                                                                                                                       aaaGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑ACCGCCCGCTTTTAACG‑AGCAGCTG‑TTTTTTGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGAACGTCGTCATTTCCGGAAAACACGATAATCCTGACCGTTTTTCCCCCGCTTCACCGCTGACACATGAAAAAGGCATTCATTTTATCGGGTATTCCACAACAAAGCCGATTTATATTGAAGTgg                                       >  1:70699/1‑250 (MQ=255)                                                                                                                                                                                                                       gATGCCATTGTAATGGCGGGGGATGC‑ATTTGTC‑AC‑CGTAAA‑CACGCCCGCTTTAGCCG‑AGCAGCTG‑TTTTTTGTAAACCTGTCTGCGCTTAGCGACCGAGGAAAGCGCCCGAACGGCGGCATTGTCGGGAATAAAGAAAAACCTGACCGTTTTTCCCCCGCTTCACCCCTGACACATTAAAAACGCATTAATTTAATCGGGTTTTCGACAAAAGAGCCGCTTCATATTTCAGTGCCTTTAGCAGgaggggt                       >  1:346795/1‑247 (MQ=255)                                                                                                                                                                                                                          gCCATTGTAAGGGCGGGGTATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGTGCCTTCAGCAGGAGAGCtttt                    <  2:162288/251‑1 (MQ=255)                                                                                                                                                                                                                                gTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGTGCCTTCAGCAGGAGAGCTTTTGGCGGt              <  1:347176/251‑1 (MQ=255)                                                                                                                                                                                                                                        gggggATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGTGCCTTCAGCAGGAGAGCTTTTGGCGGTTGGAgcgc      >  1:106826/1‑251 (MQ=255)                                                                                                                                                                                                                                            gATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGTGCCTTCAGCAGGAGAGCTTTTGGCGGTTGGAGCGCTGGc  <  2:78342/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                GAGGGCACATTCTGACATTATAGCGAGATCCATAAGCTCCGGAATTTCAGGCGGAGCGGGTCTCGTTTTCTATATAAAAGAGAGGTGAAAGCATTTGCGGATTTTACATACGGCTGACTGGCATCTTGGAAAAACGCTTGAAGGAAGAAGCAGGCTGAGTGAACAGGCCGATGTGCTGGATGAACTGAATACGATTGTAAAGGATGAGCAGATCGATGCCATTGTAATGGCGGGGGATGC‑ATTTGAC‑AC‑CGTAAA‑CCCGCCAGCTTTAGCCG‑AGCAGCTG‑TTTTATGAAAGCCTGTCTGCGCTTAGCGACAGAGGAAAGCGCCCGATCGTCGTCATTGCCGGAAATCACGATAATCCTGACCGTTTGTCCGCCGCTTCACCGCTGACACATGAAAACGGCATTCATTTAATCGGTTATCCGACAACAGAGCCGATTCATATTGAAGTGCCTTCAGCAGGAGAGCTTTTGGCGGTTGGAGCGCTGGC  >  NC_000964/1143483‑1143968

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