Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,792,977	G→A	N194N (AAC→AAT)	btuD ←	vitamin B12 ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,792,977	0	G	A	100.0%	163.4 / NA	50	N194N (AAC→AAT)	btuD	vitamin B12 ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/21);  total (29/21)

ATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGA  >  NC_000913/1792730‑1793214                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              ttaagtgatattggttgtctggagattcagggggccagtctaGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTg                                                                                                                                                                                                                                            >  1:75227‑M2/43‑251 (MQ=255)   aagcgatataggttgtctggagattcagggggccagtataACGTCGTAATTCCCCGCTTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTTTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGACCGATGCGCATGACGCAATGTGTGATTg                                                                                                                                                                                                                                            >  2:456750‑M2/41‑249 (MQ=255)         tattggttgtctggagattcagggggccagtctaGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTg                                                                                                                                                                                                                                    >  1:347913‑M2/35‑251 (MQ=255)              gttgtctggagattcagggggccagtctaGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGc                                                                                                                                                                                                                               >  1:300745‑M2/30‑251 (MQ=255)                                     agtctaGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGc                                                                                                                                                                                                        >  2:277069‑M2/7‑251 (MQ=255)                                        ctaGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGCTGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCtgg                                                                                                                                                                                                     >  2:434731‑M2/4‑250 (MQ=255)                                          aGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCAAGACGGCGAAACTTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGac                                                                                                                                                                                                    >  2:22118/1‑250 (MQ=255)                                                  aaTTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAgcgcgc                                                                                                                                                                                           >  2:319835/1‑251 (MQ=255)                                                             tttCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAtttt                                                                                                                                                                                 <  1:175787/250‑1 (MQ=255)                                                                 aaGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCt                                                                                                                                                                             >  1:266001/1‑250 (MQ=255)                                                                          aTGNAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACttt                                                                                                                                                                   <  2:414807/251‑1 (MQ=255)                                                                           tGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAAAGCACTttg                                                                                                                                                                  >  1:275268/1‑251 (MQ=255)                                                                            gAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTttgt                                                                                                                                                                 <  1:340499/251‑1 (MQ=255)                                                                                    aTGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAAc                                                                                                                                                         >  2:190777/1‑251 (MQ=255)                                                                                      ggTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTAAGCAACAACCCCCGAGGCGCTTGACCCATGGTGGGATTGGGTTCGTGATTGCCCTTCCCATTCCCAAGTCTTTGCTGACACAGCGGGCCAAGAATTTTGCTTAACCCACttttttccccaccac                                                                                                                                                       >  1:336622/1‑239 (MQ=255)                                                                                      ggTCGAAATCAGCATTATGTGACCTTCGATATCCAGACGGCGAAAATTCATACCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACAt                                                                                                                                                       >  2:17069/1‑251 (MQ=255)                                                                                       gTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATc                                                                                                                                                      <  1:183824/251‑1 (MQ=255)                                                                                               cAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTg                                                                                                                                               <  2:26775/250‑1 (MQ=255)                                                                                                aGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGtt                                                                                                                                             >  2:10940/1‑251 (MQ=255)                                                                                                  cATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCCCTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGGCACAGCGCGCTCAGAATTTTGGCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCa                                                                                                                                           <  1:431401/251‑1 (MQ=255)                                                                                                                cGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCaag                                                                                                                              <  1:92235/250‑1 (MQ=255)                                                                                                                cGATATCAAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCGCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCtcatca                                                                                                                                >  1:159579/1‑248 (MQ=255)                                                                                                                 gATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCaaga                                                                                                                             >  1:93873/1‑250 (MQ=255)                                                                                                                    atCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAagc                                                                                                                          >  1:90643/1‑250 (MQ=255)                                                                                                                            ggCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTgg                                                                                                                 <  1:328224/251‑1 (MQ=255)                                                                                                                             gCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGc                                                                                                                >  2:92061/1‑251 (MQ=255)                                                                                                                              cgaaaagttattccATAGGCCTGCGCCAGCTTTGGCGGCGTGAGCACCTTTTCCCTGCGTCTACTGGCCAGCATTTTTCCACCTTTTAGCAACCACTTTCGATGTGCACGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGTTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTCGTTGCGCAACATCAAGACTGTTCATCGGCTCATAAAGAAGATGCAATTGGcc                                                                                                               >  2:147827/13‑251 (MQ=255)                                                                                                                                      cATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGa                                                                                                        <  1:72135/250‑1 (MQ=255)                                                                                                                                                              tGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAac                                                                                >  1:335613/1‑250 (MQ=255)                                                                                                                                                              tGGCGGAGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCNac                                                                                >  2:411687/1‑250 (MQ=255)                                                                                                                                                                     gTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCgc                                                                        >  2:170056/1‑251 (MQ=255)                                                                                                                                                                              tctTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTTCAACaccacccccgcaacacc                                                               >  2:102729/1‑240 (MQ=255)                                                                                                                                                                                ttCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGt                                                              >  2:253299/1‑250 (MQ=255)                                                                                                                                                                                  cccTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGCCTGTTCATCGGCTCATCAAGAAGCAGCACTTGGCCTGCGGGa                                                                                                        >  2:263114/1‑206 (MQ=255)                                                                                                                                                                                  cccTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGa                                                                                                        <  1:263114/206‑1 (MQ=255)                                                                                                                                                                                  cccTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTAc                                                            <  2:298427/250‑1 (MQ=255)                                                                                                                                                                                           cccTGGCCAGCATTTTTCCACCCTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGGGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCAt                                                  <  2:275268/249‑1 (MQ=255)                                                                                                                                                                                               ggCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAATCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTc                                                <  1:411687/249‑1 (MQ=255)                                                                                                                                                                                               ggCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCa                                               <  1:110221/250‑1 (MQ=255)                                                                                                                                                                                                gCCAGCATTTTTTCCCCCTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCAcc                                             <  2:336622/251‑1 (MQ=255)                                                                                                                                                                                                 ccAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGGGTGATTTGCAACACCCCCGAAGCAAGACGCACGGGTTTCCATTCaacg                                            >  2:337265/1‑248 (MQ=255)                                                                                                                                                                                                        ttttCCTCCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGGGACCGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGGGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGGTCAGCGGCTCATCAAGAAGCAGCAATTGTCCTGCGGGATTGGCTTGTGGTGGGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGaa                                      <  1:22118/250‑1 (MQ=255)                                                                                                                                                                                                             cACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGa                                <  1:10940/251‑1 (MQ=255)                                                                                                                                                                                                                       cAACCACGCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCTGAAAGTTGATTGGTGCTAc                      <  2:244082/251‑1 (MQ=255)                                                                                                                                                                                                                              gCCCGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTaactccg                 >  2:114241/1‑243 (MQ=255)                                                                                                                                                                                                                                 cGATGCGCATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGttt            >  1:81703/1‑251 (MQ=255)                                                                                                                                                                                                                                       gcATGACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGtcatc      <  1:253299/251‑1 (MQ=255)                                                                                                                                                                                                                                           gACGCAATGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACTCAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAaga  <  1:190777/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaTGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCTCAACATCAAGACTGTTCATCGGCtc                                                                                                                                    <  2:420027/116‑1 (MQ=255)                                                                                                                                                                                                                                                aaTGTGTGATTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCtc                                                                                                                                    >  1:420027/1‑116 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              ATTCCGTACGGTTCCTCTGTTTTTTATTCTTGCATTAATTTAGCGTCGTAATTACCCGATTTTCAAGATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCCGGAAAGTTGATTGGTGCTACGTCCGAGTTTGTCATCAAGA  >  NC_000913/1792730‑1793214

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