Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,175,827	C→T	noncoding (74/76 nt)	thrT →	tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,827	0	C	T	98.6%	236.1 / NA	70	noncoding (74/76 nt)	thrT	tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/36);  minor base G (0/1);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GAGACGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGC  >  NC_000913/4175583‑4176066                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                gagtcGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGGAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCCGCATCActtc                                                                                                                                                                                                                                           >  2:322243/5‑251 (MQ=255)  ggaCGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCActtct                                                                                                                                                                                                                                          <  2:158034/250‑1 (MQ=255)      gATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTtc                                                                                                                                                                                                                                      <  2:239722/251‑1 (MQ=255)         aaGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGACGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTctcc                                                                                                                                                                                                                                   <  1:117675/251‑1 (MQ=255)            ttCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCtcctcc                                                                                                                                                                                                                                >  1:39311/1‑251 (MQ=255)            ttCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCCCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCtcctcc                                                                                                                                                                                                                                >  1:14656/1‑251 (MQ=255)               gCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGCGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTg                                                                                                                                                                                                                             >  1:584422/1‑251 (MQ=255)                   aCCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGtttt                                                                                                                                                                                                                         <  1:258581/251‑1 (MQ=255)                   aCCGGATCGAATAAATGGAGCTTATATCGATATTTAGTGCAGAAGGAAAATCAGGGCGCCGAGTTCCAGGATGCGGGCATAGTATAATGGCTACTGCCTCCGCATTCCAAGCTGATGATGCGGGTTCGATTCCCGCCGGCCGATCCAAGATGTGCTGAGATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTTGAATCCGCCTATCAGCATCCCTTCTTTTCTCCTCCCTGttt                                                                                                                                                                                                                          <  2:209969/250‑1 (MQ=255)                      ggCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTc                                                                                                                                                                                                                      >  1:431886/1‑251 (MQ=255)                          cGAATAAAGAGAGCCTCTCTCGATATGCAGAGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATGACCTCAGCCTTCCAAGCTGATGATGCGTGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGGTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTGTTTTCCTTc                                                                                                                                                                                                                  <  2:188727/251‑1 (MQ=255)                           gAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGACGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCt                                                                                                                                                                                                                 >  2:229755/1‑251 (MQ=255)                            aaTAAAGAGCGCTTCTCTCGATATTCATTGCAGAATGATAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTTATCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCTGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTg                                                                                                                                                                                                                <  1:308902/251‑1 (MQ=255)                            aaTAAAGAGAGGTTCTCTCGATATTCAGTGCAGAATGAGAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTg                                                                                                                                                                                                                <  1:308890/251‑1 (MQ=255)                            aaTAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTg                                                                                                                                                                                                                <  1:578999/251‑1 (MQ=255)                               aaaGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGGTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTATAGCGCACCCTTTGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACGTCTTTTCTCCTCCCTGTTTTTTCCTTCTGttt                                                                                                                                                                                                             <  1:499503/251‑1 (MQ=255)                                    gagCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGc                                                                                                                                                                                                        >  1:2360/1‑251 (MQ=255)                                      gCTTATCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTAGCTACCCGCTGCCAGCTCCAAGATGTGCTGATATAGATCAGTTGGGAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTACCTGTTTTTTCCTTCTGTTTATTGCAt                                                                                                                                                                                                      >  1:540650/1‑251 (MQ=255)                                       cTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCAtt                                                                                                                                                                                                     >  1:326741/1‑251 (MQ=255)                                         tctctcGATATTCAGTGCAGAATGAACATCAGGTAGCAGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTCGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTcc                                                                                                                                                                                                   >  2:267759/1‑250 (MQ=255)                                          ctctcGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTca                                                                                                                                                                                                   >  1:465839/1‑250 (MQ=255)                                               gATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGt                                                                                                                                                                                             <  2:502437/251‑1 (MQ=255)                                                     cAGTGCAGAATGAAAATCAGGTTGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGAGGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCAGCACTTCTTTTCGCCCCCCTGTTTGTTCCGTCGGTTTAGTGCATTCAACAAGTCGGGc                                                                                                                                                                                        <  2:258510/250‑1 (MQ=255)                                                                    aTCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                                                                                                                                                                        <  1:229755/251‑1 (MQ=255)                                                                     tCAGGTAGCCGAGTTCCCGGATGCGGGCATCGTCTAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGCTATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtg                                                                                                                                                                       >  1:199124/1‑251 (MQ=255)                                                                        ggTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtgg                                                                                                                                                                    >  1:405135/1‑251 (MQ=255)                                                                         gTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCGATCAGCATCACTTCTCTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTGGCCTGGTTGATgtggt                                                                                                                                                                   <  2:67941/251‑1 (MQ=255)                                                                               gAGTTCCAGGCTGCCGGCCTCGCATAATGGCTATTACCTCAGCCTTCCAAGCTGCTGATGCGGGTTAGATTCCCGATGCCCGCTACCAGCTGTGCTGCTATCGCTCAGTTGGTAGCGCGCACCCTTGCTAGGGGAGAGGACGGCAGTACGCATCTGCCTATCAGCATCACTTCTTTTCCCCTCCCTGTTTTCTCCTTCTGTTTATTGCATTCCACAACTCGGGAATCTTGCCTGGTTACTGTGGTGATATc                                                                                                                                                             >  2:236946/1‑251 (MQ=255)                                                                               gAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATc                                                                                                                                                             <  1:236405/251‑1 (MQ=255)                                                                                        gATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGTGTGAGGTCGGCAGTTCTAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTGCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTTCCGGGTTTATGTGGTGATATCACCGATTTa                                                                                                                                                    <  2:360391/251‑1 (MQ=255)                                                                                           gCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATcc                                                                                                                                                 >  1:397288/1‑251 (MQ=255)                                                                                            cGGGCATCGTATAATGGCTATTACCGCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCg                                                                                                                                                <  2:14656/251‑1 (MQ=255)                                                                                              ggCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgtg                                                                                                                                              <  1:175091/251‑1 (MQ=255)                                                                                                             cTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAAt                                                                                                                               >  2:418196/1‑251 (MQ=255)                                                                                                               aTTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCg                                                                                                                             <  1:485365/251‑1 (MQ=255)                                                                                                                         ccTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTTCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGGGCGCACCCTTGGTAAGGGTGATGTCGGCTGTTCTAATCTGCCTATCAGCATCACTTCTTGTCTCCGCCCTGTGGTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCGTGTTGCCTGGGTGATGTGGTGAGGGCACCGATGTATACGTGTCGTAGAGGGACAATCGATGTCTAAAg                                                                                                                   <  2:335225/251‑1 (MQ=255)                                                                                                                                 gCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGttt                                                                                                           >  2:101744/1‑251 (MQ=255)                                                                                                                                 gCTGAGGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGttt                                                                                                           <  2:431886/251‑1 (MQ=255)                                                                                                                                    gatgatGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTg                                                                                                          <  2:451201/249‑1 (MQ=255)                                                                                                                                     atgatgCGGGTTCGATTCCCGCTGGCCGCTCCAAGATGAGCTGATATAGCTCAGATGGGAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAAAAGGACAATCGATGTCTAAAGAAAAGTTTGAAc                                                                                                       >  2:502824/1‑251 (MQ=255)                                                                                                                                     atgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAAc                                                                                                       <  1:535601/251‑1 (MQ=255)                                                                                                                                         tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTg                                                                                                          >  1:273828/1‑244 (MQ=255)                                                                                                                                         tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGCTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTCAAGAAAAGTTTGAACGTAc                                                                                                   >  1:553333/1‑251 (MQ=255)                                                                                                                                         tgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGCTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTg                                                                                                          <  2:273828/244‑1 (MQ=255)                                                                                                                                          gCGGGTTCGAGTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGGAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAGGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACa                                                                                                  <  1:267759/251‑1 (MQ=255)                                                                                                                                           cGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACaa                                                                                                 <  2:69597/251‑1 (MQ=255)                                                                                                                                                tCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAAcc                                                                                             <  1:240966/250‑1 (MQ=255)                                                                                                                                                 cGATTCCCCCTGCCCGCTCCAAGATGCGCTGATATAGATCAGTTGGCAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTAGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACACGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGCGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGc                                                                                           >  2:552099/1‑251 (MQ=255)                                                                                                                                                  gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGc                                                                                           >  2:154285/1‑250 (MQ=255)                                                                                                                                                       ccGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGtt                                                                                      >  2:205704/1‑250 (MQ=255)                                                                                                                                                          cTGCCCGATCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCACTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTaa                                                                                    >  2:431065/1‑249 (MQ=255)                                                                                                                                                            gCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTc                                                                                >  1:58003/1‑251 (MQ=255)                                                                                                                                                                   ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTAc                                                                           <  1:445136/249‑1 (MQ=255)                                                                                                                                                                   ccAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATGGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTAc                                                                           <  1:358644/249‑1 (MQ=255)                                                                                                                                                                    cAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTAGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGCTGTGGAGATATCACCGATTTATACGTGTCTTAGAGGCACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAAGCGCACGTTAACGTCGGTACtct                                                                        >  2:199589/1‑249 (MQ=255)                                                                                                                                                                            gCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAGTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCAc                                                                 <  1:400727/250‑1 (MQ=255)                                                                                                                                                                            gCTGATATAGCTAAGTTGTTAGAGCGCGCCCTTGGTAAGGGTGATGTCTGCAGTTAGAATCTGCCTATCAGCATCACTGCGTTTCTCCTCCCTGGTTGGTCCTTCTGTTTATCGCATTCAACAAGTCGGGCATGTTGCCGGGTTGATGTGGTGATATCTCCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTGTGAACGTACAAAACCGCACTTTAACTTGGGTACTATCAGCCAc                                                                 <  2:297433/250‑1 (MQ=255)                                                                                                                                                                              tGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTAGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGGTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGTGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTACCGTCGGTACTATCGGCCACGtt                                                              >  1:330579/1‑251 (MQ=255)                                                                                                                                                                                         agctGGTAGAGCGCACCCTTGGTAAGGGTGAGGCCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTa                                                    >  2:306901/4‑250 (MQ=255)                                                                                                                                                                                          gTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTaaa                                                  <  1:279039/251‑1 (MQ=255)                                                                                                                                                                                             ggTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAaac                                                <  1:322243/250‑1 (MQ=255)                                                                                                                                                                                                         ccTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGTTACTATCGGCCACGGTGACCATGGTAAAACAACGCTGACCGCt                                   <  2:344900/251‑1 (MQ=255)                                                                                                                                                                                                             tgTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGGTGACCATGGTAAAACAACGCTGACCGCTGCa                                <  2:39311/249‑1 (MQ=255)                                                                                                                                                                                                                 aGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATccc                           >  2:184788/1‑249 (MQ=255)                                                                                                                                                                                                                   ggTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCATCACTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCTTGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACt                          >  1:34081/1‑250 (MQ=255)                                                                                                                                                                                                                          tCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCg                                     >  2:140354/1‑232 (MQ=255)                                                                                                                                                                                                                          tCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCg                                     <  1:140354/232‑1 (MQ=255)                                                                                                                                                                                                                           cGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTg                 <  1:205704/251‑1 (MQ=255)                                                                                                                                                                                                                                      aTCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCt       >  2:118534/1‑250 (MQ=255)                                                                                                                                                                                                                                          gCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTAcagc  >  1:467476/1‑248 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                GAGACGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGC  >  NC_000913/4175583‑4176066

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