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	100.0%	245.2 / NA	64	N51T (AAC→ACC)	truD	tRNA pseudouridine(13) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (33/31);  total (33/31)

ACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTT  >  NC_000913/2870921‑2871392                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                aCCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAAcc                                                                                                                                                                                                                                >  1:276096/1‑250 (MQ=255)     ttCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAat                                                                                                                                                                                                                           >  1:678092/1‑251 (MQ=255)      ccAGCGCTCCTAAACGCGGCTGGGCCTTGTCCCGCGCTTACTCCACCTCCGGGCAGCCTTCCAGTGAAAAGGCGCTCAGATCGGGCTTTTCTTGGCCCGTCACGCGAGCGAATAACCTCTGTT‑CC‑TTAACAGCATGTTGGTCTGTTGCCCCAGGAAAGTTAACTGCACGGGTATGATTTG‑TCAGGATTTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAat                                                                                                                                                                                                                           <  2:483893/249‑1 (MQ=255)       cAGCGCGCCTAAACGCAGCGTGCGCTGGTGCCGCGCATACCCCAGCCCCTGGCAGCCTTCCAGTGGAAAGGCGCTCAGAGCGGTCATTTCCTGGCCCGGCACGCGGGCGCATAACCAC‑GTT‑CCCGTAACAGCATTTTTGTTTTTTTGCCCCGCGAAGCCGACTTCACGGGAATGAATTT‑TCAGGATTTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAatat                                                                                                                                                                                                                         <  2:88876/251‑1 (MQ=255)          cgcgcCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCt                                                                                                                                                                                                                      >  1:400572/1‑251 (MQ=255)                         cTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATc                                                                                                                                                                                                         >  2:403917/1‑249 (MQ=255)                            tctcTTGTGCCCCGCATACTCCAGCACCGGCCTGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTCCCTTGCCCGGCACGCGAGCGCATAACCACGGTC‑CC‑GTAACCGCTTTTTTGTTTTTTCGCCCAGAGAAGCTGATTTACCGGGCATGAATTT‑TCAGAAATTTCGCCGGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGc                                                                                                                                                                                                     <  1:400407/247‑1 (MQ=255)                              gcTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTc                                                                                                                                                                                                   >  2:129530/1‑250 (MQ=255)                                 tgtgCCGCGCATCCTCCAGCCCCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGAATTTCCTTGCCGGGCACGCGAGCGCATAACCCCTGTT‑CC‑GTAACAGCATTTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCaa                                                                                                                                                                                                 <  1:596978/249‑1 (MQ=255)                                   tgCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGcc                                                                                                                                                                                             >  1:690835/1‑251 (MQ=255)                                   tgCCGAGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGcc                                                                                                                                                                                             >  1:122787/1‑251 (MQ=255)                                             aCTCCAACACATGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTcc                                                                                                                                                                                   >  2:100561/1‑251 (MQ=255)                                                          gCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGt                                                                                                                                                                      >  1:717900/1‑251 (MQ=255)                                                           cAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCAt                                                                                                                                                                                                          >  1:153664/1‑214 (MQ=255)                                                           cAGCCTTCCAGTTGAAAGGCGCTCAGAGCGGCAATTGCCTGCCCCGGCACGGGGCCTCATAACC‑CTGTTTCC‑CTAACAGAATTTTTTTTTTTTTGCCCGGCGAAGCGGACTTCACGGGGATAATTTT‑TGAGGAATTTCGCCATTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCAt                                                                                                                                                                                                          <  2:153664/214‑1 (MQ=255)                                                                  ccAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATACCCACTTTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGa                                                                                                                                                              <  1:285598/251‑1 (MQ=255)                                                                         aaaGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCCAATCTTCCACCACCACAAAGTCTGCCGGATTGGCtt                                                                                                                                                       >  2:473069/1‑251 (MQ=255)                                                                             gcgcTCAGACCGGGCATTTCCTTGCCCGGCACGGGAGCGCATAACCACTGTC‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCAAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTcag                                                                                                                                                   <  1:401537/251‑1 (MQ=255)                                                                              cccTCAGACCGGGCATTCCCTTGCCCGCAACGGGGGCGCAAAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTCCACGGGCATGAATTT‑CCAGAAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTcagc                                                                                                                                                  <  1:643566/249‑1 (MQ=255)                                                                               gctcgaatggggCATTCCCTTGCCGGCCACGCGTGCGCAAAACCATTGTT‑CC‑GTAACAGAATTTTTTTCTTTTTCCCCAGCGAAGCTTACTTCACGGGCATGATTTT‑CCAGGAATTCCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTcagca                                                                                                                                                 <  2:62447/242‑1 (MQ=255)                                                                                  cAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGcc                                                                                                                                              >  2:486605/1‑251 (MQ=255)                                                                                  cAGAGCGGGCATTTCCTTGCCCGGCCCGCGAGCGCATACCCA‑TGTT‑CCCTTAACGGCATGTTTGTCTTGTGGCCCAGCGAGGCTGACTTTTCGGGCGTGAATTT‑TCAGGATTTTCGCTA‑TGCCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGAGTGGCTTTCAGCAGcc                                                                                                                                              <  2:331649/251‑1 (MQ=255)                                                                                   aGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGccc                                                                                                                                             >  2:336211/1‑251 (MQ=255)                                                                                    gATCGGGCATTTACTTGCCCGGCACGCGAGCGAATAACCACTGTT‑CC‑GTAACAGCATGTTTGTATTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCAACACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCg                                                                                                                                            >  2:439733/1‑251 (MQ=255)                                                                                      tCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCggg                                                                                                                                          >  1:576861/1‑250 (MQ=255)                                                                                              ttCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCACAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCg                                                                                                                                  >  2:444159/1‑251 (MQ=255)                                                                                                   tccccGGCCCGCGGGCGCATAACCCCTGGT‑CC‑GAAACAGCTTGTTTGTCTTTTTGCCCGGCGAAGCTGTCTCCGCGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGttt                                                                                                                              <  2:594569/248‑1 (MQ=255)                                                                                                   tGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGttt                                                                                                                              >  1:171524/1‑250 (MQ=255)                                                                                                       cgccccgctagcatttAACCTCTGTT‑CC‑GAAACAGCTTGTTTTTCTTTTTTCCAAGCGAAGCTTACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTCC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGt                                                                                                                         <  1:528017/236‑1 (MQ=255)                                                                                                        ggcccgcgtgcgctaaacccttGTT‑CG‑GTAACAGCTTGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATTATTTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTg                                                                                                                        <  1:109461/230‑1 (MQ=255)                                                                                                          cACGCGAGCGCATAACCCCTGTT‑CC‑GTAACAGCATTTTTGTCTTTTGCCCCAGCGAAGCTGCCTTCAGGGGCATGATTTT‑TCAGGAATTTCGCCAGTTC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGa                                                                                                                      <  2:678092/251‑1 (MQ=255)                                                                                                              cgAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGt                                                                                                                   >  2:377381/1‑250 (MQ=255)                                                                                                                  cgcATACCCCCTTTT‑CC‑GTATCAACATGTTGTTCTTTTTGCCCAGCGAGGTTGACTTCCCGGCCATGAATTT‑TCAGGATTTCCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTg                                                                                                              <  1:544873/251‑1 (MQ=255)                                                                                                                     aTAACCCCTTTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATAAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTgaga                                                                                                           <  1:403917/251‑1 (MQ=255)                                                                                                                             tGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGCGGTAAGTGAGATTATCaa                                                                                                    >  2:571932/1‑250 (MQ=255)                                                                                                                                 ccc‑GTAACAACATGTTTGTTTTTTTCCCCAGCGAAGCTGACTTCACGGGCAAGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTc                                                                                                <  2:171524/250‑1 (MQ=255)                                                                                                                                     gTAACAGCTTGTTTGTTTTTTTCCCCAGGGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑AGCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAAt                                                                                             <  2:719361/251‑1 (MQ=255)                                                                                                                                        acggCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAt                                                                                          <  1:691549/248‑1 (MQ=255)                                                                                                                                           gCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAtttt                                                                                       >  1:54543/1‑251 (MQ=255)                                                                                                                                               gTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAACGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTttgttg                                                                                   >  1:97838/1‑251 (MQ=255)                                                                                                                                                tttGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGc                                                                                  >  1:416729/1‑251 (MQ=255)                                                                                                                                                  tGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCt                                                                                >  1:273284/1‑251 (MQ=255)                                                                                                                                                      tttttGCCCAGCAAAGCTTACTTCACGGGAATGAATTT‑ACAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTaa                                                                             <  1:448730/250‑1 (MQ=255)                                                                                                                                                       tttggCCCAGCAAAGCTAACTTCACGGGCATGAATTT‑TCAGGAATTCCGCCAGGGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTa                                                                              <  1:123940/244‑1 (MQ=255)                                                                                                                                                       ttgtGCCCAGCAAAGCTGACTCCACGGGCATGAATTT‑ACAGAAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAAt                                                                            <  1:641947/247‑1 (MQ=255)                                                                                                                                                              cAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAgcg                                                                    <  1:204240/251‑1 (MQ=255)                                                                                                                                                                         ccTTCACGGACAT‑AATTTGTAAGGAATTGCGCCAGTGC‑AGCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTCCAGCAGCCCGGTTCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTca                                                         <  1:453558/250‑1 (MQ=255)                                                                                                                                                                                 ggCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCCTTTTGTTGCCTTAATGTGTAGCGCCACCGCTTCACAGGCAAt                                                 >  2:150311/1‑251 (MQ=255)                                                                                                                                                                                      gAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCtt                                            <  1:550194/251‑1 (MQ=255)                                                                                                                                                                                        aTTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTcc                                          >  2:98325/1‑251 (MQ=255)                                                                                                                                                                                              cAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCAc                                      <  1:100561/250‑1 (MQ=255)                                                                                                                                                                                              cAGGAATTTCGCCAGGGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCAc                                      <  1:95345/250‑1 (MQ=255)                                                                                                                                                                                               aGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGt                                    >  2:66824/1‑251 (MQ=255)                                                                                                                                                                                                            gTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGCTTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATccc                        >  1:655824/1‑250 (MQ=255)                                                                                                                                                                                                            gTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATcccc                       >  2:165189/1‑250 (MQ=255)                                                                                                                                                                                                             tGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGTGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATccc                        >  2:511892/1‑249 (MQ=255)                                                                                                                                                                                                                 aTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGttt                   >  1:29332/1‑251 (MQ=255)                                                                                                                                                                                                                  cccGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCGTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAAGCCCAGtttt                  <  2:654490/250‑1 (MQ=255)                                                                                                                                                                                                                    agCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTc                >  2:402041/2‑251 (MQ=255)                                                                                                                                                                                                                     gCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGTCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTCTCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAGTCCCAGTTTTTcc               >  1:682115/1‑251 (MQ=255)                                                                                                                                                                                                                           aaaCGGGTATTGCAGCCGGTTTGGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCgtagta         <  2:266935/251‑1 (MQ=255)                                                                                                                                                                                                                                 gTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCGTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCt   <  2:122787/251‑1 (MQ=255)                                                                                                                                                                                                                                 gTATTGCAGCCGGTGTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCAGCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTc                                                                                                <  1:479116/158‑1 (MQ=255)                                                                                                                                                                                                                                  tATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATACCAGTTTTTCCGTAGTAGTGGCtt  >  1:33066/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                ACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTT‑CC‑GTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGACTTCACGGGCATGAATTT‑TCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTT  >  NC_000913/2870921‑2871392

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