Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,411	+C	coding (494/516 nt)	b4555 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,410	1	.	C	100.0%	304.7 / NA	94	A165P (GCT→CCT)	b4555	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (41/53);  total (41/53)

TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGC‑‑‑G‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATT‑G‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGAAA‑GCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492                                                                            |                                                                                         tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                     <  2:93170/93‑1 (MQ=255) tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                     >  2:6713870/1‑93 (MQ=255) tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                     >  2:1132533/1‑93 (MQ=255) tCTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCAGGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                     <  1:2981338/93‑1 (MQ=255) tCTGGAAGTGATGGCGGGCGGCGTGTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtt                                                                     <  1:6109615/93‑1 (MQ=255)  cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGttt                                                                    >  2:6652401/1‑92 (MQ=255)  cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGttt                                                                    >  1:990386/1‑92 (MQ=255)  cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGttt                                                                    <  1:1673509/93‑2 (MQ=255)  cTGGAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGttt                                                                    <  1:1567329/93‑2 (MQ=255)    ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttgg                                                                  >  2:6385580/1‑90 (MQ=255)    ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttgg                                                                  >  2:3952663/1‑90 (MQ=255)    ggAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttgg                                                                  >  2:5544515/1‑90 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                 <  1:5372490/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                 <  2:4370028/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                 <  2:6498533/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                 <  1:1500944/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                 <  2:4097734/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                 <  2:5650687/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggt                                                                 <  1:832008/93‑5 (MQ=255)     gAAGTGATGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGtttggg                                                                 >  2:1014114/1‑89 (MQ=255)           aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACggg                                                           <  1:6058621/93‑1 (MQ=255)           aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACggg                                                           <  2:3791411/93‑1 (MQ=255)           aTGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACggg                                                           <  1:3213125/93‑1 (MQ=255)            tGGCGGGCGGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGa                                                          <  2:1844340/93‑1 (MQ=255)             ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGAt                                                         <  2:5203299/93‑1 (MQ=255)             ggcgggcgGCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGAt                                                         <  2:3067848/93‑1 (MQ=255)              gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATg                                                        >  1:4191617/1‑93 (MQ=255)              gcgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATg                                                        >  1:5841435/1‑93 (MQ=255)               cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGCTAAGA‑‑‑AACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGa                                                       <  2:912291/93‑1 (MQ=255)               cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGa                                                       <  2:4998234/93‑1 (MQ=255)               cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGa                                                       <  1:31479/93‑1 (MQ=255)               cgggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGa                                                       >  1:2192480/1‑93 (MQ=255)                gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGaa                                                      <  2:3841035/93‑1 (MQ=255)                gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGaa                                                      <  1:47385/93‑1 (MQ=255)                gggcggCGTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGaa                                                      >  2:806217/1‑93 (MQ=255)                  gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacg                                                    >  2:274444/1‑93 (MQ=255)                  gcggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacg                                                    >  2:2993890/1‑93 (MQ=255)                   cggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacga                                                   <  2:3073690/93‑1 (MQ=255)                    ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgac                                                  <  1:5688350/93‑1 (MQ=255)                    ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgac                                                  >  2:5179197/1‑93 (MQ=255)                    ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgac                                                  <  2:2367591/93‑1 (MQ=255)                    ggcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgac                                                  <  2:1368824/93‑1 (MQ=255)                     gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgacg                                                 >  1:2108765/1‑93 (MQ=255)                     gcgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAacgacg                                                 <  2:1665924/93‑1 (MQ=255)                      cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGc                                                <  1:4782314/93‑1 (MQ=255)                      cgTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGc                                                >  2:777217/1‑93 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATTGG‑ATTGCGTTTGGTGTCGGGATGAACGACGcc                                               >  2:6280598/1‑93 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                               <  2:3562144/93‑1 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                               >  2:3784278/1‑93 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                               >  2:6297100/1‑93 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                               >  2:4869242/1‑93 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                               >  2:2100640/1‑93 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                               <  1:4708802/93‑1 (MQ=255)                       gTTTCAAAAGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGcc                                               >  1:3137376/1‑93 (MQ=255)                         ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGa                                             <  1:5667309/89‑1 (MQ=255)                         ttcaaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGa                                             <  1:6379746/89‑1 (MQ=255)                           caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                           <  1:2934883/91‑1 (MQ=255)                           caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                           <  2:5979161/91‑1 (MQ=255)                           caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                           <  2:3475704/91‑1 (MQ=255)                           caaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGaaa                                           <  2:1507336/91‑1 (MQ=255)                            aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                          <  2:5637763/92‑1 (MQ=255)                            aaaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAAt                                          <  2:3991366/92‑1 (MQ=255)                             aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                         >  1:2205579/1‑93 (MQ=255)                             aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                         >  1:4506931/1‑93 (MQ=255)                             aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                         >  1:721326/1‑93 (MQ=255)                             aaaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATg                                         <  2:6010772/93‑1 (MQ=255)                              aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                        <  1:6289377/93‑1 (MQ=255)                              aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                        <  2:5981706/93‑1 (MQ=255)                              aaGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGc                                        <  1:3895288/93‑1 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                       >  2:6020789/1‑93 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                       >  1:3236866/1‑93 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                       >  1:296548/1‑93 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                       >  2:6691263/1‑93 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                       >  2:2939531/1‑93 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                       >  2:2767946/1‑93 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCt                                       <  2:3264646/93‑1 (MQ=255)                               aGGCCATGCGCTGGAAGCGGTGGC‑‑‑GAAGAAACTGGGCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGAAGCCGAAATGcg                                       >  1:1083384/1‑92 (MQ=255)                                                                actgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg         >  1:1339363/5‑93 (MQ=255)                                                                actgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATggg         <  2:2395438/89‑1 (MQ=255)                                                                 ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGCTGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        >  2:4704086/4‑93 (MQ=255)                                                                 ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        >  1:5305292/4‑93 (MQ=255)                                                                 ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        <  2:3666627/90‑1 (MQ=255)                                                                 ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        <  1:1812278/90‑1 (MQ=255)                                                                 ctgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGc        <  2:6247674/90‑1 (MQ=255)                                                                  tgggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCa       <  1:1193761/91‑1 (MQ=255)                                                                    ggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt     <  2:1739758/93‑1 (MQ=255)                                                                    ggCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGt     <  1:4147958/93‑1 (MQ=255)                                                                     gCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTg    >  1:1354706/1‑93 (MQ=255)                                                                     gCTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTg    >  1:3847057/1‑93 (MQ=255)                                                                      cTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   >  2:3390041/1‑92 (MQ=255)                                                                      cTACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgc   >  1:2566403/1‑92 (MQ=255)                                                                       tACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg  >  1:269785/1‑91 (MQ=255)                                                                       tACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg  <  2:2558569/93‑3 (MQ=255)                                                                       tACAGCCTGAAGGATT‑GTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGGAAAGGCTGCATTATGGGCAGTgcg  <  2:5981713/93‑3 (MQ=255)                                                                            |                                                                                         TCTGGAAGTGATGGCGGGCGGCGTTTC‑AAAGGCCATGCGCTGGAAG‑GCTGGC‑‑‑G‑‑‑‑AACT‑GGCTACAG‑CTGAAGGATT‑G‑ATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGC‑GGGAAA‑GCTGCATTATGGGCAGTGGC  >  NC_000913/4008346‑4008492

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