Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	694,941	G→A	L5L (CTC→CTT)	miaB ←	isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	694,941	0	G	A	100.0%	228.1 / NA	61	L5L (CTC→CTT)	miaB	isopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/31);  total (30/31)

TTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGC  >  NC_000913/694695‑695174                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             ttGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGtt                                                                                                                                                                                                                                          >  1:584547/1‑251 (MQ=255)  tGCGATGCCAAGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGt                                                                                                                                                                                                                                           <  2:260661/249‑1 (MQ=255)          cACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCGTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCa                                                                                                                                                                                                                                 <  2:507090/251‑1 (MQ=255)                gCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACt                                                                                                                                                                                                                           >  1:74323/1‑251 (MQ=255)                       aCGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCtt                                                                                                                                                                                                                    <  1:403717/251‑1 (MQ=255)                           cGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGa                                                                                                                                                                                                                <  1:259936/251‑1 (MQ=255)                                       cTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCAGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGcc                                                                                                                                                                                                      >  2:254179/1‑249 (MQ=255)                                         ggATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAgg                                                                                                                                                                                                   <  2:326301/250‑1 (MQ=255)                                          gATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTTGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGaa                                                                                                                                                                                                 <  2:297557/251‑1 (MQ=255)                                          gATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTACTGGAACACTTTTTCCTGCGCCTTCTCGCGGCTTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCCGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGa                                                                                                                                                                                                  >  1:205209/1‑250 (MQ=255)                                             tcttctCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGc                                                                                                                                                                                              >  2:469252/1‑251 (MQ=255)                                                   cttgtagcAGTTCCCAGCGACCCAACTGATGGAAGACTTTTTCCTGGGCCTTCTCGCGGTGTGAGCAGGTGTTCAGCAGCAGCACATCCGGTTCTTCCGCCACGTCGGGCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAGCTTCTATGAAGCGTGCGCGGTCAACTGACAGCCCGAGGTTTTAAGAAGAAGTGTGTTGGTTATGGACTTGCTCTGGCGAGTTAGGAGCCAgttagtcaggg                                                                                                                                                                                          <  2:263380/242‑11 (MQ=255)                                                   ctTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAgg                                                                                                                                                                                           <  1:440575/248‑1 (MQ=255)                                                   cgtttAACAGGTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTGCTGGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTGCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGgct                                                                                                                                                                                        <  1:513831/249‑2 (MQ=255)                                                     tttAACAGTTTCCAGCGACCCAACGGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCAGCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGgcg                                                                                                                                                                                        <  1:513839/249‑1 (MQ=255)                                                              ttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCCGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAAt                                                                                                                                                                            >  2:198268/1‑251 (MQ=255)                                                              ttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTaa                                                                                                                                                                             <  1:149204/250‑1 (MQ=255)                                                              ttCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAAt                                                                                                                                                                            >  1:424718/1‑251 (MQ=255)                                                                               aTGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAAtagtag                                                                                                                                                                                                        >  2:173313/1‑207 (MQ=255)                                                                               aTGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGAAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAAtagtag                                                                                                                                                                                                        <  1:173313/207‑1 (MQ=255)                                                                                   aaGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGAcc                                                                                                                                                        >  2:556249/1‑250 (MQ=255)                                                                                    aGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGAcc                                                                                                                                                        >  2:471521/1‑249 (MQ=255)                                                                                            tCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGc                                                                                                                                               <  1:211136/250‑1 (MQ=255)                                                                                              cTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑tt                                                                                                                                           >  1:33356/1‑251 (MQ=255)                                                                                                   ccTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATc                                                                                                                                        >  1:507084/1‑249 (MQ=255)                                                                                                                  aGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTaaa                                                                                                                        >  1:36322/1‑250 (MQ=255)                                                                                                                  aGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTaaa                                                                                                                        >  1:537833/1‑250 (MQ=255)                                                                                                                      ggTGTTCAGCAGCAGCACAGCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAGCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGGGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATcc                                                                                                                   <  1:531350/251‑1 (MQ=255)                                                                                                                         gtTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGt                                                                                                                >  1:306322/1‑251 (MQ=255)                                                                                                                            cagcagcagCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTaa                                                                                                              <  2:238492/250‑1 (MQ=255)                                                                                                                               cagcagcaCATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACtt                                                                                                          >  2:457008/1‑251 (MQ=255)                                                                                                                                   agcaCATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑g                                                                                                     <  2:158612/251‑1 (MQ=255)                                                                                                                                              ttcttcCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTa                                                                                           >  2:219820/1‑250 (MQ=255)                                                                                                                                                 ttcCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAAc                                                                                         >  1:247360/1‑249 (MQ=255)                                                                                                                                                  tcCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACa                                                                                        <  1:295791/249‑1 (MQ=255)                                                                                                                                                    cGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATTTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGgat                                                                                    <  1:544755/251‑1 (MQ=255)                                                                                                                                                       cACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGgatgat                                                                                 <  1:36222/251‑1 (MQ=255)                                                                                                                                                          gTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGt                                                                                                                <  2:479890/218‑1 (MQ=255)                                                                                                                                                          gTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGt                                                                                                                >  1:479890/1‑218 (MQ=255)                                                                                                                                                                    tGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACaa                                                                     <  2:32215/250‑1 (MQ=255)                                                                                                                                                                     gATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACaaa                                                                    <  1:556249/250‑1 (MQ=255)                                                                                                                                                                     gATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAAt                                                                   >  1:219248/1‑251 (MQ=255)                                                                                                                                                                      aTAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATc                                                                  <  2:538067/251‑1 (MQ=255)                                                                                                                                                                           cGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATcaacc                                                              >  2:12915/1‑250 (MQ=255)                                                                                                                                                                           cGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAAT‑CTGGGGCCATTAGGGTAAT‑CTTTTCCGGCGGTGGGGCCAGGATGGGCGGTTATCACCCCTTGGGGGAAAAATTCCT‑TAAAATTAAAACCATTTTGCAAACAGGTTGATTAACCAT‑ACAAATcaacca                                                             >  2:404598/1‑250 (MQ=255)                                                                                                                                                                              gggtggCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACg                                                          <  1:183412/251‑1 (MQ=255)                                                                                                                                                                                  ggCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGaa                                                       >  1:405753/1‑250 (MQ=255)                                                                                                                                                                                            agATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGACCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGACAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGt                                            >  1:435558/1‑251 (MQ=255)                                                                                                                                                                                             gATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGTCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCAGATTGCTAACAGGATTATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTc                                           <  1:285121/251‑1 (MQ=255)                                                                                                                                                                                              aTCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTcg                                          >  2:471689/1‑251 (MQ=255)                                                                                                                                                                                              aTCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTcg                                          <  1:7846/251‑1 (MQ=255)                                                                                                                                                                                               acGGCCATCTTCGCTGAATCGTACTAGTTCATATGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGAAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTCTCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACACACGGAAATTGCCATTGTcgg                                         >  2:549989/2‑251 (MQ=255)                                                                                                                                                                                                 ggCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTcggcg                                       >  1:501537/1‑251 (MQ=255)                                                                                                                                                                                                   ccATCTTCGATGAATCGTACTCGTTCATCTGACAGCTCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCgga                                     <  2:67952/251‑1 (MQ=255)                                                                                                                                                                                                    cATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCggag                                    <  1:386083/251‑1 (MQ=255)                                                                                                                                                                                                         tCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCggagga                                  <  1:430776/248‑1 (MQ=255)                                                                                                                                                                                                                    tACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGgcgc                    >  1:404815/1‑251 (MQ=255)                                                                                                                                                                                                                            cATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGgcgc            <  1:226268/251‑1 (MQ=255)                                                                                                                                                                                                                                 gACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAAAAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCTGAGGAATGGTCGGCGGCGCACTGGCGCTgggg       >  1:291182/1‑251 (MQ=255)                                                                                                                                                                                                                                    aGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTg    <  1:526588/251‑1 (MQ=255)                                                                                                                                                                                                                                       cccAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGc  >  2:232376/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             TTGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAT‑AGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCG‑TTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAA‑GCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGC  >  NC_000913/694695‑695174

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