Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,887,917	T→A	6.9%	I25I (ATA→ATT)	rnd ←	ribonuclease D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,887,917	0	T	A	6.9%	194.1 / 4.1	62	I25I (ATA→ATT)	rnd	ribonuclease D
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base A (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGT  >  NC_000913/1887631‑1888204                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                 tGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTa                                                                                                                                                                                                                                                                                                 >  2:352641/1‑288 (MQ=255) tGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTACGCGTACGAACAATTTCAGCATCCAGGGCTa                                                                                                                                                                                                                                                                                                 >  2:225412/1‑288 (MQ=255)       cGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCgg                                                                                                                                                                                                                                                                                          >  1:549894/1‑289 (MQ=255)               cATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCac                                                                                                                                                                                                                                                                                  >  1:521126/1‑289 (MQ=255)                   aCATCGGGGGTCCGCAGAAGGCAACAAGGATTTGCGTGTCAATCAAGGGTTGGGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCacgga                                                                                                                                                                                                                                                                               <  2:414294/288‑1 (MQ=255)                   aCATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCacgga                                                                                                                                                                                                                                                                               <  2:493044/288‑1 (MQ=255)                         ggCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATCAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAACGGCCCGGACGGCtt                                                                                                                                                                                                                                                                         >  2:747168/1‑288 (MQ=255)                           cGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTca                                                                                                                                                                                                                                                                       <  2:421861/288‑1 (MQ=255)                            gTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTcac                                                                                                                                                                                                                                                                      <  2:661749/288‑1 (MQ=255)                            gTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTcac                                                                                                                                                                                                                                                                      <  1:147745/288‑1 (MQ=255)                              ccGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAACTTCAGTATCCAGGGCTATCGCCGGAAAGGCacgga                                                                                                                                                                                                                                                                               >  2:597086/1‑277 (MQ=255)                              ccGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCacgga                                                                                                                                                                                                                                                                               <  1:597086/277‑1 (MQ=255)                               cGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGCTGCCCGCCATCGAAAAGTTGAATCAACCCCAGCTGCCGGTAATAAGTGTGCGTACCAACAAATGCAGTATCCAGGGCCATCGCCGGAAAGGCAAGGAAGGTTCAACACaa                                                                                                                                                                                                                                                                  >  2:18606/1‑289 (MQ=255)                                            aaGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCa                                                                                                                                                                                                                                                     <  1:542103/289‑1 (MQ=255)                                                 gttGCGTGTCAAGCAAGGGGTGGGGGAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTc                                                                                                                                                                                                                                                <  2:179088/288‑1 (MQ=255)                                                   tGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCg                                                                                                                                                                                                                                              <  2:549894/289‑1 (MQ=255)                                                        gtCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGt                                                                                                                                                                                                                                          <  1:638850/288‑1 (MQ=255)                                                         tcaatcaaGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTGCACACAAAGAAGCCAGCGCATCGTCCGTGGTaa                                                                                                                                                                                                                                        >  1:296407/1‑289 (MQ=255)                                                                   gTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGt                                                                                                                                                                                                                               <  2:296407/288‑1 (MQ=255)                                                                     tgtgGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCCGGTAATAAGTGCGCGTACGAACAAATTCAGTCTCCAGGGCTATCGCCGGAAAGGCACGGCCGGCTTCACACAAAAAAGCAAGCGCATACTACGTGGTTATCAATTGGTAAt                                                                                                                                                                                                                            >  1:738179/1‑289 (MQ=255)                                                                       tgGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGACAGGCACGGACGGCTTTACACAAAGAAGCCAGCGCACTGTACGTGGTAATCATTTGGGAAGTc                                                                                                                                                                                                                          >  2:213033/1‑289 (MQ=255)                                                                       tgGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACCCAAAGAAGCCAGCGCATCGGCCGGGGTAATCATTTGGTAAtt                                                                                                                                                                                                                           >  2:638502/1‑288 (MQ=255)                                                                                           tGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGCAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTg                                                                                                                                                                                                      >  1:322356/1‑289 (MQ=255)                                                                                           tGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTg                                                                                                                                                                                                      >  1:215671/1‑289 (MQ=255)                                                                                            gAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGAACGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTgg                                                                                                                                                                                                     >  1:178412/1‑289 (MQ=255)                                                                                                   acTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑Caa                                                                                                                                                                                             >  2:268651/1‑289 (MQ=255)                                                                                                               tCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCgg                                                                                                                                                                                 >  1:536813/1‑289 (MQ=255)                                                                                                               tCACTGCCTGCATGGAGAAATTTCGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACCAAGAAGCCACCGCAACGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCgg                                                                                                                                                                                 >  2:442459/1‑289 (MQ=255)                                                                                                                     cctgcATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAAt                                                                                                                                                                                                                            >  1:365265/1‑241 (MQ=255)                                                                                                                     cctgcATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAAt                                                                                                                                                                                                                            <  2:365265/241‑1 (MQ=255)                                                                                                                               gaAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATAAACCCCAGCTGCGGGTAATAAGTGCGCGTCCGAACAAATTCAGTATCCAGGGCAATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCACAGCATCGTCCGTGGTAATCCTTTGGTAATTGAAATTGTTTTATATTTAGTGGGCGT‑ACAAAAACACGCCGGATTAACCGGCGTCtt                                                                                                                                                                  >  2:543413/1‑287 (MQ=255)                                                                                                                               gaAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATAAACCCCAGCGGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCAATCGCCGGAAAGGAACGGACGGCTTCACACCAAGACGCAAGAGCATCGTCCGTGGTAATCATTTGGTACTTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAAAGCAGGTTATCCCGGGGTCTga                                                                                                                                                                 >  2:543409/1‑289 (MQ=255)                                                                                                                                         gATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACtt                                                                                                                                                       <  2:215671/289‑1 (MQ=255)                                                                                                                                               cGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGc                                                                                                                                                  <  2:35254/288‑1 (MQ=255)                                                                                                                                                ggATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTc                                                                                                                                                <  1:225412/289‑1 (MQ=255)                                                                                                                                                ggATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTc                                                                                                                                                <  1:268651/289‑1 (MQ=255)                                                                                                                                                  aTCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAACGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGCCCGTGGTAATCATTTGGTAATTCTAATTGTTTTCTCTTTGGTTGGCGG‑CAAAAAAAACGCCGTATTAACCGGCCTCTGACGACTGCCTTAACGCTCAg                                                                                                                                              >  2:236268/1‑289 (MQ=255)                                                                                                                                                  aTCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACtt                                                                                                                                                       <  1:7387/280‑1 (MQ=255)                                                                                                                                                  aTCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACtt                                                                                                                                                       >  2:7387/1‑280 (MQ=255)                                                                                                                                                  aTCGCGCAGGAGCGCTTTCAGCGGGGACCAGGCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCCGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCCCGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCAGGTGGTCATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAACAAACGCCGGATTACACGGCGTCTGACCACTGACTTAACGCTCAg                                                                                                                                              >  2:552085/1‑289 (MQ=255)                                                                                                                                                                cTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTcc                                                                                                                                  >  1:592533/1‑287 (MQ=255)                                                                                                                                                                cTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTcc                                                                                                                                  >  2:334476/1‑287 (MQ=255)                                                                                                                                                                           gACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCcgctc                                                                                                                      >  1:591047/1‑286 (MQ=255)                                                                                                                                                                             ccAGTCGGTGATCCCGAGTGGATCGATGAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCCATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAATCGCCGGATTACCCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCt                                                                                                                   <  2:260862/289‑1 (MQ=255)                                                                                                                                                                               aGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCtt                                                                                                                  >  2:454538/1‑288 (MQ=255)                                                                                                                                                                                 tCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGTCCAAAAAAAACGCGGGTTTAACCGCCGCTTGAGGATTG‑CCTAACGGTCAAGCCTTAttcgcccatttttctgccgccttc                                                                                                                >  1:632357/1‑266 (MQ=255)                                                                                                                                                                                  cGGTGATCCCGAGTGGTTCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAACCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCg                                                                                                                <  1:395043/287‑1 (MQ=255)                                                                                                                                                                                              gTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCtcg                                                                                                  >  2:64786/1‑289 (MQ=255)                                                                                                                                                                                                                             gAAAAGTGGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGt                                                                   <  1:82555/289‑1 (MQ=255)                                                                                                                                                                                                                                         cAACCCCAGCGGCTGGTAATAAGTGCCGGTACGAACAAATTCAGTATCCAGGGCTATCGCCAGCAAGGAACGGACGGCTTAACAGAAAGGAGCCAGCGTATCGTCCGTGGTAATCATTCGGTAATTCAAATTGTTTTCTCTTTAGGCAGCGT‑CAAAACAAACGCCGGATTAACCGGCGTCTGACGACTTACTTAACGCTCAGGCTTTATTGTCCGCTTTGCCGCGCGCTTCGTCCCGTAACTCTCGTCGCAAAATTTTTCCGCCGTTAGAGTTCGGTAACTCATCACGa                                                       <  2:669546/289‑1 (MQ=255)                                                                                                                                                                                                                                               cAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTcc                                                 <  1:638502/289‑1 (MQ=255)                                                                                                                                                                                                                                                 gCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTcca                                                >  1:23558/1‑288 (MQ=255)                                                                                                                                                                                                                                                   tgagggTAAAAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACAAAAAGAAGCAAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                                              <  1:554540/285‑1 (MQ=255)                                                                                                                                                                                                                                                   tGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTccacc                                              <  1:352641/288‑1 (MQ=255)                                                                                                                                                                                                                                                       ggTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCt                                          >  2:197186/1‑288 (MQ=255)                                                                                                                                                                                                                                                          aataaGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCgg                                      <  1:197186/289‑1 (MQ=255)                                                                                                                                                                                                                                                           ataagttggccgaCGAACAAATTAAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCAAACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGt                                     <  2:322356/277‑1 (MQ=255)                                                                                                                                                                                                                                                           ataaGTGCGCGTACGAACAAATTCAGTATCAAGGGCAATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGt                                     <  1:543413/289‑1 (MQ=255)                                                                                                                                                                                                                                                           ataaGGGCGCGTACGAACAAATTCATTATCCAGGTCAATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCATCTTCGGt                                     <  1:543409/289‑1 (MQ=255)                                                                                                                                                                                                                                                                   gcgTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGt                              <  2:592533/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                ttCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGa                <  1:64786/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                       tCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACggcggc         <  1:213033/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                               tATCGCCGGAAAGGCACGGACGGCTTCACCCAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGt  <  1:767335/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                               tATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGt  <  1:454538/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                 TGGAAGCGAAACCCCATGACATCGGGCGTCCGCAGAAGGCAGCAAGGATTTGCGTGTCAATCAAGGGTTGTGGTAATTCGCCAAAGACATTGAGGAACACTTCCAGATCTTCACTGCCTGCATGGAGAAATTTTGTGATGGACGGATCGCGCAGGATCGCTTTCAGCGGTGACCAGTCGGTGATCCCGAGTGGATCGATTAGCGCCAGATGCTCGCCATCGAAAAGTTGAATCAACCCCAGCTGCGGGTAATAAGTGCGCGTACGAACAAATTCAGTATCCAGGGCTATCGCCGGAAAGGCACGGACGGCTTCACACAAAGAAGCCAGCGCATCGTCCGTGGTAATCATTTGGTAATTCAAATTGTTTTCTCTTTAGTGGGCGT‑CAAAAAAAACGCCGGATTAACCGGCGTCTGACGACTGACTTAACGCTCAGGCTTTATTGTCCACTTTGCCGCGCGCTTCGTCACGTAATTCTCGTCGCAAAATTTTTCCGACGTTAGATTTCGGTAACTCATCACGAAACTCCACCAGCTTCGGTACTTTGTATCCCGTGAGCTGACGGCGGCAAAAAGT  >  NC_000913/1887631‑1888204

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