Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	109,960	G→A	R561H (CGT→CAT)	secA →	protein translocation ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	109,960	0	G	A	100.0%	122.9 / NA	43	R561H (CGT→CAT)	secA	protein translocation ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/26);  total (17/26)

CCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATTCG  >  NC_000913/109727‑110203                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    ccACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTTGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCcgtcgt                                                                                                                                                                                                                                     <  1:201300/250‑1 (MQ=255)    cGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCg                                                                                                                                                                                                                                 <  1:386478/251‑1 (MQ=255)                   gATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGcgcg                                                                                                                                                                                                                  <  1:272464/251‑1 (MQ=255)                      ttttgCTCAGGCAGGTTACCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGAGAGTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCGGAGCAAAGTGAAAAAAGTAGAGCCGACTGGCAGGTACGGCACGATGCGGTACTGGAAGCGGGGGGCCTGCATATCATCGGTAACGAGCATCACGAAACCAGTCGTATCGACAACCAGTTGCGCGGt                                                                                                                                                                                                                <  2:59578/248‑1 (MQ=255)                       gttgCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTc                                                                                                                                                                                                               <  2:111993/250‑1 (MQ=255)                       ggtgCTCAGGCAGGTTAGCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCg                                                                                                                                                                                                              <  2:239461/249‑1 (MQ=255)                                         ccGGCTGCGGTGACTATCGCGACCAATAGTGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGa                                                                                                                                                                                            <  1:416266/251‑1 (MQ=255)                                                              aCCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTa                                                                                                                                                                       <  1:397578/251‑1 (MQ=255)                                                                                 gTACAGATTTTGGGCTCGGTGGTAGCTGGCAGGCAGAAGTGGCCGCGCTGGAAAAGCCGACCGCAGAGCAAAGTGAAAAAATTAAAGCCGACTGGCAGGTACGGCACGATGCGGGACTGGAAGCAGGTGGCCTGCATATCAGCGGTACCGAGCAGCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGGTCTTCCCGTTTCTACCTGTCGATGGAAGATgcg                                                                                                                                                    <  2:203034/251‑1 (MQ=255)                                                                                    cAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTGCTACCTGTCGATGGAAGATGCGCTg                                                                                                                                                 <  1:132067/251‑1 (MQ=255)                                                                                     aGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTTGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGa                                                                                                                                                <  1:191646/251‑1 (MQ=255)                                                                                             tgCTCGGAGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGAAGGGGTACTGGACGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAAACCCGTCGTATCGGTAACCAATTGCGCGGGCGTTCTGGACGTCAGGGGGGTGCTGGTTCTTCCCGTTTTTACCTTTCGATGGAAGATGCGCTGATGCGtgtt                                                                                                                                        >  2:337523/1‑248 (MQ=255)                                                                                              tcTCGGAGGTAGCTGGCAGGCAGAAGTTGGCGCGCTGGAAAATCCGCACGCAGAGCAAATTTAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGAATATGATCGGTACGGAGCATCACGAATCCCGTCGTTTCGGTAACCATGGGCGCGGTCGGTGTGGTCGTCAGGTGGATGCTGGTTCTTCCCGTTTTTACCGGTGGATGGAAGATGCGCTGATGCGTAtt                                                                                                                                        <  2:352682/249‑1 (MQ=255)                                                                                                     ggtAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCATAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAAGTGGCCTGCATATCATCGGTACCGAGCATCACGAATACCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTACGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGTTGGGAGATGCGCTGATGCGGATTTTTGCTTc                                                                                                                                >  1:262006/1‑251 (MQ=255)                                                                                                        aGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGGACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGATTccg                                                                                                                              >  1:30207/1‑250 (MQ=255)                                                                                                            ggaaggcagAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTccgacc                                                                                                                           <  2:13437/246‑1 (MQ=255)                                                                                                              caggcagAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGGCCCGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTAGCGATAACCAGTTGCGCGGTCGTTCTGgt                                                                                                                                                                                                       <  2:325282/171‑1 (MQ=255)                                                                                                              caggcagAAGTTGCCGCGCTGGAAAATCAGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCCGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTggc                                                                                                                                                                                                       >  1:325282/1‑170 (MQ=255)                                                                                                                ggcagAAGTTGCCGCGCTGTAAATACCGACCGCAGAGCAAATTGAAAAAAGTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGTTGGCCTGTATTGCAGCGGTTCCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTAt                                                                                                                     <  2:247591/251‑1 (MQ=255)                                                                                                                  cagAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACGGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATAGCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTTTCGATGGAAGATGCGCTGAGGCGTATTTTTGCTTCCGACCGAGTATcc                                                                                                                   <  2:233920/251‑1 (MQ=255)                                                                                                                       gTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATCGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCat                                                                                                              <  1:348152/251‑1 (MQ=255)                                                                                                                           ccGCGCTGGAAAATCCTACCGCAGAGCAAATGGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACGGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCTTCGTATCGATAGCCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCatgatg                                                                                                          <  1:151133/251‑1 (MQ=255)                                                                                                                             gcgcTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCAGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATAATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCCGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGCAGATGCGCTGGTGCGTATTTTTGCTTCCGCCCGAGAATCCGGCATGATGCg                                                                                                        >  1:22933/1‑251 (MQ=255)                                                                                                                               gcTGGAAAATCCGACCGCAGAGCAAAGTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTa                                                                                                      <  2:323215/251‑1 (MQ=255)                                                                                                                                       aTCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGTTAACCAGTTGCGCGGTCGTTCTGGTCGTCCGGGGGGTGCTGGTTCTTCCCGTTTCTACCTGTCGATGGGAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGCTGCGTTAACTgggg                                                                                              >  1:288604/1‑250 (MQ=255)                                                                                                                                         ccGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGGACGTCACGATGCGGGACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCCGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTCTTTTTGCTTCCGACCGAGTATCCGGCATGATGCGCAAACTGGGTAt                                                                                            >  1:413936/1‑251 (MQ=255)                                                                                                                                           gACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTTCGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGGCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGa                                                                                          <  1:43454/251‑1 (MQ=255)                                                                                                                                            aCCGCAGAGCAAATTGAAAAAATTAAAGCCGTCTGGCAGGTACGTCACGATGAGGGAATGGAAGCAAGGGGCCTGCATATCATCGGTTCCGAGAAACACGAAACCCGTCGCATAGATAACCAGTTGCGCGGTCGTTCTGGTCGACAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGGCGGTGGGAGATGCGCTGATGAGAATTTTTGCTTACGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGa                                                                                          >  2:14647/1‑250 (MQ=255)                                                                                                                                                    gCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCGCGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGAGGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCg                                                                                 <  1:381656/251‑1 (MQ=255)                                                                                                                                                                aaTTAAAGCCGACTGGCAGGTACGTCACGAGGCGGTACTGGAATCCAGTGGCCTGCATATCATCGGTACCGAACATCACGAATCCCCTCGGGTCGATAACCAGTTGCGCGGTCGGTCTGGGCGTCCGGGGGGTGCTGGATCTTCCCGTTTCTTCCTGTCGGTGGAGGATGCGATGTGGCGTATTTTTGCTTCCGACCGAGTTTACGGCATGGTGCGTAAACTGGGTATTAAACCAGGCGAAGccctttcac                                                                     >  2:35808/1‑244 (MQ=255)                                                                                                                                                                 aTTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGGATCGATAACCAGTTGCGCGGTCGTTCTGGTCGGCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAaca                                                                    <  2:372441/251‑1 (MQ=255)                                                                                                                                                                 aTTAAAGCCGACTGGCAGGAACGTCACGCTGCGGGACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGCGCATCACGAATCCCGTCGTATCGATAACCAGTTTCGAGGTCGTTCTGGTCGTCAGGGGGATGCTTGTTCTTCCCCGTTCTACCTTTCGATGGAAGATGCGCTGGTGCGTATTTTTTCTTCCGACCCAGGATCCGGCCTGATGGGTAAACTGGGTATGAAGCCCGGCgcaacccttgaaa                                                                     >  2:358263/1‑238 (MQ=255)                                                                                                                                                                   tAAAGCCGACTGGCAGGTACGGCACGATGCGGTACTGGAAGCAGGGGGCCTGCATATCCTCGGGACCGAGCATCCCGAATCCCGTCGTATCGATAACCAGTTGCGCGGGCGTTCTGGTCGTCAGGGGGATGCTGGGTCTTCCCGTTTCTACCTGTCGATGGAAGGTGCGCTGATGCGCTTTCTTGCTTCCGACCCAGTATCCGGCCTGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCCTTGAACAcc                                                                  >  2:103180/1‑251 (MQ=255)                                                                                                                                                                    aaaGCCGCCGGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTTCCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCTCGGTCGGTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACAccc                                                                 <  2:376682/251‑1 (MQ=255)                                                                                                                                                                        ccGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTgg                                                                                                <  1:318646/216‑1 (MQ=255)                                                                                                                                                                         cGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATAGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGGTGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGACATTGAACACCCGTggg                                                            >  1:49849/1‑251 (MQ=255)                                                                                                                                                                                     aCGTCACGATGCGGTACTGGAAGCAGGTGGCCCGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCg                                                 <  1:35808/250‑1 (MQ=255)                                                                                                                                                                                          aCGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGCCCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGcc                                           >  1:16509/1‑251 (MQ=255)                                                                                                                                                                                            gATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCa                                          >  1:406316/1‑250 (MQ=255)                                                                                                                                                                                                                   ccTGCATATCATCGGTATCGAGCATCACGAATCCCGGCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATGGCCAACGCCCAGCGTAAAGTTGAAAGcc                  <  1:358263/251‑1 (MQ=255)                                                                                                                                                                                                                                ggTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACAt     >  1:25422/1‑251 (MQ=255)                                                                                                                                                                                                                                   aCCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGACGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATTCg  >  1:266041/1‑251 (MQ=255)                                                                                                                                                                                                                                   aCCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATTCg  >  2:116391/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    CCACGCCAACGAAGCGGCGATTGTTGCTCAGGCAGGTTATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAGCCGTAACTTCGACATTCG  >  NC_000913/109727‑110203

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