Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,677,896	G→C	5.8%	N226K (AAC→AAG)	ldcC ←	lysine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,677,896	0	G	C	5.8%	224.8 / 3.7	69	N226K (AAC→AAG)	ldcC	lysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base C (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

ACCCGTCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAT  >  NC_002947/4677749‑4678041                                                                                                                                                    |                                                                                                                                                  aCCCGTCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGt                                                                                                                                                  <  3:851/149‑1 (MQ=255)      tCATGATGCGCGCGTGCACCACCGACTTGGGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTgcgc                                                                                                                                             <  8:77394/149‑1 (MQ=255)         tgatgaTCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGcgcc                                                                                                                                          >  6:292089/1‑149 (MQ=255)             gatCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgat                                                                                                                                       >  4:128534/1‑146 (MQ=255)              atCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgatc                                                                                                                                      <  3:128534/148‑4 (MQ=255)               tCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgatca                                                                                                                                     >  8:247140/1‑144 (MQ=255)               tCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTgc                                                                                                                                               <  1:288737/138‑1 (MQ=255)               tCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTgc                                                                                                                                               >  2:288737/1‑138 (MQ=255)                 gcgTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgatcaca                                                                                                                                   <  7:247140/148‑7 (MQ=255)                 gcgTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGc                                                                                                                                  <  3:148524/149‑1 (MQ=255)                    tGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTc                                                                                                                               <  5:243496/149‑1 (MQ=255)                      caccacCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                             >  8:218683/1‑149 (MQ=255)                           ccGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCg                                                                                                                        <  8:46139/149‑1 (MQ=255)                           ccGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCg                                                                                                                        <  5:120733/149‑1 (MQ=255)                             gACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCGATTAATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                             >  2:156995/1‑142 (MQ=255)                             gACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTAATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                             <  1:156995/142‑1 (MQ=255)                                        gTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCtg                                                                                                           >  2:186623/1‑149 (MQ=255)                                             ggTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTgg                                                                                                      >  2:324120/1‑149 (MQ=255)                                                 caccaggaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                             <  8:270295/122‑1 (MQ=255)                                                 caccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGt                                                                                                     <  8:304872/146‑1 (MQ=255)                                                 caccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGt                                                                                                     >  7:304872/1‑146 (MQ=255)                                                  accagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGag                                                                                                 >  2:211052/1‑149 (MQ=255)                                                   ccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCt                                                                                                            <  2:294929/137‑1 (MQ=255)                                                   ccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCt                                                                                                            >  1:294929/1‑137 (MQ=255)                                                   ccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGag                                                                                                 >  6:224639/1‑148 (MQ=255)                                                    cagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTc                                                                                                    >  3:55150/1‑144 (MQ=255)                                                    cagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTc                                                                                                    <  4:55150/144‑1 (MQ=255)                                                             gtcgtcACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAAcc                                                                                        >  8:20254/1‑147 (MQ=255)                                                                 tcACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTc                                                                                  <  5:224639/149‑1 (MQ=255)                                                                    cGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGc                                                                                                              <  1:99359/118‑1 (MQ=255)                                                                    cGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGc                                                                                                              >  2:99359/1‑118 (MQ=255)                                                                    cGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCgg                                                                               <  1:324120/149‑1 (MQ=255)                                                                         gACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACgg                                                                          <  7:20254/149‑1 (MQ=255)                                                                             ccgatctGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGaaa                                                                       <  1:59693/142‑1 (MQ=255)                                                                                   tGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAg                                                                <  4:169276/149‑1 (MQ=255)                                                                                     ccAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGAGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGaaa                                                                       >  2:59693/1‑140 (MQ=255)                                                                                      cAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTc                                                             >  7:85797/1‑149 (MQ=255)                                                                                       aGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGAGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTcc                                                            >  7:306949/1‑149 (MQ=255)                                                                                          cGATCATGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGACTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAg                                                         >  6:324310/1‑149 (MQ=255)                                                                                               ttgttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAg                                                    >  1:84260/1‑149 (MQ=255)                                                                                                 gttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAgg                                                   <  1:186623/148‑1 (MQ=255)                                                                                                  ttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACgg                                                                    <  7:203253/130‑1 (MQ=255)                                                                                                  ttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACgg                                                                    >  8:203253/1‑130 (MQ=255)                                                                                                  ttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGgt                                                 >  1:26205/1‑149 (MQ=255)                                                                                                    ggCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGc                                                                              >  2:166461/1‑118 (MQ=255)                                                                                                       cGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGtttt                                            >  7:18561/1‑149 (MQ=255)                                                                                                           gAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                         <  2:84260/148‑1 (MQ=255)                                                                                                                  cATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGGCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAAc                                 <  7:218683/149‑1 (MQ=255)                                                                                                                   aTTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACt                                >  4:22347/1‑149 (MQ=255)                                                                                                                       atgaCGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGt                             >  5:315599/1‑148 (MQ=255)                                                                                                                       atgaCGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGt                             <  3:22347/148‑1 (MQ=255)                                                                                                                         gaCGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCAc                                                                            <  2:169939/99‑1 (MQ=255)                                                                                                                            gAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAgg                       >  1:160099/1‑149 (MQ=255)                                                                                                                            gAAAAAGGTGTGATCGGCACCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAg                        >  5:279098/1‑148 (MQ=255)                                                                                                                               aaaGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCtt                    <  8:18561/149‑1 (MQ=255)                                                                                                                                aaGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTg                   >  1:187810/1‑149 (MQ=255)                                                                                                                                  ggTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGc                 <  2:26205/149‑1 (MQ=255)                                                                                                                                    tgtgATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGcc                >  5:174818/1‑148 (MQ=255)                                                                                                                                            gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt       <  6:355555/149‑1 (MQ=255)                                                                                                                                            gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt       >  5:355555/1‑149 (MQ=255)                                                                                                                                            gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt       >  2:217090/1‑149 (MQ=255)                                                                                                                                                cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCGTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTa   <  5:324310/149‑1 (MQ=255)                                                                                                                                                cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                         >  4:198416/1‑111 (MQ=255)                                                                                                                                                cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                         <  3:198416/111‑1 (MQ=255)                                                                                                                                                cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTa   <  1:217090/149‑1 (MQ=255)                                                                                                                                                 gAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAt  <  6:315599/149‑1 (MQ=255)                                                                                                                                                 gAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAt  <  8:85797/149‑1 (MQ=255)                                                                                                                                                  aaGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCggg         <  7:19110/141‑1 (MQ=255)                                                                                                                                                  aaGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCggg         >  8:19110/1‑141 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ACCCGTCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAT  >  NC_002947/4677749‑4678041

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: ‑