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

GGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACG  >  NC_000913/2870915‑2871397                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     ggCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTc                                                                                                                                                                                                                                          >  2:442227/1‑251 (MQ=255)  gCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGGTTCt                                                                                                                                                                                                                                         >  2:347430/1‑251 (MQ=255)   cGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCt                                                                                                                                                                                                                                         >  2:68915/1‑250 (MQ=255)        ccTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAAcc                                                                                                                                                                                                                                     >  2:559688/1‑249 (MQ=255)                 gcgcCAAAACGCAGCTTGCTCTTGCCCCGCGCATACCCCGGCCCCTGGCAGCCTTCCAGTTGAAAGGCGCCCAGAGCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATACCCCCTGTTCCGTAACAGCATTTTTGTCTTTTTGCCCAGCAAAGCTGAC‑TCCACGGGCAGGAATTTTCAGGAATTTCGCAAGTGCTTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTc                                                                                                                                                                                                                          <  2:149298/251‑1 (MQ=255)                 gcccaaaaacgcgtcTTGCGTTTGGGCCGCGCATACTCCAGCACCGGGCACCCTGCCAGTGGAAAGGCGCTAATGTCGGTCTTTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAAAAGCATGTTTGTTTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCTGGAATTTTCAGGAATTCCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTc                                                                                                                                                                                                                          <  1:314381/237‑1 (MQ=255)                    ccTAAACGAGGCTGGCGCTTGTGCCGCGCATACCCCACCCCCTGGCAGCCTTCCAGTTGAAAGGCCCTCAGTTCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATACCCCCTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGAGAAGCTGCC‑TCCACGGGCATGAATTTCCAGAATTTTCGCCGGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCAcc                                                                                                                                                                                                                       <  2:9747/251‑1 (MQ=255)                          cGCAGCTTGCGCTTGTGCCGCGCATACTCCCGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCAcc                                                                                                                                                                                                                 >  2:366715/1‑251 (MQ=255)                             aGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATc                                                                                                                                                                                                              >  1:175435/1‑251 (MQ=255)                                      ttgggCCGCGCATTCTCCAGCCCCGGGAGGCTTTCCAGTTAAAAGGCCCTAAGATCGGGCATTTCCTTGCCCGGC‑CGCCGAGCGAATACCCACTGTTCCGTAACAGAATGTTTGTTTTTTTGCCCAGCGAAGC‑GACCTTCACGGGCATGATTTTTCGGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATGGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCaa                                                                                                                                                                                                      <  2:201653/246‑1 (MQ=255)                                       tgtgCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACTTTCACCATCAGGCTCaa                                                                                                                                                                                                      >  1:68500/1‑249 (MQ=255)                                        gtcccGCGCATCCTCCACCCCCGGGCACCCTCCCCGTTGAAAGGCGCCCAGATCGGGCTTTTCCGTGCCCGGCGCGC‑GAGCCCATAACCCCTGTTCCGTACCAGCATGTTTGTCTTTTTGCCCAGCAAAGCTGAT‑TTCCCGGCCATGAATTTCCAGGAATTTCCCCTGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAg                                                                                                                                                                                                    <  1:6575/247‑1 (MQ=255)                                                         gCACCTGGCAGCCTTCCAGTTGAACGGAGCTCACATCGGGCATTTCCTTGCCAGGCACGC‑GAGCGCATCACCCCTGTTCCGTAACAGCATGTTTGTCTTTTTGCCAAGCGACGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCAGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTcccccacc                                                                                                                                                                                  >  2:297818/1‑251 (MQ=255)                                                               ggCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGc                                                                                                                                                                                                                            >  1:231048/1‑203 (MQ=255)                                                                                     gcTCAGATCGGGCATTTCCTTGCCCGACACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGACAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCACAAATCTTCCACCACCACAAAGTCTTCCGAACTGACTTTcagc                                                                                                                                                       >  2:25287/1‑250 (MQ=255)                                                                                           aTCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTATGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCgg                                                                                                                                                >  1:87630/1‑251 (MQ=255)                                                                                                         ggCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGtt                                                                                                                                    <  1:559688/248‑1 (MQ=255)                                                                                                          gCCCGGCACGC‑GAGCGCATAACCACGTTTCCGTAAAAGCTTGTTTGTCTTTTTGCCCAGCGAAGCTGCC‑TTCACGGGCAGGAATTTGCAGGAATTCCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTAc                                                                                                                                 <  2:470540/251‑1 (MQ=255)                                                                                                             cGGCACGG‑GAGCGCATAACCACTGTTCAGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGAGCATGAATTTTCAGGCATTTCGCCAGTGCATCCGCCACAACACGGGTATTGAGGCCGGTTTTGAGGATTCTAACCAg                                                                                                                                                                                                                                   >  1:486387/1‑150 (MQ=255)                                                                                                               gCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGGATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTgg                                                                                                                            >  2:12153/1‑251 (MQ=255)                                                                                                                          gcATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCATAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTg                                                                                                                   >  1:588450/1‑250 (MQ=255)                                                                                                                                ccACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTAt                                                                                                            >  2:555014/1‑251 (MQ=255)                                                                                                                                    ggTTCCGTAACAGCTTGTTTTTTTTTTTGCCCAGCGAAGCTGCC‑TTCCCGGTCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCaaa                                                                                                        <  1:196313/250‑1 (MQ=255)                                                                                                                                          gTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAAt                                                                                                  >  1:72350/1‑251 (MQ=255)                                                                                                                                                  gtGTTTGTCTTGTCGCCCAGCAAAGCTGCC‑TTCAGGGGCATGAATTGTCAGGAATTTCGCCGGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGTATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTttgt                                                                                          <  1:479544/250‑1 (MQ=255)                                                                                                                                                    gttttgcTTTTTGCCCAGCGAAGCGGAT‑GCCACGGTAATGAATTTGCAGAATTTTCGCAAGGGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTGTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATgttgtt                                                                                         <  2:588450/244‑1 (MQ=255)                                                                                                                                                           tttttGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCCGTGCATCCGCCACAAAACGGGTCTTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGGGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAAt                                                                                 >  1:284496/1‑251 (MQ=255)                                                                                                                                                                       gaaccTGAT‑TTCACGGGCAGTATTTTTCAGTAATGTCGCCAGTGCAGCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTAGCAAACGCAATCATTTTGTTGCCTCAATGAGTAGCTCCAc                                                                     <  2:208291/247‑1 (MQ=255)                                                                                                                                                                       gAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCAc                                                                     >  1:64154/1‑251 (MQ=255)                                                                                                                                                                              cgctttacgggagtaaaTTTCCAGGAATTTCTCCGGTCCCCCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAACATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTAGCAAACTCAATCATTTTGTTGCCTTAATTAGGAGCGCCACCGCTTc                                                               <  2:80468/236‑1 (MQ=255)                                                                                                                                                                                          aTGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATccc                                                   >  1:205409/1‑251 (MQ=255)                                                                                                                                                                                               ttttcaagattctcacctggtcacccGCCACAAAACGGGTATTGCAGCCTGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTGGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTcc                                               <  2:537168/226‑1 (MQ=255)                                                                                                                                                                                                     ggAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTc                                        <  1:95991/251‑1 (MQ=255)                                                                                                                                                                                                     gaaactccacccgTGCAACCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAAGCTTCCACCACCACAAAGTCGTCCGGATTGGCTTTCATCAGCCCGGTGCCTTGCGGTTTACCGTTGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTGGCGCCAACGCTTCACGGGCAATTGCTTCCCCACGTc                                        <  2:27629/239‑1 (MQ=255)                                                                                                                                                                                                                         gCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTcc                    >  1:80173/1‑251 (MQ=255)                                                                                                                                                                                                                             caAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCgtag                <  1:221368/251‑1 (MQ=255)                                                                                                                                                                                                                               aaaCGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCCTTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTACATCCCAGTTTTTCCgtagt               >  2:1136/1‑250 (MQ=255)                                                                                                                                                                                                                                aaCGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCgtagt               >  1:407657/1‑249 (MQ=255)                                                                                                                                                                                                                                 aCGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCTGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGGGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCACCACGTCCGGTAAATCCTAGTTTTTCCGTagtagt            <  1:25287/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGc         >  1:269823/1‑251 (MQ=255)                                                                                                                                                                                                                                         tGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCCTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTGTACCGTGGAGGTAAGTGAGATGATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCTCCACCGCTTCACAGGCAAGCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTc     <  1:566526/250‑1 (MQ=255)                                                                                                                                                                                                                                           cAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACg  <  1:68915/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     GGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGC‑GAGCGCATAACCACTGTTCCGTAACAGCATGTTTGTCTTTTTGCCCAGCGAAGCTGAC‑TTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGCATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGGCTTTCACG  >  NC_000913/2870915‑2871397

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