Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	92,810	T→C	T466T (ACT→ACC)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	92,810	0	T	C	95.1%	175.1 / 1.0	61	T466T (ACT→ACC)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (32/26);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

GACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGC  >  NC_000913/92563‑93052                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      gACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGTGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTAGATTACCAAAGTTGACCCACCGCGTCCCGGTGCACGTGTCTTCCCGGAACCCATTGGCCGCACggtg                                                                                                                                                                                                                                                    >  2:43160/1‑251 (MQ=255) gACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGCTTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACcgtg                                                                                                                                                                                                                                                    >  2:75782/1‑251 (MQ=255)             cGACCAATTTGGGGTTGGTCGGCGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGACCGCACCGTGGGGCATatgatt                                                                                                                                                                                                                                        >  2:347006/1‑250 (MQ=255)             cGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGTGCATatgatg                                                                                                                                                                                                                                        >  2:129521/1‑251 (MQ=255)                ccAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCAAATGATGGa                                                                                                                                                                                                                                      >  2:361780/1‑250 (MQ=255)                 cAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACAGAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGaaa                                                                                                                                                                                                                                    <  2:351896/251‑1 (MQ=255)                  aaTTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGaa                                                                                                                                                                                                                                     >  1:174636/1‑249 (MQ=255)                       ggggTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTg                                                                                                                                                                                                                               <  1:129521/250‑1 (MQ=255)                            tGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCt                                                                                                                                                                                                                          <  1:449784/250‑1 (MQ=255)                                cGGAGAACGCAGTGGCTTATGTCCTCAAAAATAACGGTGGTCTGACATGGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGTGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTAcc                                                                                                                                                                                                                       >  1:146391/1‑249 (MQ=255)                                     aaCGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTAAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGgg‑cggc                                                                                                                                                                                                               >  2:414171/1‑251 (MQ=255)                                        gaagtagCTTATGTCCTCAACAACAACGGTGGACTGACATAGAGAGGGCCACATTCTCTTTCGGCTACGAGCTAATGGTCACACCGTTACAGTTAGAGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCCAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGACCGCACCGAGGTGCACATGATGGCAAACGAGGCGCGACCAgcgggggac                                                                                                                                                                                                             >  2:452830/7‑243 (MQ=255)                                                 tataTCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGATACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGa                                                                                                                                                                                                     >  1:113239/1‑249 (MQ=255)                                                      ccTCAAAAACAACGGTGGTCTGAAATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAggcgg                                                                                                                                                                                               <  1:96474/250‑1 (MQ=255)                                                       cacAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCGTCCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCGCCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGgcggcg                                                                                                                                                                                             <  2:158866/249‑1 (MQ=255)                                                           aaaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTa                                                                                                                                                                                         <  1:41040/251‑1 (MQ=255)                                                            aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCCACTATCGGCAGCTACGGCCTTTATCGCCCCCTGTCGCTTACCAAAGTTCACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGGGCATTTGATGGAAAGCGGGGCGCTACCGGG‑ATGCGGCGGCGTGAcggcggggctta                                                                                                                                                                                         >  2:90969/1‑244 (MQ=255)                                                                 acGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGACCGCACCGTGGTGCCTATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCt                                                                                                                                                                                   >  2:446989/1‑251 (MQ=255)                                                                      ggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgt                                                                                                                                                                              >  1:364789/1‑251 (MQ=255)                                                                                  atagGGCCACCCTCTCTTTCGGCTACGGGCTAATGGTAACCCCGTTACAGTTAGCGCGAGTCTACGCAACTATCGTCAGCTACGGCATTTAGCGCCCACTGGCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaa                                                                                                                                                                  <  2:352400/249‑1 (MQ=255)                                                                                   gagGGCCACCTTCTCTTTCGGCTACGCGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTGCCCGGTGAACGTGTCTGCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaaa                                                                                                                                                                 <  1:19411/251‑1 (MQ=255)                                                                                       gCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTATACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGCCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCgg                                                                                                                                                             >  2:406496/1‑251 (MQ=255)                                                                                        ccACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGAGAGTCTAAGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCACGGTGAAAGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGAGCCTATGATGGAAAGCCTGGCGCT‑CCAGGGAGGCGGCGGCGTGAAGGCGGCGATTAAAGGAGACCGTATCACCCTTAAAAcagg                                                                                                                                                             >  2:133663/1‑247 (MQ=255)                                                                                          cccTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAc                                                                                                                                                          <  1:414171/250‑1 (MQ=255)                                                                                                       gCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTc                                                                                                                                           <  2:362128/251‑1 (MQ=255)                                                                                                        cTACGGGCTACTGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAGGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAGAGGCTATCGTATCGCCCTTAAAACCGGTAC‑C‑GCGAAACAGGTCg                                                                                                                                          >  2:4484/1‑251 (MQ=255)                                                                                                        cTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCg                                                                                                                                          >  1:314107/1‑251 (MQ=255)                                                                                                            gggCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGcc                                                                                                                                      <  1:363148/251‑1 (MQ=255)                                                                                                            gggCTAAGTGTAACACCGGTACTGTTAGCGCGAGGCTACGCAACTATCGTCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGcc                                                                                                                                      <  2:20532/251‑1 (MQ=255)                                                                                                              gCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTCTCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCACGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGAATAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑A‑GCGAAAAAGGTCGGGCCgg                                                                                                                                    >  1:62683/1‑251 (MQ=255)                                                                                                                tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGCTTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGGGGTGCATATGATGGAAAGAGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTAGCGTATCGCCATTAAAACCGG‑AC‑CAGCGAAAAAGGGCGGGCCGGAc                                                                                                                                  >  2:1803/1‑251 (MQ=255)                                                                                                                tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGCATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTCAAACCGGTAC‑C‑GCGGAACAGGGCGGGCCgtcc                                                                                                                                  >  2:206291/1‑248 (MQ=255)                                                                                                                   tGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCACACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGGGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTGCCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGCAC‑C‑GCGAAAAAgggcgggcagggaggg                                                                                                                               >  2:175186/1‑243 (MQ=255)                                                                                                                       aaCACCGTTTCAGTTAGCGCGAGTCTACGCAACTACCGGCAGCTACGGCATTTATCGCCCACTGTCGATTTCCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCt                                                                                                                           <  2:354930/251‑1 (MQ=255)                                                                                                                       aaCACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGc                                                                                                                            <  1:160286/250‑1 (MQ=255)                                                                                                                        acacCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAATTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCTGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTAC‑C‑GCGAAAAAGGTAGGGCCCGACgggccctc                                                                                                                          >  2:179159/1‑245 (MQ=255)                                                                                                                         caccgtttcaggtAGCGCGAGTCTACGCACCTATCGGCAGCTACGGCAGTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTGCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGGGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGTTCGGGCCGGACGGGCGCTAc                                                                                                                         <  2:416168/239‑1 (MQ=255)                                                                                                                           ccgttACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACAt                                                                                                                       <  1:203751/248‑1 (MQ=255)                                                                                                                                   gTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGAATTTATCGCACACTGGCGATTACCACAGTTGGCCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAATGCGTGGCGCTACAGGG‑CGGAGACGGCGTGGAGGCGGCGATTAAGGGCTATCGTATCGCCATTAGAACCGGTACAC‑G‑GAAAAAGGTCGGAACGGAAGGTCGGTACCTCaataaat                                                                                                               >  2:406084/1‑251 (MQ=255)                                                                                                                                     tAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGGCGGGCCGGACGGTCGCTACATCAATAAatat                                                                                                             >  2:98664/1‑251 (MQ=255)                                                                                                                                                       cTATCGGCAGCTACGGCATTTATCGACCACTGTCGATTACCACAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCTCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGc                                                                                           >  1:133950/1‑251 (MQ=255)                                                                                                                                                                 cTACGGCATTTATCGCCCACTGTCGATTCCCAAAGTTGACCCCCCGGTTCCCGGTGCACGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGTGCGTATGATGGAAAGCGTGGCGCTACCGGA‑CGGCGCCGGCGTGAAGGCGGCGATTAACGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGACCACGGTCGAGCCGGACGGTCGTTACATCATTATAAACATTGATTATACCGCAGGCGTTGCGCCTg                                                                                 >  2:107860/1‑251 (MQ=255)                                                                                                                                                                  tACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCg                                                                                                                                          <  1:38795/193‑1 (MQ=255)                                                                                                                                                                  tACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCg                                                                                                                                          >  2:38795/1‑193 (MQ=255)                                                                                                                                                                       cATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCTCTACCGGT‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCAGTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTc                                                                           <  1:305758/251‑1 (MQ=255)                                                                                                                                                                                   aCTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAAACCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTcg                                                               >  2:468299/1‑251 (MQ=255)                                                                                                                                                                                    atGTCGATTACCACAGTTGACCCCCCGGATCCCGCGGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCACATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTACAGGCTATCGAATCGCCATTAGAACCGGAAC‑C‑GCGTGACGGGTAGGGCCCGACGGGGGCTACATCCGTAAATATATAGGTTAAAACGCAGGCCATGCGccgtggggtccacagcacgtcac                                                              >  2:259890/2‑230 (MQ=255)                                                                                                                                                                                      gTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCgcgc                                                            <  1:243068/251‑1 (MQ=255)                                                                                                                                                                                                aaaGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGgttgttgt                                                  >  1:436574/1‑251 (MQ=255)                                                                                                                                                                                                           cccGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATc                                       <  1:361780/251‑1 (MQ=255)                                                                                                                                                                                                            ccGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATcc                                      >  1:144835/1‑251 (MQ=255)                                                                                                                                                                                                             cGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTTTCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTTCTTATCCCGCAGGCGTTGCGCCTGCGCGTCCGCCGCGCATCGCGCTGGTTGTTGTTATCAACGATCCg                                     >  2:111370/1‑251 (MQ=255)                                                                                                                                                                                                                 tCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAgg                                 >  1:283990/1‑251 (MQ=255)                                                                                                                                                                                                                     ggTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGGGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCggg                             <  1:4484/251‑1 (MQ=255)                                                                                                                                                                                                                              gtCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                      <  2:24739/185‑1 (MQ=255)                                                                                                                                                                                                                              gtCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTgc                                                                                      >  1:24739/1‑185 (MQ=255)                                                                                                                                                                                                                                  tCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGTCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAtactac                   <  1:198396/248‑1 (MQ=255)                                                                                                                                                                                                                                     cGGAAGCCAGTGGCCGCACCGTGGGGAATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGTCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTAGACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACggcggc             <  2:174636/251‑1 (MQ=255)                                                                                                                                                                                                                                       gAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTAGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAATGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGACAGCCGCGCTTCGCGCTGGTTGTTGTTTTCAACGATCCGCAGGCggggacataccaacgccgggg           >  2:338162/1‑233 (MQ=255)                                                                                                                                                                                                                                        aaTCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGcc          <  1:206291/251‑1 (MQ=255)                                                                                                                                                                                                                                            cATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGACTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTaaaatcttcggaggagcccgtt      >  1:19009/1‑232 (MQ=255)                                                                                                                                                                                                                                                gTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAAttataaggagggggcgttttccg  >  1:405883/1‑229 (MQ=255)                                                                                                                                                                                                                                                gTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCcgc  <  1:394048/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      GACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCCCCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAC‑C‑GCGAAAAAGGTCGGGCCGGACGGTCGCTACATCAATAAATATATTGCTTATACCGCAGGCGTTGCGCCTGCGAGTCAGCCGCGCTTCGCGCTGGTTGTTGTTATCAACGATCCGCAGGCGGGTAAATACTACGGCGGCGCCGTTTCCGC  >  NC_000913/92563‑93052

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