Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,116,500	T→G	7.1%	L29R (CTG→CGG)	PP_4502 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,116,500	0	T	G	7.1%	259.2 / 9.9	84	L29R (CTG→CGG)	PP_4502	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (40/38);  new base G (2/4);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

CCTGGGCATGACCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTCGAATCG  >  NC_002947/5116353‑5116648                                                                                                                                                    |                                                                                                                                                      ccTGGGCATGACCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTg                                                                                                                                                      <  5:83594/149‑1 (MQ=255)  cTGGGCATGACCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTgg                                                                                                                                                     <  2:141877/149‑1 (MQ=255)   tGGGCATGACCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGgc                                                                                                                                                    >  3:315922/1‑149 (MQ=255)       cATGACCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGtccgcc                                                                                                                                                <  3:131394/149‑1 (MQ=255)        aTGACCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGccgccc                                                                                                                                               <  6:305184/149‑1 (MQ=255)           aCCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTg                                                                                                                                                      >  6:349123/1‑139 (MQ=255)           aCCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTg                                                                                                                                                      <  5:349123/139‑1 (MQ=255)            ccGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCGGgcgggcggc                                                                                                                                             >  6:148700/1‑140 (MQ=255)                  tCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCGGGCGGGCGGCCAGgcgc                                                                                                                                      <  4:334747/148‑1 (MQ=255)                     tCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACAt                                                                                                                                  >  4:15021/1‑149 (MQ=255)                          gtggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGaa                                                                                                                              <  4:92552/148‑1 (MQ=255)                          gtggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAAc                                                                                                                             <  3:97719/149‑1 (MQ=255)                          gtggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCGGGCGGGCGGCCAGgcgcttttcca                                                                                                                               <  5:148700/147‑8 (MQ=255)                          gtggCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCGGGCGGGCGGCCAGgcgcttttcca                                                                                                                               <  2:316522/147‑8 (MQ=255)                              cgcAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCt                                                                                                                         <  5:359029/149‑1 (MQ=255)                                         tACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTg                                                                                                                        >  2:281842/1‑139 (MQ=255)                                         tACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTg                                                                                                                        <  1:281842/139‑1 (MQ=255)                                         tACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATCGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTc                                                                                                               <  2:23064/148‑1 (MQ=255)                                            cTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaagg                                                                                                           <  2:8731/149‑1 (MQ=255)                                              gTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATc                                                                                                                                 <  2:408973/125‑1 (MQ=255)                                              gTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATc                                                                                                                                 >  1:408973/1‑125 (MQ=255)                                                   gcCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggccgaagg                                                                                                    <  5:399889/149‑1 (MQ=255)                                                    ccACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCaaggccgaaggc                                                                                                   >  3:364485/1‑149 (MQ=255)                                                       ccTGCACGATGGCCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCcgc                                                                                                <  5:409842/149‑1 (MQ=255)                                                         tGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCgt                                                                                              >  8:128197/1‑149 (MQ=255)                                                          gCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGCCCGAAGGCCGCgtg                                                                                             >  4:242309/1‑149 (MQ=255)                                                                  aCCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggc                                                                                     >  2:346334/1‑149 (MQ=255)                                                                  aCCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTgg                                                                                      <  5:30399/148‑1 (MQ=255)                                                                   cccTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggcc                                                                                    >  3:396426/1‑149 (MQ=255)                                                                   cccTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggcc                                                                                    >  7:143795/1‑149 (MQ=255)                                                                    ccTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggccg                                                                                   >  4:338266/1‑149 (MQ=255)                                                                    ccTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggccg                                                                                   >  7:363291/1‑149 (MQ=255)                                                                        cTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggccggcc                                                                                <  8:319146/148‑1 (MQ=255)                                                                        cTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTggccggc                                                                                 <  4:364485/147‑1 (MQ=255)                                                                          aCCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCg                                                                              <  8:143795/148‑1 (MQ=255)                                                                          aCCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGc                                                                             >  6:202237/1‑149 (MQ=255)                                                                               cAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCg                                                                        >  3:48791/1‑149 (MQ=255)                                                                                aaaCCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCgg                                                                       <  5:265733/149‑1 (MQ=255)                                                                                    ccAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTgg                                                                                      <  6:338976/130‑1 (MQ=255)                                                                                    ccAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTgg                                                                                      >  5:338976/1‑130 (MQ=255)                                                                                    ccAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑t                                                                  >  5:34579/1‑149 (MQ=255)                                                                                     cgaGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGTGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑Tc                                                                 <  7:383507/147‑1 (MQ=255)                                                                                      aaGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTg                                                                                       <  4:268854/127‑1 (MQ=255)                                                                                      aaGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTg                                                                                       >  3:268854/1‑127 (MQ=255)                                                                                       aGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGa                                                               <  5:202237/149‑1 (MQ=255)                                                                                         gATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGa                                                               >  1:406170/1‑147 (MQ=255)                                                                                              aCATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAAc                                                        >  4:296626/1‑149 (MQ=255)                                                                                                  gCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGAcc                                                    >  2:247972/1‑149 (MQ=255)                                                                                                    tctACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAACGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACCGCAACGACCCg                                                  >  4:79006/1‑149 (MQ=255)                                                                                                       aCGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCCAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGc                                               >  7:275944/1‑149 (MQ=255)                                                                                                        cGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCggg                                                <  3:338266/147‑1 (MQ=255)                                                                                                        cGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGc                                               <  8:363291/148‑1 (MQ=255)                                                                                                         gCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCg                                              >  2:19946/1‑148 (MQ=255)                                                                                                           cATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAg                                           >  7:351303/1‑149 (MQ=255)                                                                                                            aTTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGc                                          >  2:219/1‑149 (MQ=255)                                                                                                            aTTATCGCCAGCGACGCCGCAAACTCCCTGGAACCCCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCCGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGc                                          <  2:341688/149‑1 (MQ=255)                                                                                                               aTCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCggg                                       >  3:111605/1‑149 (MQ=255)                                                                                                               aTCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCgg                                        <  3:296626/148‑1 (MQ=255)                                                                                                               aTCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCgg                                        <  6:192462/148‑1 (MQ=255)                                                                                                                tCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCACTCGACAGCAACGACCCGGGCGAAGCgg                                        <  1:19946/148‑1 (MQ=255)                                                                                                                tCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCggg                                       <  3:15021/148‑1 (MQ=255)                                                                                                                tCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGt                                      <  4:315922/149‑1 (MQ=255)                                                                                                                   ccAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCTGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTc                                   >  5:244778/1‑149 (MQ=255)                                                                                                                      gCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGc                                <  8:351303/149‑1 (MQ=255)                                                                                                                       cGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCGGAAGGCCGGGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCg                               <  3:242309/149‑1 (MQ=255)                                                                                                                       cGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCg                               <  4:396426/149‑1 (MQ=255)                                                                                                                            tCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGACGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGc                          >  5:163593/1‑149 (MQ=255)                                                                                                                             cGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAAc                                                                                                                             <  3:273116/50‑1 (MQ=255)                                                                                                                             cGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAAc                                                                                                                             >  4:273116/1‑50 (MQ=255)                                                                                                                              gCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGcc                         <  7:128197/148‑1 (MQ=255)                                                                                                                              gCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCt                        >  1:292933/1‑149 (MQ=255)                                                                                                                                aaaCTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGa                      >  7:36174/1‑149 (MQ=255)                                                                                                                                  aCTCCCTGGAAAAGCGGCTGGGCGCCCGGCCGGCGGAAATGGTACGCCTGCAGCAGGTCAATGCCGTAGGCCGGGGGGTGTTGGCCGGCCCGGACCCGGCCA‑TCGAAATCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATc                    <  4:48791/149‑1 (MQ=255)                                                                                                                                      ccTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTg                >  6:426467/1‑149 (MQ=255)                                                                                                                                      ccTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTg                >  7:198613/1‑149 (MQ=255)                                                                                                                                          gAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACCGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCg             >  2:300169/1‑148 (MQ=255)                                                                                                                                          gAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCg             >  8:231340/1‑148 (MQ=255)                                                                                                                                            aaaGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAAt          >  7:128546/1‑149 (MQ=255)                                                                                                                                             aaGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAAtt         >  3:86120/1‑149 (MQ=255)                                                                                                                                             aaGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAAtt         <  4:111605/149‑1 (MQ=255)                                                                                                                                               gcgcCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTc        >  3:209012/1‑148 (MQ=255)                                                                                                                                               gcgcCGGGCGGGCGGCCAGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTg                >  1:148976/1‑140 (MQ=255)                                                                                                                                               gcgcCGGGCGGGCGGCCAGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTg                <  2:148976/140‑1 (MQ=255)                                                                                                                                                   cTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTCGAATc   >  6:425341/1‑149 (MQ=255)                                                                                                                                                    tGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTCGAATCg  >  1:409359/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCTGGGCATGACCGTGATCTTCCTGGTGGCGCAAGGCGTGTACCTGTCGCGCCACCTGCACGATGACCCTTCTACCTCCAAACCCAAGGATTGACATGCTCTACGCCATTATCGCCAGCGACGTCGCAAACTCCCTGGAAAAGCGCCTGGCCGCCCGCCCGGCGCACATCGAACGCCTGCAGCAGCTCAAGGCCGAAGGCCGCGTGGTGTTGGCCGGCCCGCACCCGGCCA‑TCGACAGCAACGACCCGGGCGAAGCGGGTTTCAGCGGTAGCCTGATCGTTGCCGAATTCGAATCG  >  NC_002947/5116353‑5116648

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