Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	365,742	Δ1 bp	coding (564/3075 nt)	lacZ ←	beta‑D‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	365,741	0	C	.	100.0%	319.8 / NA	65	coding (565/3075 nt)	lacZ	beta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (35/30);  total (35/30)

GAAATCATCATTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGC  >  NC_000913/365598‑365886                                                                                                                                                |                                                                                                                                                  gAAATCATCATTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCAcc                                                                                                                                            <  1:852217/150‑1 (MQ=255) gAAATCATCATTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACcg                                                                                                                                           <  2:996453/151‑1 (MQ=255)     tcatcaTTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACcgcg                                                                                                                                         >  1:762900/1‑149 (MQ=255)     tcatcaTTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACcgcg                                                                                                                                         >  1:471896/1‑149 (MQ=255)      catcatTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAgg                                                                                                                                      <  2:681967/151‑1 (MQ=255)      catcatTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAgg                                                                                                                                      <  1:147490/151‑1 (MQ=255)             aaaGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGtttt                                                                                                                               >  2:484063/1‑151 (MQ=255)              aaGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTtc                                                                                                                              <  1:890988/151‑1 (MQ=255)                  gAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCg                                                                                                                          >  2:210934/1‑151 (MQ=255)                              ggAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAAt                                                                                                              <  2:745076/151‑1 (MQ=255)                               gAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATg                                                                                                             <  1:555203/151‑1 (MQ=255)                                 aaTCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATgcg                                                                                                           >  1:491446/1‑151 (MQ=255)                                       tGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAgg                                                                                                     >  1:1061831/1‑151 (MQ=255)                                           ttGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCa                                                                                                  >  2:750979/1‑150 (MQ=255)                                             gtgtAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCaaa                                                                                                <  2:1028010/150‑1 (MQ=255)                                                  gTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGa                                                                                          <  2:921297/151‑1 (MQ=255)                                                     ggTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACgg                                                                                       <  1:280769/151‑1 (MQ=255)                                                     ggTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACgg                                                                                       >  1:1082117/1‑151 (MQ=255)                                                      gTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGc                                                                                      >  1:163681/1‑151 (MQ=255)                                                        ttATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCa                                                                                     >  1:210763/1‑150 (MQ=255)                                                             cagcaACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGAc                                                                               <  1:1065757/151‑1 (MQ=255)                                                                caACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGt                                                                            >  2:1632/1‑151 (MQ=255)                                                                caACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGt                                                                            <  1:856155/151‑1 (MQ=255)                                                                caACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGt                                                                            <  2:213789/151‑1 (MQ=255)                                                                 aaCGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTc                                                                           <  2:606001/151‑1 (MQ=255)                                                                   cGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCt                                                                         <  1:704986/151‑1 (MQ=255)                                                                   cGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCt                                                                         <  2:705667/151‑1 (MQ=255)                                                                        cGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGTAAACGACTGTCCTGGCCg                                                                    <  1:231701/151‑1 (MQ=255)                                                                          tCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGt                                                                   >  2:1054390/1‑150 (MQ=255)                                                                           cACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTaa                                                                 <  2:846394/151‑1 (MQ=255)                                                                           cACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTa                                                                  <  1:647833/150‑1 (MQ=255)                                                                            aCGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTTAGACGGCAAACGACTGTCCTGGCCGTaa                                                                 >  1:881306/1‑150 (MQ=255)                                                                            aCGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAAc                                                                <  2:113872/151‑1 (MQ=255)                                                                             cGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAAcc                                                               >  1:491738/1‑151 (MQ=255)                                                                              ggAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAAcc                                                               >  1:418820/1‑150 (MQ=255)                                                                                   aTGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCa                                                         <  2:142532/151‑1 (MQ=255)                                                                                         cTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGccc                                                   >  2:875/1‑151 (MQ=255)                                                                                         cTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGccc                                                   >  2:920636/1‑151 (MQ=255)                                                                                          tCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGccc                                                   >  1:1044686/1‑150 (MQ=255)                                                                                          tCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGACAAACGACTGTCCTGGCCGTAACCGACCCAGCGccc                                                   >  1:280811/1‑150 (MQ=255)                                                                                           cATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGccc                                                   >  1:512720/1‑149 (MQ=255)                                                                                           cATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCg                                                  <  1:780594/150‑1 (MQ=255)                                                                                              ccGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTg                                               <  2:17163/150‑1 (MQ=255)                                                                                              ccGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGc                                              <  1:781057/151‑1 (MQ=255)                                                                                                gCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGcac                                            >  2:309989/1‑151 (MQ=255)                                                                                                 ccACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGcacc                                           <  2:230055/151‑1 (MQ=255)                                                                                                     atatCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGcaccac                                         >  2:412058/1‑149 (MQ=255)                                                                                                       atCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGcaccac                                         >  2:320992/1‑147 (MQ=255)                                                                                                       atCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGcacca                                          >  1:834885/1‑146 (MQ=255)                                                                                                       atCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGAt                                     <  2:171422/151‑1 (MQ=255)                                                                                                         ccTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATg                                    <  2:262128/150‑1 (MQ=255)                                                                                                         ccTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATg                                    >  2:139480/1‑150 (MQ=255)                                                                                                          cTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGGCCAGGCCAAATTCAGACCGCAACCGACTGTCCTGGCCGCAACCGACCCAGGCCCCGTTGCACCACAGGTg                                    >  1:699930/1‑149 (MQ=255)                                                                                                               cTTCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGcc                             >  2:726016/1‑151 (MQ=255)                                                                                                                ttCCAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGcc                             >  1:385808/1‑150 (MQ=255)                                                                                                                  ccAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAAc                                <  2:461235/145‑1 (MQ=255)                                                                                                                  ccAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCg                            >  2:672256/1‑149 (MQ=255)                                                                                                                   cAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGTTGAAACGCCGAGt                         >  2:613919/1‑151 (MQ=255)                                                                                                                   cAGATAACTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAg                          >  1:922547/1‑150 (MQ=255)                                                                                                                        aaCTGCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAAc                     >  2:50967/1‑150 (MQ=255)                                                                                                                            gCCGTCACTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAACGCCAt                <  2:322053/151‑1 (MQ=255)                                                                                                                                  aCTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCaaaa           >  2:1063057/1‑150 (MQ=255)                                                                                                                                  aCTCCAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCaaaa           >  1:1102228/1‑150 (MQ=255)                                                                                                                                     ccAGCGCAGCA‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAaataa       <  2:29125/151‑1 (MQ=255)                                                                                                                                          gcagca‑CATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCACATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTcgc  >  1:6422/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                  GAAATCATCATTAAAGCGAGTGGCAACATGGAAATCGCTGATTTGTGTAGTCGGTTTATGCAGCAACGAGACGTCACGGAAAATGCCGCTCATCCGCCACATATCCTGATCTTCCAGATAACTGCCGTCACTCCAGCGCAGCACCATCACCGCGAGGCGGTTTTCTCCGGCGCGTAAAAATGCGCTCAGGTCAAATTCAGACGGCAAACGACTGTCCTGGCCGTAACCGACCCAGCGCCCGTTGCACCACAGATGAAACGCCGAGTTAACGCCATCAAAAATAATTCGC  >  NC_000913/365598‑365886

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