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

CAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGG  >  NC_000913/2870906‑2871389                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             cAGGGTAAGGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGGGCCGCGCATACTCCAGCCCCGGGCAGCCTTCCATTTGAAAGGCGCTCAGATGGGGCATTCCCTTGCCCGGCACGCGGGCGCATACCCACTGTTCCGTAACAGCAGGTTT‑GTCTTTTT‑GCCCAGCGAAGCGGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTCC‑ATCCGCCACAAAACGGGTATTGCAGCCGgttt                                                                                                                                                                                                                                           <  1:387005/251‑1 (MQ=255)  aGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGgtttt                                                                                                                                                                                                                                          >  2:411428/1‑251 (MQ=255)   gggTAAAGGCGTTACCTTTCAGCGCGACTAAACGCAGCTTGCGCTTGTGCCGCACATACTCCAGCACCTGGCAGCCTTCCAGTTGAAATGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCCCTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCCCAAAACGGGTATTGCAGCCGGTTTTg                                                                                                                                                                                                                                         >  2:313653/1‑251 (MQ=255)               tACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCCCGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAAcc                                                                                                                                                                                                                             >  2:389349/1‑251 (MQ=255)               tACCTTCCAGCGCGCTAAAACGCATCTTGCTCTTGTCCCGCGCAGACTCCAGCACCTGCCAGCCTTCCAGGTGAAAGGCGCTAAGATCGGGCTTTTCCTTCCCCGGCACGCGAGCGAATACCCTCTGTTCCGTAACAGTATGTTT‑GTTTTTTT‑GCCCAGCGAAGTTGACTTCACGGGCAGGAATTTTCAGGATTTTCGCTAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAAc                                                                                                                                                                                                                              <  1:427770/250‑1 (MQ=255)                      cAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAatat                                                                                                                                                                                                                      >  2:4139/1‑251 (MQ=255)                                            cgcTTGTGCCGCGCATACGCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGAATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAACACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACATTCACCATCAGGCTc                                                                                                                                                                                                >  1:298012/1‑251 (MQ=255)                                              cTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGCTTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCa                                                                                                                                                                                               >  1:236568/1‑250 (MQ=255)                                                 gggCCGGGCATACTCCGGCACCGGGCAGCCTCCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGGGAGCCCATAACCACGGTTCCTAAACAGCATGTTT‑GTCTTTTT‑GCCCGGCGAATCGGACTTCACGGGCATGAATTTTCAGGAATTGCGCTAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTGGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGc                                                                                                                                                                                           <  1:310519/249‑1 (MQ=255)                                                             ccccgcgcccgggCAGCCTCCCTGTTTAAAGGCGGTCAGACTGGCCATTTCCTTGCCGCGCGCGCGGGGGCAGAACCCATGTTCGGTACCA‑TATTTTT‑GTTTTTTTTCCCCGGCGAAGCTGACTACACGGGCATGAATTTTCAGGAATTCCGTCGTTGA‑ACCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTcca                                                                                                                                                                               <  2:86890/240‑1 (MQ=255)                                                               ccAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTccacc                                                                                                                                                                             >  1:379404/1‑251 (MQ=255)                                                                             cctcccAGTGGAAAGCCGCTCAGACCGGGCATTTCCTTGCCCGGCACGCGAGCTCATAACCACTGTCCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCAAAGCTGACTTCACGGGAATAAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTc                                                                                                                                                               <  2:290903/247‑1 (MQ=255)                                                                              cTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTcc                                                                                                                                                              >  1:61637/1‑251 (MQ=255)                                                                                     gtGAAAGGCGCCCAGATCGGGCATTTCCTTGCCCGGCCCGGGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCGGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTCCGCCAGGGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTgg                                                                                                                                                       <  1:4139/250‑1 (MQ=255)                                                                                      tGAAAGGCGCTCAGATCGGGCATTTCCTTGCCGGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTgg                                                                                                                                                       <  1:392665/250‑1 (MQ=255)                                                                                              gctgaGATCCGGCATTTCCTGTCCCGGCCCGCGAGCGCAAACCCACTGTTCCGACACGTCTTTTTT‑TTTTTTTT‑GCCCAGCGAATCTGCCTACAGGGGCAGAAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTcagca                                                                                                                                              <  2:184866/247‑1 (MQ=255)                                                                                                 cAGATCGGGCATTCCCTGGCCCGGCACGCGGGCGCATAACCGCTGGTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGATTTCACGGGATTGAATTTTCAGGAATTTCGCCAGTGG‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCaa                                                                                                                                                                                       <  1:241664/207‑1 (MQ=255)                                                                                                 cAGATCGGGAATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑CTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCaa                                                                                                                                                                                       >  2:241664/1‑207 (MQ=255)                                                                                                           attcctttGCCCGGCCCGGGAGCGCATAACCACTTTCCCGTAACAGCATGTTT‑TTATTTTT‑GCCCAGCGAAGCTGACTTCACGGTATTGAATTTTCAGGATTTTCGCAAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTg                                                                                                                                 <  2:347897/245‑1 (MQ=255)                                                                                                            tttCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGc                                                                                                                                >  1:140558/1‑251 (MQ=255)                                                                                                                                  gcATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGCGGTAAGTgg                                                                                                          >  2:110927/1‑250 (MQ=255)                                                                                                                                       ccctcttttCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCAAACCTGACTCCACGGGCATTATTTTTCAGTATTTTCGCAAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTa                                                                                                     <  2:393601/244‑1 (MQ=255)                                                                                                                                       aCCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGGGCCTTGCGGTTTACCGTGGAGGTAAGTGAGAtt                                                                                                      >  1:364497/1‑250 (MQ=255)                                                                                                                                            ggTTCCGTAACAG‑ATGTTTTGTCTTTTT‑GCCCAGCAAAGCGGACTTCACGGGCAGGAATTTCCAGGAATGTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGAGTATCaaa                                                                                                <  2:64965/250‑1 (MQ=255)                                                                                                                                                   aaaCAGCATGTTT‑GTCTTTTG‑GCCCAGCGAAGCTGCCTTCCCGGGAATGAATTTTCAGGAATTCCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATc                                                                                         <  1:245160/250‑1 (MQ=255)                                                                                                                                                    gaaacagagtttttttttttgtccccAGCTAAGCTGACCTCACGGGTATTATTTTTCATGAATTTCGCCATTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAAt                                                                                          <  2:351928/228‑1 (MQ=255)                                                                                                                                                      catcATGTTT‑TTCTTTTT‑CCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTCCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAtt                                                                                      <  1:216519/248‑1 (MQ=255)                                                                                                                                                      cAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGTCTCCACGGGCATGATTTTTCGGGAGTTTCGCAAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCAtt                                                                                      <  2:119863/251‑1 (MQ=255)                                                                                                                                                          aTGTTT‑GTCTTTTT‑GCCCAGCGCAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTttgt                                                                                  >  1:181432/1‑251 (MQ=255)                                                                                                                                                                                  aaccgtactaaaCGGGCATGAATTTTCGGGATGTTC‑CCAGGGCCAACCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCTCCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCAcc                                                            <  2:247/240‑1 (MQ=255)                                                                                                                                                                                        aCTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTca                                                      <  1:111027/251‑1 (MQ=255)                                                                                                                                                                                                  cATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATcc                                            >  1:189487/1‑251 (MQ=255)                                                                                                                                                                                                         tttCAGGAATTCCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTccc                                      <  1:111539/250‑1 (MQ=255)                                                                                                                                                                                                              ggAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTc                                >  2:94727/1‑251 (MQ=255)                                                                                                                                                                                                               gattttCGCCAGGGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTcc                               <  1:336330/248‑1 (MQ=255)                                                                                                                                                                                                               gAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTcc                               >  2:399031/1‑251 (MQ=255)                                                                                                                                                                                                                ccctttccgctttacgtCCGACACAAAACGGTTATTGCAGCCGGTTTTTAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTcc                               <  2:260107/235‑1 (MQ=255)                                                                                                                                                                                                                 aTTTCGCCAGTGC‑CTCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCCCCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGATGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCATTCATTTTGTTGCCTTAATGTGCAGCGCAACCGCTTCACAGGCAATCCCTTCCCCACGTCCgg                             >  2:60363/1‑251 (MQ=255)                                                                                                                                                                                                                            gC‑ATCCGCCACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTcc                               >  1:260107/1‑238 (MQ=255)                                                                                                                                                                                                                                    ccACAAAACGGGTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCg           <  1:422744/251‑1 (MQ=255)                                                                                                                                                                                                                                             gggTATTGCAGCCGGTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTgg  >  1:93298/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             CAGGGTAAAGGCGTTACCTTTCAGCGCGCCTAAACGCAGCTTGCGCTTGTGCCGCGCATACTCCAGCACCTGGCAGCCTTCCAGTTGAAAGGCGCTCAGATCGGGCATTTCCTTGCCCGGCACGCGAGCGCATAACCACTGTTCCGTAACAGCATGTTT‑GTCTTTTT‑GCCCAGCGAAGCTGACTTCACGGGCATGAATTTTCAGGAATTTCGCCAGTGC‑ATCCGCCACAAAACGGGTATTGCAGCCGTTTTTGAGGATTCTAACCAGAATATGCTCACCTTCACCATCAGGCTCAAAGCCCAAATCTTCCACCACCACAAAGTCTTCCGGATTGGCTTTCAGCAGCCCGGTGCCTTGCGGTTTACCGTGGAGGTAAGTGAGATTATCAAACTCAATCATTTTGTTGCCTTAATGAGTAGCGCCACCGCTTCACAGGCAATCCCTTCCCCACGTCCGGTAAATCCCAGTTTTTCCGTAGTAGTGG  >  NC_000913/2870906‑2871389

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