Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,677,899	G→C	6.1%	R225R (CGC→CGG) ‡	ldcC ←	lysine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,677,899	0	G	C	6.1%	207.2 / 3.7	65	R225R (CGC→CGG) ‡	ldcC	lysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base C (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

TCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAT  >  NC_002947/4677754‑4678041                                                                                                                                                  |                                                                                                                                               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                                                                                                     >  7:304872/1‑146 (MQ=255)                                            caccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGt                                                                                                     <  8:304872/146‑1 (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                                                                                                              >  2:99359/1‑118 (MQ=255)                                                               cGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGc                                                                                                              <  1:99359/118‑1 (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)                                                                                                                                                  |                                                                                                                                               TCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAT  >  NC_002947/4677754‑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: ‑