Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,884,734	C→T	P81L (CCG→CTG)	rnpA →	RNase P protein component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,884,734	0	C	T	98.5%	219.7 / NA	67	P81L (CCG→CTG)	rnpA	RNase P protein component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/31);  minor base G (0/1);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTGGAGTGATAAAAGGCAGT  >  NC_000913/3884487‑3884977                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    cccTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTTTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACAGGAACTCCTGGc                                                                                                                                                                                                                                                  <  2:438982/251‑1 (MQ=255)  ccTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCCGTAAATTCACATTCGTCTTCCAGCAGCAACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAACCGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCt                                                                                                                                                                                                                                                 >  1:117438/1‑251 (MQ=255)  ccTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATCCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCGTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGTATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCt                                                                                                                                                                                                                                                 <  1:293898/251‑1 (MQ=255)   cTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCt                                                                                                                                                                                                                                                 >  1:140629/1‑250 (MQ=255)   cTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCt                                                                                                                                                                                                                                                 >  1:567806/1‑250 (MQ=255)     gAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATg                                                                                                                                                                                                                                              >  1:555042/1‑251 (MQ=255)      aGTGGGTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTACCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGGCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGTTTAAACGTCTGCCGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCATGGCTATgg                                                                                                                                                                                                                                             >  2:441947/1‑251 (MQ=255)          gTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGAttt                                                                                                                                                                                                                                         <  1:297757/251‑1 (MQ=255)                   gCATTTCCCAGGGAGTTCCGCTTGTTATCTCTCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCgtggtggt                                                                                                                                                                                                                                <  2:472629/251‑1 (MQ=255)                          ccAGGGAGTTACGCTTGTTAACTCCCCGTCAATTCACATTCGTCTTCCAGCCGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGaa                                                                                                                                                                                                                          >  2:415340/1‑250 (MQ=255)                           cAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGaaaa                                                                                                                                                                                                                        <  1:287178/251‑1 (MQ=255)                                    ccGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTg                                                                                                                                                                                                               <  2:211908/250‑1 (MQ=255)                                    aCGCTTGTTAACCCCCATCCAATTCACATTCGTCTTCCAGCGGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTg                                                                                                                                                                                                               <  2:94382/251‑1 (MQ=255)                                     cGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCAAAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGc                                                                                                                                                                                                              >  1:340634/1‑251 (MQ=255)                                              aCTCCCAGTCAATTCACATTCGTGTTCCAGCAGCCACAACGGGCTGGCACACCGCAAATTACCATTCTCGGCTGCCTGAATTCCCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGCACGCACTCGGATTAAACGTATGACGCGTGAAAGCTTCCGTCTGCGACAACATAAACTCCTGGCTATGGATTTCGTGGAGGTGGCGAAAAAAGGGGTTGCCGACCTCg                                                                                                                                                                                                      >  2:67228/1‑250 (MQ=255)                                                ccccAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATa                                                                                                                                                                                                   <  2:577717/250‑1 (MQ=255)                                                     gTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGACACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGATATGGATTTCGTGGTGGTGGCGGAAAAAGGGGTTGCCGACCTCGATAACCGt                                                                                                                                                                                              >  1:563830/1‑251 (MQ=255)                                                       cAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATg                                                                                                                                                                                                                                              <  2:552241/201‑1 (MQ=255)                                                       cAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACCCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATg                                                                                                                                                                                                                                              >  1:552241/1‑201 (MQ=255)                                                           tCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGctctc                                                                                                                                                                                        >  1:463325/1‑251 (MQ=255)                                                               aTTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCgg                                                                                                                                                                                    >  2:541461/1‑251 (MQ=255)                                                               aTTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCgg                                                                                                                                                                                    >  2:241947/1‑251 (MQ=255)                                                                       tCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGt                                                                                                                                                                              <  2:612315/249‑1 (MQ=255)                                                                               ccACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAAtt                                                                                                                                                                    >  1:167432/1‑251 (MQ=255)                                                                                       ggCTGGCACGCGGCAAATTACCATTCTCGGCTGCCTGACTTAGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAACCGCTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTATGCGCCAACATGAACTCCTGGCTATGGCTTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGcgcc                                                                                                                                                            >  2:255571/1‑251 (MQ=255)                                                                                           ggCACGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTGCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCa                                                                                                                                                        <  2:167432/251‑1 (MQ=255)                                                                                              aCGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTg                                                                                                                                                     >  1:465128/1‑251 (MQ=255)                                                                                              aCGCCGCAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTg                                                                                                                                                     >  1:151455/1‑251 (MQ=255)                                                                                                 ccgcAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCg                                                                                                                                                  <  2:299964/251‑1 (MQ=255)                                                                                                 ccgcAAATTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCg                                                                                                                                                  <  1:73085/251‑1 (MQ=255)                                                                                                     aaaTTACCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTg                                                                                                                                              >  2:287604/1‑251 (MQ=255)                                                                                                          aCCATTCTCGGCTGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTc                                                                                                                                          <  1:8470/250‑1 (MQ=255)                                                                                                                 ccggctgcttGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCGGCGCCAACATGAACTCCTGCCTATGGATTGCGTGGTGGTGGCGAAAAAAGGGGTTCCCGCCCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTcc                                                                                                                                  <  2:129731/242‑1 (MQ=255)                                                                                                                          tGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGccc                                                                                                                         >  2:237988/1‑251 (MQ=255)                                                                                                                                cGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAACCGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTc                                                                                                                   <  1:604016/251‑1 (MQ=255)                                                                                                                                 gCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCg                                                                                                                  >  1:357900/1‑251 (MQ=255)                                                                                                                                                      tCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGCCTATCACCGCCTGAtt                                                                                               >  2:376480/1‑249 (MQ=255)                                                                                                                                                          aCAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCTGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTcc                                                                                          <  1:102087/250‑1 (MQ=255)                                                                                                                                                          aCAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTcc                                                                                          <  1:443303/250‑1 (MQ=255)                                                                                                                                                             gTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCt                                                                                       >  1:590006/1‑250 (MQ=255)                                                                                                                                                                gCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCt                                                                                                                                                                                                        >  1:260402/1‑134 (MQ=255)                                                                                                                                                                gCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCt                                                                                                                                                                                                        <  2:260402/134‑1 (MQ=255)                                                                                                                                                                     gAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGcc                                                                              <  1:217248/251‑1 (MQ=255)                                                                                                                                                                      aaaaaCTTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCGGGCTTTGGATTTCGTGGTGGTGGGGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCg                                                                             <  2:555042/251‑1 (MQ=255)                                                                                                                                                                          aCGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAgggg                                                                                                                                                                                                                  >  1:372861/1‑114 (MQ=255)                                                                                                                                                                          aCGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAgggg                                                                                                                                                                                                                  <  2:372861/114‑1 (MQ=255)                                                                                                                                                                          aCGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCAtt                                                                         >  1:127401/1‑251 (MQ=255)                                                                                                                                                                                     gcCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACt                                                              >  1:228092/1‑251 (MQ=255)                                                                                                                                                                                             aCGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCTTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTg                                                      >  1:534621/1‑251 (MQ=255)                                                                                                                                                                                              cGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGACGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTTTCGTTTCACTCCAACCTGt                                                     >  2:548379/1‑251 (MQ=255)                                                                                                                                                                                               acAATCGGACTACACGACTGACGCGTAAAAGGTGCCGTTTGCCCCAACAGGAACGCCTGGCTATGGATTTGGTGGGGGGGGTGAAAAAAGGGGTGGCAGTCTTCGATAACCGTGCGCCATCGACAGCGTTGGAAAAATTATGCCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGGTTCACTCCAACCTGtt                                                    <  2:45248/250‑1 (MQ=255)                                                                                                                                                                                                cAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTc                                                   <  2:357900/251‑1 (MQ=255)                                                                                                                                                                                                         tAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACg                                           <  1:241947/250‑1 (MQ=255)                                                                                                                                                                                                                     gcgTGAAAGCTGCCGTCTGCGCAAACATGAACTCCTGGCTATGGATTTCGTGGTGGGGGCGAAAAAAGGGGTTGCCGACCTCGATACCCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCAtt                              <  2:373276/251‑1 (MQ=255)                                                                                                                                                                                                                      cgTGAAAGCTTCCGTCTGCGCCAACATGAAATCCTGGCTATGGATTTCGTGGGGGTGGCGAAAAAAGGGGTTGCCGACCTAGATAACCATGCTCTCTCGGAAGCATTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTACTGATAGCCCTCATTCGGGTCTATCAACGCCGGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTg                             >  2:581852/1‑251 (MQ=255)                                                                                                                                                                                                                        tGAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTgcg                           >  1:145098/1‑251 (MQ=255)                                                                                                                                                                                                                         gAAAGCTTCCGTCTGCGCCAACATGAACTCCTGGCTATGGCTTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTgcgc                          >  1:592270/1‑251 (MQ=255)                                                                                                                                                                                                                         gAAAGCGTCCGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGGGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTAGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTACGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTgcgc                          >  2:205956/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccGTCTGCGCCAACATGACCTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTgg                  <  1:136709/251‑1 (MQ=255)                                                                                                                                                                                                                                  cGTCTGCGCCAACATGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTggg                 >  2:481239/1‑250 (MQ=255)                                                                                                                                                                                                                                    gcagcgcCAACATGAACTCCTGGCTATGGATTTCGTGGGGGTTGGGAAAAAAGGGGTGGCCGACGTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTAGGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCGGATTAGGCCGCTACTCGGGCCGCATGGTCGTGTCACTGCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTGGAGt               <  1:613242/248‑1 (MQ=255)                                                                                                                                                                                                                                       gcgcGAACTTGAACTCCTGGATATGGATTGCATAGTGGTGCCGAAGAAAGGGTTATTCGGCCTCGATAACCTTGCTGTGTCGGACCCTTTGGAAAAATTAGGTCTACGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGCTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGGCGTTTCACTCCAGCCTTTTAAAGCTACGGGATGGAGGCATTGCGCAGGTTTGGAGTGa             <  2:632542/250‑1 (MQ=255)                                                                                                                                                                                                                                          ccAACATGAACTCCTGGCTATGGATTTCGGGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATACCCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTGGAGTGATaaa         <  1:560600/251‑1 (MQ=255)                                                                                                                                                                                                                                           caacaaGAGCTCCTGGCGATGGATTTCGGGGTGGTGGCGAAAAAAGGGGTTGCGGCCCTCGATAACCGTGCTCGCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGCTTTATCAACGCCTGATTAGGCCGCTACTCGGGCCGCATTGTCGTGTCACTCCAACCTGTTCAAGCTACGGGATTGAGGCATTGCGCAGGTTTGGAGTGATaaaa        <  1:226770/251‑1 (MQ=255)                                                                                                                                                                                                                                           cAACATGAACTCCTGGCTATGGTTTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAGCCTGTTCAAGCTACGGAATTGAGGCATTGCGCGGGTTTGGAGTGATaaaa        <  2:79253/251‑1 (MQ=255)                                                                                                                                                                                                                                                tGAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTGGAGTGATAAAAGGc     >  1:14639/1‑249 (MQ=255)                                                                                                                                                                                                                                                 gAACTCCTGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTGGAGTGATAAAAGGCAgt  >  1:247099/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CCCTGAGTGGTTAAGCTCGCATTTCCCAGGGAGTTACGCTTGTTAACTCCCAGTCAATTCACATTCGTCTTCCAGCAGCCACAACGGGCTGGCACGCCGCAAATTACCATTCTCGGCCGCCTGAATTCGCTGGGGCATCCCCGTATCGGTCTTACAGTCGCCAAGAAAAACGTTCGACGCGCCCATGAACGCAATCGGATTAAACGTCTGACGCGTGAAAGCTTCCGTCTGCGCCAACATGAACTCCCGGCTATGGATTTCGTGGTGGTGGCGAAAAAAGGGGTTGCCGACCTCGATAACCGTGCTCTCTCGGAAGCGTTGGAAAAATTATGGCGCCGCCACTGTCGCCTGGCTCGCGGGTCCTGATAGCCCTCATTCGGGTCTATCAACGCCTGATTAGTCCGCTACTCGGGCCGCATTGTCGTTTCACTCCAACCTGTTCAAGCTACGGAATTGAGGCATTGCGCAGGTTTGGAGTGATAAAAGGCAGT  >  NC_000913/3884487‑3884977

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