Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,871,153	T→G	N51T (AAC→ACC)	truD ←	tRNA pseudouridine(13) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,871,153	0	T	G	97.4%	136.0 / ‑3.7	38	N51T (AAC→ACC)	truD	tRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (20/17);  total (21/17)
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.33e-01

AGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTT  >  NC_000913/2870907‑2871399                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aGGGTAAAGGCGTTAACTTTCAGAGCGCCTAAACGAAGCTTGCGCTTGTGAAGAGCATACTCCAGCAACTGGCAGCCTCCCAGTTGAAAGGCGCTCAGATCGGGCATTTACTTGACCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTAGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGgtttt                                                                                                                                                                                                                                                    >  1:416299/1‑251 (MQ=255)         ggCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTc                                                                                                                                                                                                                                            >  1:83316/1‑251 (MQ=255)                                 aCGCAGCTTGCGCTTAAACCGCGCAAACTCCAGCACCTGGCTGCCTTCCAGTTGAAAGGCGCTCAGCTAGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCTGCGAAGCTGACTTCACGGCCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGCATATGCTCACCTTCAc                                                                                                                                                                                                                    >  2:90989/1‑251 (MQ=255)                                     aGCTTGCGCTTGTGCCGCGCATACCCCAGCCCCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTCCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTTTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGATTTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATc                                                                                                                                                                                                                <  1:294710/251‑1 (MQ=255)                                         tGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCCGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGc                                                                                                                                                                                                            >  2:307744/1‑251 (MQ=255)                                               ggggCCGCGCACACCCCAGCACCGGCCAGCCTCCCAGTTGAAAGGCGCTCATATGGGGCATTTCCTTGCCCGGCACGCGAGCTCATAACCCCTGTTCCGTAACTGCATGTTTGTCTTTGT‑GCCCAGCGAAGCTGATTTCTGGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGATGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAg                                                                                                                                                                                                      <  1:447769/248‑1 (MQ=255)                                                gtgCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGc                                                                                                                                                                                                     >  1:47012/1‑251 (MQ=255)                                                    cgcgCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCa                                                                                                                                                                                                  >  1:384845/1‑250 (MQ=255)                                                     gcgcATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCaa                                                                                                                                                                                                 >  1:106968/1‑250 (MQ=255)                                                               cAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTccaccc                                                                                                                                                                                      >  2:68684/1‑250 (MQ=255)                                                                           gCCTTCCAGTTGAAAGGCGCTCAGATGGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCGGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCtt                                                                                                                                                                          <  1:201895/251‑1 (MQ=255)                                                                                        aaGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCttt                                                                                                                                                             >  1:16731/1‑251 (MQ=255)                                                                                          ggCGCTCAGAGCGGGCCTTCCTTTGCCCGCCACGCGGGCGCATAACCTCCGGTCTGTAACAGCTTGTTTGTATTTTA‑GCCCAGCGAATCGCATTACCCGGGCATGAATTTACAGGAATTCCTCCGGCCCCACCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCCCAAAGCCCAAATCTCCCACCACCACAAAGTCTTCCGGATTGGCTTTca                                                                                                                                                           <  1:385950/251‑1 (MQ=255)                                                                                             gcTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTcagc                                                                                                                                                         >  1:23043/1‑250 (MQ=255)                                                                                                 aGATCGGGCATTTCCATGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGCCTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACACAGTCTTCCGGATTGGCTTTCAGCAGccc                                                                                                                                                    >  2:42166/1‑251 (MQ=255)                                                                                                              ccgtccccGGCACGCGAGCCCAAAACCTCTGTCCCGTAACAGCATGTTTTTCTTTTT‑CCCCAGCGAAGCTGACTTCAGGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGt                                                                                                                                       <  1:160239/246‑1 (MQ=255)                                                                                                                 ggCCCGGCACGCGAGCAAATAACCTCTGTCCCGTAACAGCATGTTTGTTTTTTT‑GCCCAGCGAAGCTGACTCCAGGGGCATAATTTTTCAGAAATTTCCCAAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGttt                                                                                                                                     <  1:454118/249‑1 (MQ=255)                                                                                                                            cgAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAACTGACTTCACGGGCATGAATTTTCAGGAATTTCGCAAGTGCATCCGCCACACAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTgtcggt                                                                                                                          >  1:417851/1‑245 (MQ=255)                                                                                                                                            gtgccGTAACAGCATGTTTGTTTTTTC‑GCCCAGCTAATCTGACGGCTCGGGCATGATTTGTCAGGAATTTCGCCTGTGCATCCGCCACAAAACGTGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGACTTCCGGATTGGCTTTCAGGATCCCGGTTCCTTGCGTTTTACCGTGGAGGTAAGTGAGATTATCaaa                                                                                                          <  1:409449/247‑1 (MQ=255)                                                                                                                                                        aaTGTTT‑TCTTTTTCGCCCAGCGAAGCTCTCTTCCGGGGCTGGAATTTTCAGGATTGTCGCCAGTGCACCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAtttt                                                                                              <  2:344640/249‑1 (MQ=255)                                                                                                                                                            tttGTCTTTTT‑GCCCAGCAAAGCTGACTTCACGGGCATGATTTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTttgttg                                                                                          <  1:304022/250‑1 (MQ=255)                                                                                                                                                              tGTCTTTTT‑GCCCAGCGAAGCTCACTTCACGGGCATGAATTTACAGGATTTGCGCCAGAACATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGcc                                                                                        <  1:10629/250‑1 (MQ=255)                                                                                                                                                                ttttttg‑GCCCTGCGACGCTGACTTCAGGGGCATGAATTTTCAGGAATTTCGCCTGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAATCTCAATCATTTTGTTGCCtt                                                                                      <  1:94024/248‑1 (MQ=255)                                                                                                                                                                    ttt‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATg                                                                                  >  1:192847/1‑250 (MQ=255)                                                                                                                                                                          ccAGCGAAGCTGACTTAACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATCTGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAgc                                                                            >  2:79647/1‑251 (MQ=255)                                                                                                                                                                             gcggagctgtctccaggggaagtatttttAAGGATTTTCGTCCGTGCAGCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGcc                                                                         <  2:106968/226‑1 (MQ=255)                                                                                                                                                                              gaaaGCTGACTTCACGGGCAGGAATTTTCAGGAATTCCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCa                                                                        <  1:68684/249‑1 (MQ=255)                                                                                                                                                                               gAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCAc                                                                       >  1:158272/1‑251 (MQ=255)                                                                                                                                                                                                        ttCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTcccc                                               >  2:407250/1‑250 (MQ=255)                                                                                                                                                                                                        ttCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTAAATCATTTTGTTGCCTTAATGAGTAGCGCCAACGCTTCACAGGCAATCCCTTcccc                                               >  2:147258/1‑250 (MQ=255)                                                                                                                                                                                                         tCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCa                                                                                                      >  1:88011/1‑194 (MQ=255)                                                                                                                                                                                                                      cAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGAAAATcc                                >  2:392230/1‑251 (MQ=255)                                                                                                                                                                                                                       aGTGCATCCACCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATccc                               <  2:92478/251‑1 (MQ=255)                                                                                                                                                                                                                            aTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGttt                          >  2:351762/1‑251 (MQ=255)                                                                                                                                                                                                                            aTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGttt                          >  1:168135/1‑251 (MQ=255)                                                                                                                                                                                                                              ccGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGtttt                         <  1:351762/250‑1 (MQ=255)                                                                                                                                                                                                                                    caAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCgt                    >  1:151532/1‑249 (MQ=255)                                                                                                                                                                                                                                                    gCCGGTTTTTAGGATTCTAACCAGAATATGCTCACCTTCACCAACAGGCTCAAAGCCCAAATCTTCCCCCCCCACAAAGTCTTCCGGATTGGCTTTCAGCAGCTCGGTGCCTGGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGTAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGtt  <  1:234414/251‑1 (MQ=255)                                                                                                                                                                                                                                                    gCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGtt  <  1:277409/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACGTT  >  NC_000913/2870907‑2871399

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