Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,381,504	T→C	E272G (GAA→GGA)	frdA ←	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,381,504	0	T	C	100.0%	275.2 / NA	74	E272G (GAA→GGA)	frdA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/41);  total (33/41)

CCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGC  >  NC_000913/4381357‑4381650                                                                                                                                                    |                                                                                                                                                   ccGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCg                                                                                                                                                 <  1:607026/151‑1 (MQ=255)   gAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCggg                                                                                                                                               >  1:985691/1‑151 (MQ=255)      gTGACACAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGATGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTATCTCGGCTCGCCCAGCGGAGTTCCCgggtc                                                                                                                                             <  2:1093715/150‑3 (MQ=255)        gACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCAt                                                                                                                                          >  2:559448/1‑151 (MQ=255)        gACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCAt                                                                                                                                          <  2:174551/151‑1 (MQ=255)            cAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCg                                                                                                                                      <  2:185531/151‑1 (MQ=255)              aGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTa                                                                                                                                    <  2:1016269/151‑1 (MQ=255)              aGTCGAGATAAACCACAGCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGcc                                                                                                                                       <  2:858707/148‑1 (MQ=255)               gTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTaa                                                                                                                                   <  1:328240/151‑1 (MQ=255)                     aTAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTg                                                                                                                              <  2:1194027/150‑1 (MQ=255)                       aaaCCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGc                                                                                                                             >  1:661027/1‑149 (MQ=255)                       aaaCCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAg                                                                                                                           <  1:976874/151‑1 (MQ=255)                         aCCACACCGCCACGCGGCGTGGAGATGGTGTTGCCTTTCCGCCATTCGTGCCAGAGGGCCTGAGAGACTTTGTCGCGTGAACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGAt                                                                                                                         <  2:707342/151‑1 (MQ=255)                           cacaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATa                                                                                                                        >  1:1080945/1‑150 (MQ=255)                             caTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACg                                                                                                                     <  2:887509/151‑1 (MQ=255)                                cGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGt                                                                                                                   <  2:1051742/150‑1 (MQ=255)                                      gcggcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAtt                                                                                                            <  1:1212050/151‑1 (MQ=255)                                        ggcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAtttt                                                                                                          >  2:408726/1‑151 (MQ=255)                                          cgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAttttt                                                                                                         >  2:585976/1‑150 (MQ=255)                                           gTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgt                                                                                                       <  2:830049/151‑1 (MQ=255)                                           gTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgt                                                                                                       <  2:668084/151‑1 (MQ=255)                                           gTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgt                                                                                                       >  1:102580/1‑151 (MQ=255)                                            tGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgt                                                                                                       <  1:921215/150‑1 (MQ=255)                                              gagaTGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgttgg                                                                                                    <  2:1206331/151‑2 (MQ=255)                                                  tGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGAcc                                                                                                  >  1:1104003/1‑149 (MQ=255)                                                   ggTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGGGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGATTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGa                                                                                               <  2:759672/151‑1 (MQ=255)                                                   ggTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGATTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGa                                                                                               <  2:793198/151‑1 (MQ=255)                                                   ggTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGa                                                                                               <  2:256801/151‑1 (MQ=255)                                                            tttACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAAt                                                                                             <  1:765452/144‑1 (MQ=255)                                                            tttACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAAt                                                                                             >  2:765442/1‑144 (MQ=255)                                                            tttACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCt                                                                                      >  1:477560/1‑151 (MQ=255)                                                                   cATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACcgc                                                                               >  2:375771/1‑151 (MQ=255)                                                                       cGTGCCAGAAGGCCGGAGAGGCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGATTCCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCGATTTTTGTTGACCAGAATACCGCCTTCACCGCgg                                                                             <  2:806263/149‑1 (MQ=255)                                                                         tGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGGAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAAc                                                                         >  2:176840/1‑151 (MQ=255)                                                                             aGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGTGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGATTACCGCCTTCACCGCGGCAACCtt                                                                      <  2:205029/150‑1 (MQ=255)                                                                              gAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCg                                                                    <  1:18895/151‑1 (MQ=255)                                                                               aaGGCCTGAGAGGCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCgg                                                                   <  2:186086/151‑1 (MQ=255)                                                                                  gCCTGAGAGACTTTGTCGAGTGGACCCAGTTCCATATATTTGTTTTTGGGGTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGACTTCACCGCGGCAACCTTCGGtca                                                                <  2:360600/151‑1 (MQ=255)                                                                                    cTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtcatc                                                              >  1:257051/1‑151 (MQ=255)                                                                                    cTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtcat                                                               >  2:135028/1‑150 (MQ=255)                                                                                     tGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtcatc                                                              >  1:1095832/1‑150 (MQ=255)                                                                                        gagaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAgg                                                           >  2:319275/1‑150 (MQ=255)                                                                                        gagaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGa                                                          <  1:401344/151‑1 (MQ=255)                                                                                         agaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGa                                                          <  2:683706/150‑1 (MQ=255)                                                                                               tGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCgg                                                    <  1:771635/150‑1 (MQ=255)                                                                                               tGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCgg                                                    >  1:195933/1‑150 (MQ=255)                                                                                                  cgcgTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGaa                                                  <  2:892786/149‑1 (MQ=255)                                                                                                  cgcgTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAAcc                                                <  2:290969/151‑1 (MQ=255)                                                                                                     gtgtACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTgg                                             <  2:49281/151‑1 (MQ=255)                                                                                                     gTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTg                                              <  1:1098083/150‑1 (MQ=255)                                                                                                        gACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGc                                            >  1:301693/1‑149 (MQ=255)                                                                                                         aCCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGa                                         <  1:179308/151‑1 (MQ=255)                                                                                                          cccAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGAc                                        >  1:1029916/1‑151 (MQ=255)                                                                                                           ccAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGAcc                                       >  1:819052/1‑151 (MQ=255)                                                                                                            cAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCg                                      >  2:185700/1‑151 (MQ=255)                                                                                                                     atatTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTg                              >  2:309364/1‑150 (MQ=255)                                                                                                                          tGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTg                        >  1:596340/1‑151 (MQ=255)                                                                                                                           gTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTTGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTg                        >  1:170484/1‑150 (MQ=255)                                                                                                                              tttCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAAc                     <  1:19132/150‑1 (MQ=255)                                                                                                                                tCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGa                   >  2:545755/1‑150 (MQ=255)                                                                                                                                 cGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGaa                  <  1:1187044/150‑1 (MQ=255)                                                                                                                                 cGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAAt                 <  1:717202/151‑1 (MQ=255)                                                                                                                                 cGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAAt                 <  1:527392/151‑1 (MQ=255)                                                                                                                                  ggCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAAtt                <  2:110546/151‑1 (MQ=255)                                                                                                                                   gCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTc               <  2:25689/151‑1 (MQ=255)                                                                                                                                      cGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCAt            >  1:23000/1‑151 (MQ=255)                                                                                                                                        cccAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGt          >  2:1038016/1‑151 (MQ=255)                                                                                                                                           aGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCAc       >  2:516119/1‑151 (MQ=255)                                                                                                                                           aGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCAc       >  1:993021/1‑151 (MQ=255)                                                                                                                                             cGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc     >  1:506812/1‑151 (MQ=255)                                                                                                                                              ggAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc     >  1:88135/1‑150 (MQ=255)                                                                                                                                               gAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc     >  1:433612/1‑149 (MQ=255)                                                                                                                                                aGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGc  <  1:1150277/151‑1 (MQ=255)                                                                                                                                                aGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGc  >  1:787546/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGC  >  NC_000913/4381357‑4381650

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