Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,986,686	C→T	*399* (TAG→TAA)	hemY ←	putative protoheme IX synthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,986,686	0	C	T	100.0%	234.7 / NA	73	*399* (TAG→TAA)	hemY	putative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (41/32);  total (41/32)

CATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGTTGTTT  >  NC_000913/3986551‑3986817                                                                                                                                        |                                                                                                                                    cATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAattac                                                                                                                                  >  1:438350/1‑139 (MQ=255)   tGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtg                                                                                                                                >  2:321747/1‑139 (MQ=255)    ggACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACtgt                                                                                                                               >  2:445742/1‑139 (MQ=255)        acaACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCg                                                                                                                           <  1:321747/139‑1 (MQ=255)          aaCGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                                         <  1:369595/139‑1 (MQ=255)                    gAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCa                                                                                                               >  1:878344/1‑139 (MQ=255)                    gAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCa                                                                                                               >  2:647444/1‑139 (MQ=255)                                       aTACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTc                                                                                            >  1:419086/1‑139 (MQ=255)                                           aGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGa                                                                                        >  2:111473/1‑139 (MQ=255)                                            gACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGAc                                                                                       >  2:659199/1‑139 (MQ=255)                                            gACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGAc                                                                                       >  1:457969/1‑139 (MQ=255)                                            gACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGAc                                                                                       >  2:1081978/1‑139 (MQ=255)                                             aCCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                                         >  2:133723/1‑104 (MQ=255)                                             aCCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCggg                                                                                                                         <  1:133723/104‑1 (MQ=255)                                             aCCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                      >  1:30151/1‑139 (MQ=255)                                               cTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCa                                                                                    <  2:878344/139‑1 (MQ=255)                                                       aCGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCg                                                                                              >  1:1016529/1‑121 (MQ=255)                                                       aCGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCg                                                                                              <  2:1016529/121‑1 (MQ=255)                                                       aCGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                      <  1:107448/129‑1 (MQ=255)                                                       aCGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACg                                                                                      >  2:107448/1‑129 (MQ=255)                                                        cGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAg                                                                           <  1:133681/139‑1 (MQ=255)                                                         gCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGc                                                                          <  1:1081978/139‑1 (MQ=255)                                                            tGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGcttc                                                                       >  2:427188/1‑139 (MQ=255)                                                              cTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCg                                                                                              <  1:332850/114‑1 (MQ=255)                                                              cTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCg                                                                                              >  2:332850/1‑114 (MQ=255)                                                                  gTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCgg                                                                 <  2:30151/139‑1 (MQ=255)                                                                   taataaGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGc                                                                >  2:852529/1‑139 (MQ=255)                                                                    aataaGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCt                                                               <  1:659199/139‑1 (MQ=255)                                                                       aaGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTtgt                                                            <  2:438350/139‑1 (MQ=255)                                                                         gagCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGc                                                          <  1:445742/139‑1 (MQ=255)                                                                            cAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGt                                                       >  1:506743/1‑139 (MQ=255)                                                                             aGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCa                                                                                                               <  2:799726/82‑1 (MQ=255)                                                                             aGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCa                                                                                                               >  1:799726/1‑82 (MQ=255)                                                                                    tttttATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACAt                                                                                                     >  2:1037358/1‑85 (MQ=255)                                                                                    tttttATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACAt                                                                                                     <  1:1037358/85‑1 (MQ=255)                                                                                         aTGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGAAGCTTCTTCCGGGTTGTGCAGTCTGTGCAGCGCGt                                          >  2:275142/1‑139 (MQ=255)                                                                                             aTCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGc                                      <  2:457969/139‑1 (MQ=255)                                                                                              tCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCa                                     >  1:98298/1‑139 (MQ=255)                                                                                               cAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGt                                                                                                                        <  1:220238/55‑1 (MQ=255)                                                                                               cAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGt                                                                                                                        >  2:220238/1‑55 (MQ=255)                                                                                               cAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGa                                                                                        >  1:781436/1‑87 (MQ=255)                                                                                               cAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGa                                                                                        <  2:781436/87‑1 (MQ=255)                                                                                                  gAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGcc                                 <  2:98298/139‑1 (MQ=255)                                                                                                      gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCa                                                                                    >  1:786288/1‑84 (MQ=255)                                                                                                      gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCa                                                                                    <  2:786288/84‑1 (MQ=255)                                                                                                      gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGc                                                                          <  2:502411/94‑1 (MQ=255)                                                                                                      gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGc                                                                          >  1:502411/1‑94 (MQ=255)                                                                                                      gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGt                                                       >  1:106599/1‑113 (MQ=255)                                                                                                      gCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGt                                                       <  2:106599/113‑1 (MQ=255)                                                                                                       ccccGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCg                            >  2:1024232/1‑139 (MQ=255)                                                                                                          cGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                           >  1:200420/1‑73 (MQ=255)                                                                                                          cGGAGGTGCTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTcg                                                                                           <  2:200420/73‑1 (MQ=255)                                                                                                                  cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGt                                                                                                           <  1:640345/49‑1 (MQ=255)                                                                                                                  cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGt                                                                                                           >  2:640345/1‑49 (MQ=255)                                                                                                                  cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGt                                                       >  2:42250/1‑101 (MQ=255)                                                                                                                  cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGt                                                       <  1:42250/101‑1 (MQ=255)                                                                                                                  cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGc                 >  1:345113/1‑139 (MQ=255)                                                                                                                  cTTGCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGc                 >  1:698120/1‑139 (MQ=255)                                                                                                                     gCCTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTc              >  1:75251/1‑139 (MQ=255)                                                                                                                       cTCCGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCg            >  1:896364/1‑139 (MQ=255)                                                                                                                         ccGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGa          >  2:911577/1‑139 (MQ=255)                                                                                                                         ccGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGa          >  1:396274/1‑139 (MQ=255)                                                                                                                          cGGGTGAGAAGGAATTTCTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTc                                                      <  2:704115/94‑1 (MQ=255)                                                                                                                          cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                             <  2:566777/119‑1 (MQ=255)                                                                                                                          cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGc                             >  1:566777/1‑119 (MQ=255)                                                                                                                          cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGAc         >  2:578267/1‑139 (MQ=255)                                                                                                                          cGGGTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTACAGTc                                                      >  1:704115/1‑94 (MQ=255)                                                                                                                           gggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACg        <  2:345113/139‑1 (MQ=255)                                                                                                                            ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTtgt                                                            <  2:27687/86‑1 (MQ=255)                                                                                                                            ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTtgt                                                            >  1:27687/1‑86 (MQ=255)                                                                                                                            ggTGAGAAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACgt       <  1:911577/139‑1 (MQ=255)                                                                                                                                 gaAGGAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGTTGttt  <  1:427188/139‑1 (MQ=255)                                                                                                                                     gAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTcc                                                                   >  2:809670/1‑70 (MQ=255)                                                                                                                                     gAATTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTcc                                                                   <  1:809670/70‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    CATGGACACAACGTTGAGTGAAGCACCCACTTGTTGTCATACAGACCTGTTTTAACGCCTGCTCCGTAATAAGAGCAGGCGTTTTTTTATGTATCAGGAAGGCCCCGGAGGTGCTTGCCTCCGGGTGAGAAGGAACTACTGTGGCGGGTTATTCTGCAACGTTAACATCAAACCGTCGCGACGCATAGCTGCAGCTTCTTCCGGCTTGTGCAGTCTGTCCAGCGCGTCGGCAAGCCATGCGTAATCGTAGGCGTCCGGACGTTGTTT  >  NC_000913/3986551‑3986817

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