Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,509,703	T→C	N120S (AAC→AGC)	serS ←	serine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,509,703	0	T	C	100.0%	245.7 / NA	69	N120S (AAC→AGC)	serS	serine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/33);  total (36/33)

GCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGTTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAG  >  NC_002947/4509566‑4509844                                                                                                                                          |                                                                                                                                              gcagcaCAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCa                                                                                                                                    >  4:125895/1‑149 (MQ=255) gcagcaCAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCa                                                                                                                                    <  3:303940/149‑1 (MQ=255)   agcacagcaAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACt                                                                                                                                  >  2:22866/1‑149 (MQ=255)        agcaAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGgcgc                                                                                                                             <  6:311061/149‑1 (MQ=255)        agcaAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGgcgc                                                                                                                             <  2:119720/149‑1 (MQ=255)          caAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTgtc                                                                                                                                           >  5:357904/1‑133 (MQ=255)          caAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTgtc                                                                                                                                           <  6:357904/133‑1 (MQ=255)            aaaCGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGAc                                                                                                                         >  7:125676/1‑149 (MQ=255)                 ggCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCa                                                                                                                    >  4:99511/1‑149 (MQ=255)                  gCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCAc                                                                                                                   >  1:126330/1‑149 (MQ=255)                      cGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTg                                                                                                               <  2:327818/149‑1 (MQ=255)                       ggACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTgg                                                                                                              <  7:16415/149‑1 (MQ=255)                            gCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTcc                                                                                                         >  5:51796/1‑149 (MQ=255)                             cTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCg                                                                                                        >  1:187909/1‑149 (MQ=255)                             cTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCg                                                                                                        <  5:111345/149‑1 (MQ=255)                                     ccGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTg                                                                                                <  1:22866/149‑1 (MQ=255)                                        cTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTggg                                                                                              >  5:177590/1‑148 (MQ=255)                                           cGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGaa                                                                                            <  8:133353/147‑1 (MQ=255)                                            gAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATgg                                                                                         >  5:20046/1‑149 (MQ=255)                                            gAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATg                                                                                          <  2:176602/148‑1 (MQ=255)                                            gAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCACCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGACCCCAACGGCACCCAGGCGACCGGCAGCTTGCGaa                                                                                            >  7:147311/1‑146 (MQ=255)                                              aGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTc                                                                                       <  8:125676/149‑1 (MQ=255)                                                tCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTcag                                                                                     <  3:125895/149‑1 (MQ=255)                                                tCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTcag                                                                                     <  2:126330/149‑1 (MQ=255)                                                        aCCGCTGATTTCACCGAGGGCCACATGGTCATTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGTCGCACTTCGACGCTGTAGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGc                                                                             <  6:63303/149‑1 (MQ=255)                                                         ccGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGcc                                                                            <  6:62828/149‑1 (MQ=255)                                                           gCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCAt                                                                          >  4:265156/1‑149 (MQ=255)                                                           gCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCAt                                                                          <  6:51796/149‑1 (MQ=255)                                                            cTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATc                                                                         >  8:291180/1‑149 (MQ=255)                                                            cTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATc                                                                         >  3:96706/1‑149 (MQ=255)                                                                tttCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGc                                                                     >  4:90722/1‑149 (MQ=255)                                                                    aCCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCgg                                                                 >  5:174280/1‑149 (MQ=255)                                                                      cGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTc                                                                                                                                     <  2:245049/79‑1 (MQ=255)                                                                      cGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTc                                                                                                                                     >  1:245049/1‑79 (MQ=255)                                                                            ccACGTGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATg                                                         <  6:20046/149‑1 (MQ=255)                                                                                aTGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATc                                                                         >  7:185574/1‑129 (MQ=255)                                                                                aTGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATc                                                                         <  8:185574/129‑1 (MQ=255)                                                                                aTGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGt                                                     >  4:131022/1‑149 (MQ=255)                                                                                   gTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATg                                                                                          >  4:48718/1‑109 (MQ=255)                                                                                   gTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATg                                                                                          <  3:48718/109‑1 (MQ=255)                                                                                    tCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCt                                                              <  8:316887/136‑1 (MQ=255)                                                                                    tCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCt                                                              >  7:316887/1‑136 (MQ=255)                                                                                      cTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGc                                                >  4:196414/1‑148 (MQ=255)                                                                                      cTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCt                                               >  2:38204/1‑149 (MQ=255)                                                                                       ttGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTc                                                                   <  1:193042/128‑1 (MQ=255)                                                                                       ttGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTc                                                                   >  2:193042/1‑128 (MQ=255)                                                                                       ttGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGaa                                                           >  4:24338/1‑136 (MQ=255)                                                                                       ttGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGaa                                                           <  3:24338/136‑1 (MQ=255)                                                                                       ttGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCCTCTTCACTGGCGCCCACCGGCACGCAGGCGTCCGGCAGGTTGGGAAAGGTCAGCAGGATGCCATCCAGCTCGCACTGAATGGCGTCCAGCt                                               >  4:287775/1‑148 (MQ=255)                                                                                              cGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTg                                       >  6:205901/1‑149 (MQ=255)                                                                                              cGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTg                                       >  6:210949/1‑149 (MQ=255)                                                                                               gAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTg                                       <  1:354557/148‑1 (MQ=255)                                                                                                 aGTCGAACGCGTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCg                                    <  6:177590/149‑1 (MQ=255)                                                                                                  gTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCgg                                   <  3:131022/149‑1 (MQ=255)                                                                                                     gAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAg                                >  2:311875/1‑149 (MQ=255)                                                                                                     gAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAg                                <  3:5532/149‑1 (MQ=255)                                                                                                     gAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCCTACGGCAGGTTGGGAATGGTCAGCAGGATCCCATCCAGATCGGCCTGAATGGCGTCCAGCTCATCCTTgcaggcag                                >  8:327566/1‑149 (MQ=255)                                                                                                        cGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCa                             <  2:122826/149‑1 (MQ=255)                                                                                                           cTTCGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGtt                          >  1:270378/1‑149 (MQ=255)                                                                                                              cGGCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTc                         <  3:196414/147‑1 (MQ=255)                                                                                                                gCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTg                     <  3:265156/149‑1 (MQ=255)                                                                                                                gCGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTg                     >  8:237878/1‑149 (MQ=255)                                                                                                                 cGTGCCCCAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTgg                    <  3:99511/149‑1 (MQ=255)                                                                                                                     ccccAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGaa                                                                                            >  8:226929/1‑73 (MQ=255)                                                                                                                        cAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGTCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATgcg             >  3:324039/1‑149 (MQ=255)                                                                                                                        cAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTcag                                                                                     <  3:332435/77‑1 (MQ=255)                                                                                                                        cAACGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTcag                                                                                     >  4:332435/1‑77 (MQ=255)                                                                                                                           cGGCGCACTTCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTc          <  3:90722/149‑1 (MQ=255)                                                                                                                                   ttCGACGCTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATACAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAg  >  2:341651/1‑149 (MQ=255)                                                                                                                                          |                                                                                                                                              GCAGCACAGCAAAACGGGCACCGGACAGCTTGGCAGCCGCTTCGAAGTCCAGGCCACCGCTGATTTCACCGAGGGCCACATGGTCCTTGATTTCGAAGTCGAACGCCTTCGGCGTGCCCCAACGGCGCACTTCGACGTTGTCGTCTTCACTGGCGCCGACCGGCACGCTGGCGTCCGGCAGGTTGGGAATGGTCAGCAGGATGCCATCCAGCTCGGCCTGAATGGCGTCCAGCTCAGCCTTGCCGGCAGCCAGTTCGTTGGCCATGCGCTCGACGTCAG  >  NC_002947/4509566‑4509844

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