Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	74,872	C→T	L203L (CTG→CTA)	thiB ←	thiamine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	74,872	0	C	T	100.0%	131.4 / NA	38	L203L (CTG→CTA)	thiB	thiamine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/21);  total (17/21)

GATACATCCAGTTGCCGGTTGGGATCGCATTCTGGAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAA  >  NC_000913/74634‑75097                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   gATACATCCAGTTGCCGGTTGGGATCGCATTCTGGAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGTCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCt                                                                                                                                                                                                                        <  2:264403/251‑1 (MQ=255)      aTCCAGTTGCCGGTTGGGATCGCATTCTGGAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCt                                                                                                                                                                                                                   >  2:237673/1‑251 (MQ=255)         cAGTTGCCGGTTGGGATCGCATTCTGGAAAGCCGGAGAAACCATAAACGGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGGAGGCTTCGCTCCa                                                                                                                                                                                                                <  2:324345/251‑1 (MQ=255)                    tGGGATCGCATTCTGGAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGgtg                                                                                                                                                                                                      <  2:335997/250‑1 (MQ=255)                        aTCGCATTCTGGAAAGCCGGAGAAACCAGAAACTGGAGGAATGTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCTTAGGCTTCGCTCCAGCCTTTGgtgaccgt                                                                                                                                                                                                 <  2:128557/251‑1 (MQ=255)                         tCGCATTCTGGAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTCGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGgtgaccgtg                                                                                                                                                                                                >  2:148804/1‑251 (MQ=255)                             aTTCTGGAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGtgaccgtgaccg                                                                                                                                                                                            <  2:396097/251‑1 (MQ=255)                                   gAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTttct                                                                                                                                                                                      <  2:164735/251‑1 (MQ=255)                                             gaAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCATATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGttt                                                                                                                                                                            <  1:260892/251‑1 (MQ=255)                                                       aCTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTCCTGGCAGCGGTGCAGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGCATGTGATAAGACGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGTCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCt                                                                                                                                                                  >  1:253982/1‑251 (MQ=255)                                                          ggaggaATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTtgt                                                                                                                                                               <  1:375497/251‑1 (MQ=255)                                                                  ttttGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCGGCCAGGCTTTTGGGGCGt                                                                                                                                                        <  1:362542/250‑1 (MQ=255)                                                                    ttGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATATTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCTCGCCAGCCTGTGGGGGCCGTGACCGTTTGCTTCGCCAGTTTCTGCCAGGCTTCTGGGGCGTCa                                                                                                                                                      <  2:62536/250‑1 (MQ=255)                                                                    ttGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCAt                                                                                                                                                     >  1:48542/1‑251 (MQ=255)                                                                                cGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTGCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACtt                                                                                                                                         <  1:38040/251‑1 (MQ=255)                                                                                            ggCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATa                                                                                                                             <  1:168170/251‑1 (MQ=255)                                                                                                                  tCCATTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGGGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGGTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCAGCTATAGCAACAGACCCAGCCCCGgt                                                                                                         <  2:178725/249‑1 (MQ=255)                                                                                                                  tCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGgtg                                                                                                        <  1:58336/250‑1 (MQ=255)                                                                                                                          cAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTgc                                                                                               <  2:87883/251‑1 (MQ=255)                                                                                                                           aGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTgcg                                                                                              <  1:3190/251‑1 (MQ=255)                                                                                                                                gTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGAt                                                                                         >  2:95771/1‑251 (MQ=255)                                                                                                                                  gACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATcc                                                                                       >  2:294502/1‑251 (MQ=255)                                                                                                                                               gTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTAACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCAccc                                                                          <  2:280459/251‑1 (MQ=255)                                                                                                                                                       cggcgTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGttt                                                                  >  1:114404/1‑251 (MQ=255)                                                                                                                                                                   ctttcttCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTTGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGa                                                                                    >  2:281353/1‑221 (MQ=255)                                                                                                                                                                   ctttcttCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGa                                                                                    <  1:281353/221‑1 (MQ=255)                                                                                                                                                                              tcGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCttt                                           >  1:924/1‑251 (MQ=255)                                                                                                                                                                                  gaATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTgcg                                                                                              <  1:66772/196‑1 (MQ=255)                                                                                                                                                                                  gaATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTgcg                                                                                              >  2:66772/1‑196 (MQ=255)                                                                                                                                                                                        gATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGc                                                                                                                                                                        >  2:173944/1‑116 (MQ=255)                                                                                                                                                                                        gATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGc                                                                                                                                                                        <  1:173944/116‑1 (MQ=255)                                                                                                                                                                                           aaGCCGGAGAGGCGGTGTAACTCAGTACCAGACCGCTTTCCCCTTTTAACAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGc                              >  2:21874/1‑251 (MQ=255)                                                                                                                                                                                              ccGGAGTGGTGGTGTAACTCAGTACCAGTTCGCTTTCCCCTTTAAAAAATGGGCCGTGGGCTTCGCTCCAGCTTTTGGTGACCGTGCCCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTCCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGGGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCggg                           <  2:48542/251‑1 (MQ=255)                                                                                                                                                                                                ggAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCggggg                         >  2:406055/1‑249 (MQ=255)                                                                                                                                                                                                                taaGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTACGCTCCAGCCTTTGGTGACCGTGACAGTTTTCTACGCCAGTATCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGGGTACTGGTGCGCGGATCCTGAGAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑tttgtttaattgttcttt        >  1:161354/3‑236 (MQ=255)                                                                                                                                                                                                                tCAGTACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGATCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGGTTTTTcctttgtttttt         >  1:161360/1‑240 (MQ=255)                                                                                                                                                                                                                    tACCAGATCGCTTTCACCTTTTAAAAATAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCTTCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑tgtttttgtttttttttgtttt    >  2:286959/1‑234 (MQ=255)                                                                                                                                                                                                                       cAGATCGCTTTCAACTTTTAAAAATAGGCAGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑CTTTTttgtttgtgtttttcctaa  >  1:131749/1‑239 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                   GATACATCCAGTTGCCGGTTGGGATCGCATTCTGGAAAGCCGGAGAAACCATAAACTGGAGGAATTTTTGCGCCAGCTCCGGCTGCTTGCTGGCAGCGGTGCGGGCGGCGACTTCCACTTGCAGATAGTGACCTTCGCTGAAGTTCGCGGCGGCGTAGTTATCTTTCTTCTCTTCGAGAATGTGATAAGCCGGAGAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGG‑TTTTTCAGTTTGTTCTTGTCATAA  >  NC_000913/74634‑75097

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