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	95.8%	148.3 / NA	48	noncoding (74/76 nt)	thrT	tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/28);  minor base G (0/2);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCT  >  NC_000913/4175587‑4176072                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      cGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCtttt                                                                                                                                                                                                                                             <  1:331411/251‑1 (MQ=255)           gctaaCCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAATATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTg                                                                                                                                                                                                                                   <  2:350102/247‑1 (MQ=255)            cTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTtc                                                                                                                                                                                                                                            <  2:320066/241‑1 (MQ=255)            cTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTtc                                                                                                                                                                                                                                            >  1:320066/1‑241 (MQ=255)                      cGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCtt                                                                                                                                                                                                                         <  1:361528/250‑1 (MQ=255)                      cGAATAAAGAGAGCTTCTCTCGAGATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCtt                                                                                                                                                                                                                         <  1:384531/250‑1 (MQ=255)                        aaTAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTATCCTCCCTGTTTTTTCCTTCTg                                                                                                                                                                                                                      >  1:339809/1‑251 (MQ=255)                             agagagCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTGGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAt                                                                                                                                                                                                                 <  1:110487/251‑1 (MQ=255)                                            atatTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAAGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTc                                                                                                                                                                                                  >  1:371397/1‑251 (MQ=255)                                                                    ggTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtgg                                                                                                                                                                          <  1:291435/251‑1 (MQ=255)                                                                    ggTAGCCGAGTTCCAGGATGCGGGCATCGCATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGGGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtgagg                                                                                                                                                                          >  2:299698/1‑248 (MQ=255)                                                                       aGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCGATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGTCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTAATGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGa                                                                                                                                                                       <  2:8438/251‑1 (MQ=255)                                                                               tCCAGGATGCGGGCATCGTATAATGGCTATTACCTCATCCTTGCAAGCTGATGATGCGGGTTCGATTCCCGCGGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTATAGCGCACCCTTGGTAAGTGTGAGGTCGGCAGTTCGAATCTGCCGATCAGCAGCACGTCGTTTCTCCTCCCTGTTTTTTCCTGCTGTTTATTGCATTCAACAAGTCGGGCTTGTTGCCTGGTTTATGTCGTGATATCAcc                                                                                                                                                                <  2:27695/250‑1 (MQ=255)                                                                                          ggCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGCTGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCgtg                                                                                                                                                    >  2:378213/1‑251 (MQ=255)                                                                                               cGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTGATCCGTGTCTTa                                                                                                                                               <  2:45814/251‑1 (MQ=255)                                                                                               cGTATAATGGCTATTACCGCAGCCTTCCAAGCTGAGGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCATCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTa                                                                                                                                               <  1:12174/251‑1 (MQ=255)                                                                                                        gCTATTACCTCAGCCTTCCAAGCTGATGATGCGTGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCTCCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCAGCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCGGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGATGGACa                                                                                                                                       <  1:343678/250‑1 (MQ=255)                                                                                                           aTTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATc                                                                                                                                    <  1:165056/250‑1 (MQ=255)                                                                                                           aTTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATc                                                                                                                                    <  1:401106/250‑1 (MQ=255)                                                                                                           aTTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATAGCACCGATTTATCCGTGTCTTAGAGGGACAATc                                                                                                                                    <  2:104405/250‑1 (MQ=255)                                                                                                              aCCTCAGCGTTTGAAGCTTATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGTTATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGACTCTGCCTATCAGCATCActt                                                                                                                                                                                                                                                  >  2:353245/1‑137 (MQ=25)                                                                                                              aCCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGCCGGCAGTTCGAATCTGCCTATCAGCATCActt                                                                                                                                                                                                                                                  <  1:353245/137‑1 (MQ=255)                                                                                                                  cAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTa                                                                                                                            <  1:25774/251‑1 (MQ=255)                                                                                                                         ccAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGGTATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGaaa                                                                                                                      >  1:244358/1‑250 (MQ=255)                                                                                                                               tgatgaTGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGAAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCAt                                                                                                                                                                                            >  1:383869/1‑174 (MQ=255)                                                                                                                               tgatgaTGCGGGCTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCCCCCTGGGTAAGGGTGAGGCCGGCAGTTCGAATCTGTCTAGCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGGTTATGGCATTCATCAAGTCGGGCAt                                                                                                                                                                                            <  2:383869/174‑1 (MQ=255)                                                                                                                                gatgatGCTGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGGCTAAAGAAAAGTTTGaa                                                                                                              <  2:137781/251‑1 (MQ=255)                                                                                                                                gatgatGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGa                                                                                                               >  2:292236/1‑250 (MQ=255)                                                                                                                                  tgatgCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGTTGTCTAAAGAAAAGTTTGAACg                                                                                                            >  2:320545/1‑251 (MQ=255)                                                                                                                                         ggTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACaaa                                                                                                      >  2:128159/1‑250 (MQ=255)                                                                                                                                           ttCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCCGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAAcc                                                                                                   >  2:322670/1‑251 (MQ=255)                                                                                                                                            tCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAAcc                                                                                                   <  1:98149/250‑1 (MQ=255)                                                                                                                                              gATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGCAGAGCGCACCCTTGGTAAGGGCGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGACCGTACAAAACCGc                                                                                                 >  1:241944/1‑250 (MQ=255)                                                                                                                                               aTTCCCGCAGCCCGCTCCAAGATGTGCTGCTATCGCTCCGTTGGTAGAGCGCACCCTTGGTAAGGGGGAGGCCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCACCAAGTCGGGCATGTTGCCTGGTTGCTGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGCTGTCTAAAGAAAAGTTTGAACGTACAAAACAGCAc                                                                                               >  2:374977/1‑251 (MQ=255)                                                                                                                                                   ccgttgctcgatCCAAGATGTGCTGAAATAGATCAGTTGGTATAGCGCTCCCTTGGTAAGGGTGAGGTCGGAAGTTCGAAGCTGCCTATACGCATCACTTCTTTACTCCTACCTGTTTTTACCTGCTGTTTATTGCATTCAACAAGTCGTGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCTTGTCTTAGAGGGACAATCGATCTCTAAAGAAAAGTTTGAACGTACAAAACCGAACGtga                                                                                           <  2:66705/240‑3 (MQ=255)                                                                                                                                                        gccggCTCCAAGATGGGCTGAGATAGCTCAGTTGGTAGAGCGCACCCCTGGTAAGTGTGATGTCGGCAGTTCGAATCTGTCTATCAGCATCACGTCTTTTCTCCTTCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGGCGGGGATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTATAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACACAACAGCACGTTAACGTc                                                                                      <  2:200546/247‑1 (MQ=255)                                                                                                                                                                            atatAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGa                                                                  <  1:128159/251‑1 (MQ=255)                                                                                                                                                                              ataGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGAc                                                                 <  2:74281/250‑1 (MQ=255)                                                                                                                                                                                     aGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTa                                                                                                                            >  1:87311/1‑184 (MQ=255)                                                                                                                                                                                     aGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTa                                                                                                                            <  2:87311/184‑1 (MQ=255)                                                                                                                                                                                                         ggTAAGGGTGAGACCGGCAGTTCGAATCTGCCTATCTGCATCTCTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAATTCGGGCATTTTGCCTGGTTGATGTGGTGATATCACTGAGTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTTCAAAACCGCACTTTAACGTCGGTTCTATCGGCCTCGTTGACCATGGTAAAACAACGCTGTCCGCTGCa                                      <  1:339584/250‑1 (MQ=255)                                                                                                                                                                                                            aaGGGGGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTTCAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATc                                   <  2:232098/250‑1 (MQ=255)                                                                                                                                                                                                             aGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCAc                                 >  1:100684/1‑251 (MQ=255)                                                                                                                                                                                                                                 aaTCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTgagatccgga              >  2:213989/1‑241 (MQ=255)                                                                                                                                                                                                                                 aaTCTGCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTg                       <  1:213989/241‑1 (MQ=255)                                                                                                                                                                                                                                      gCCTATCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTAcgg         <  1:292236/250‑1 (MQ=255)                                                                                                                                                                                                                                          aTCAGCATCACTTCTTTTCTCCTCCCTGTTTTTTTCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGtg    <  2:133398/251‑1 (MQ=255)                                                                                                                                                                                                                                            cAGCATCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGtgct  >  1:127938/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      CGATAAGTTCGCTTACCGGCTCGAATAAAGAGAGCTTCTCTCGATATTCAGTGCAGAATGAAAATCAGGTAGCCGAGTTCCAGGATGCGGGCATCGTATAATGGCTATTACCTCAGCCTTCCAAGCTGATGATGCGGGTTCGATTCCCGCTGCCCGCTCCAAGATGTGCTGATATAGCTCAGTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCT  >  NC_000913/4175587‑4176072

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: ‑