Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	71,569	T→C	I73I (ATT→ATC)	yabI →	DedA family protein YabI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	71,569	0	T	C	97.4%	110.2 / ‑6.1	38	I73I (ATT→ATC)	yabI	DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (11/26);  total (12/26)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCG‑ACTGGATTTCTTTCTGGCTGGGTTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCC  >  NC_000913/71337‑71808                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  ttGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTgggctg                                                                                                                                                                                                                                  <  1:37213/250‑4 (MQ=255)    aCAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGt                                                                                                                                                                                                                              <  1:626357/251‑1 (MQ=255)                   aaGCATTGCTGGAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTg                                                                                                                                                                                                              <  1:701454/251‑1 (MQ=255)                         tGCTGGAACACTTTATTACCCAGTCTACCGTGTATTCATTGATGGCGGTGGTTTTGGTGGCCTTTCTGGAGTCGCTGTCGCTGGTCGGTTTGATTCTACCCGGTACTGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCGGCGGCGAGTTAATTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTATCTGCTGGTCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCg                                                                                                                                                                                                         <  1:108024/250‑1 (MQ=255)                              gAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTc                                                                                                                                                                                                   <  1:603236/251‑1 (MQ=255)                                     ttattaCCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTTGCAGCGGCGAGGTAAGTTTTTGGCACGCCTGGCTGGCCGGGGTTGTCGGGTGTCTGCTGGGCG‑GCAGGGTCTCCTTCTTGCTGGGCTGGCGGTTT‑TAAAAGCCCTTGCATCGCTGGTCCTTTCtt                                                                                                                                                                                            >  1:340538/1‑250 (MQ=255)                                         tcaccaatcttccGTGTATTCATTGATGGCGGGGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGGCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTg‑aag                                                                                                                                                                                        <  2:2508/239‑1 (MQ=255)                                            ccAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGaaaa                                                                                                                                                                                    <  1:498181/251‑1 (MQ=255)                                                    gtgtggtttCATTGAGGGCGGGGGGTTTGGAGTCCTTTCTTGAGGTGCTGGTGGTGGTGGGTTTGATTCTGCCCGGGATGGTGCTGATGTTGGGGCCGGGTGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTTGCACGCCTGGCTTGCGGGGATTGGCGGTTGGTTGCTGGGCG‑CCTCGATCTCTTTCTGGCTGTGAAGGCGTTTT‑AAAAAGCCGGTGCATCGCTGGTCATTACTG‑AAAAAAAACTAAGc                                                                                                                                                                             <  2:24723/243‑1 (MQ=255)                                                        gtATTTATTGATGTCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCTGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGTAGGGATTGTCGGTTGTCTGCTGGGCG‑TCTGGATCTCCTTTTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCattac                                                                                                                                                                        <  1:82098/251‑1 (MQ=255)                                                          aTTCATTGATGGCGGTGGGGTTGGTGGCCTTTCGGGAGTCGTGGGCGCTGGTCGGTTTGATTCTGCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCTGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCattactc                                                                                                                                                                      <  2:4773/251‑2 (MQ=255)                                                              aTTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCTCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATa                                                                                                                                                                  <  1:375760/251‑1 (MQ=255)                                                              aTTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGGTTGGCAGCGGCGAGTTAAGTTTCTGGCAAGCCTGGCTGGCAGGGGTTTTCGGGTGTCTGCTGGGGG‑ACTGGGTCTCCTTTTGGCGGGGCTGGCGTTTTTAAAAAACCGTT‑CATCCCTGGGCCTTTTTG‑AAGAAAAACaaaaaattaatcgat                                                                                                                                                                   >  1:289168/1‑238 (MQ=255)                                                                ggATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTTGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AATAAAAACAAAGCATTACTCGATaaa                                                                                                                                                                <  1:465011/250‑1 (MQ=255)                                                                      cGGTGGTGTTGGTGGCCTTTCTGGAGGCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGaa                                                                                                                                                          <  1:692147/251‑1 (MQ=255)                                                                                 gtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCa                                                                                                                                               <  1:132109/251‑1 (MQ=255)                                                                                           tGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGGTTGTCGGTTGTCTGCTGGGCG‑GCTGGGTCTCCTTCTGGCTGGGCTGGCGTTTT‑TAAAAGCCGGTGCATCGCTGGTCATTTCTG‑AAGAAAAAAAAAGCATTACTCGATAAAACTGAACATGCCGTGCATCAACACAGc                                                                                                                                     >  1:427604/1‑251 (MQ=255)                                                                                           tGGAGTCGCTGGCGCTGGCCGGTTTGATTCTCCCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGTTTGGCAGCGGCGGGTTAAGTTTCTGGCACGCCTGGCTGGCGGGGATTGTCGGTTGTCTGCTGGGGG‑ACTGGATTTCCTTCTGGCGGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCTTCAaccacgc                                                                                                                                     >  1:686860/1‑246 (MQ=255)                                                                                             gAGTCGCTGGCGCTGGTCGGTTTGATTCTACACGGTACGGTGCCGATGGCGGGGATGGGGACGCTTATTGGCAGCGGCGAGGTAAGGTTTTGGCACGCCTGGCTGGCAGGGATTGTCGGGTGTCTGGTGGGCGGAC‑GGATCTACTTCTGGCTGGGCTGGAGTTTT‑AAAAAACCGTTGGATCGCTTGGCATTTCTGCAGGAAAAACAA‑GAATTTCTCGATAAAACTGCACAAGCGATGCATCaaaccaacat                                                                                                                                   >  1:629131/1‑243 (MQ=255)                                                                                              ggTCGCTGGCGCTGGTCGGTTTGATGCTACCCGGTACGGTGCTGATGGCGGGGCTGGGGGCGCTGATGGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATg                                                                                                                                  <  1:434101/250‑1 (MQ=255)                                                                                                 cgctggcgctggTCGGTTTGATTCTGCCCGGGACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAATCTGAACATGCGTTGCATCAACACAGCATGtt                                                                                                                                <  1:530345/250‑1 (MQ=255)                                                                                                   ctggcgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGGAGGGATTGTCGGGTGGCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGGTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAAAAAAGCATTAATCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCAc                                                                                                                             >  1:126273/1‑251 (MQ=255)                                                                                                            gTCGGTTTGATTCTACCCGGGACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGGTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACGCGATAAAACGGAACATGCGTTGCATCAACACAGCATGTTCACCAGTCTggg                                                                                                                    <  2:261483/251‑2 (MQ=255)                                                                                                                 tttGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGCACGCCTGGCTGGCAGGGGTTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTTTGGCTGGGCTGGCGTTTT‑TAAAAGCCGTTGCCTCGCTGGTCATTTCTG‑GAGAAAAACAAAGCATTACTCGGTAAAAATGAAAATACGTTGCATCAACACAGCAAGTTCCCCCTTCTggtcggtc                                                                                                               >  1:104046/1‑251 (MQ=255)                                                                                                                   tGATTCTTCCCGGTACGGTGCTGATGGGGGGGCTGGGGGCGCTGATTGGCAGCGGCGAGTTAAGTTTCGGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTGCACCATTCTGgtcggtcg                                                                                                              <  2:384902/250‑1 (MQ=255)                                                                                                                                      gCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGt                                                                                                                                                                                                    <  1:322288/146‑1 (MQ=255)                                                                                                                                              gggggCTGGGAGCGCTTGTTGGCTGCGGCGAGTTAAGTTTCTGGGACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTGCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTg                                                                                  <  2:289168/250‑1 (MQ=255)                                                                                                                                              cGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTg                                                                                  <  1:392368/251‑1 (MQ=255)                                                                                                                                                                 tGGCAGCGGCGAGATAAGTTTCTGGCACGCATGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGGG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGTATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATACTTGGCGGGCGCTTTGGTGGCCCGACGCGACCGCTGGTGCCAATGGGGGCGGGGATgg                                                               >  1:483648/1‑250 (MQ=255)                                                                                                                                                                     aGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTc                                                                                                                        <  2:403397/190‑1 (MQ=255)                                                                                                                                                                     aGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTc                                                                                                                        >  1:403397/1‑190 (MQ=255)                                                                                                                                                                        ggcgAGTTAAGTTGCGGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGGTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCt                                                        <  2:340538/251‑1 (MQ=255)                                                                                                                                                                                        ggCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAttc                                        <  2:126273/251‑2 (MQ=255)                                                                                                                                                                                         gCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTAGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAAttc                                        >  1:267865/1‑249 (MQ=255)                                                                                                                                                                                                    tggcagggattgtcggttgtctgctGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGGTCTTCCGGTGGCAAAATTCATTACGCCGAat                            >  1:177059/25‑251 (MQ=255)                                                                                                                                                                                                              ggtcggttgtctgctGGGCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACGCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCGGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGc                   <  1:319090/236‑1 (MQ=255)                                                                                                                                                                                                                             gggCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgccac   >  1:248788/1‑249 (MQ=255)                                                                                                                                                                                                                              ggCG‑ACTGGATCTCCTTCTGGCTGGGCTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCAAAATTCATTACGCCGAATATTATCGGCTGCCTGCTGTGgccacc  <  1:432338/251‑4 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  TTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCG‑ACTGGATTTCTTTCTGGCTGGGTTGGCGTTTT‑AAAAAGCCGTTGCATCGCTGGTCATTTCTG‑AAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAAATTTATTACGCCGAATATTATCGGCTGCCTGCTGTGGCCGCC  >  NC_000913/71337‑71808

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