Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,508,373	A→C	5.2%	L105R (CTT→CGT)	PP_2202 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,508,373	0	A	C	5.2%	291.4 / 4.6	91	L105R (CTT→CGT)	PP_2202	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (44/42);  new base C (2/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACCACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCAGCAGC  >  NC_002947/2508227‑2508519                                                                                                                                                    |                                                                                                                                                    aCCACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGc                                                                                                                                                   <  8:213089/149‑1 (MQ=255)  ccACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCg                                                                                                                                                  >  7:1226/1‑149 (MQ=255)  ccACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCg                                                                                                                                                  >  7:6389/1‑149 (MQ=255)   cACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCgg                                                                                                                                                 >  1:205527/1‑149 (MQ=255)     cgcgAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGt                                                                                                                                               <  4:338801/149‑1 (MQ=255)        gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGc                                                                                                                                            >  5:122123/1‑149 (MQ=255)        gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAAGCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGc                                                                                                                                            <  2:58835/149‑1 (MQ=255)         aatccTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCAGGCATTATCCCGGCGATCAGCGTAAGCGGGTCGcc                                                                                                                                           <  2:313830/146‑1 (MQ=255)             cTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCACCAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCgg                                                                                                                                                 >  8:135940/1‑139 (MQ=255)             cTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCgg                                                                                                                                                 <  7:135940/139‑1 (MQ=255)              tGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCa                                                                                                                                      >  5:85391/1‑149 (MQ=255)                      gtTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATc                                                                                                                              <  3:188830/149‑1 (MQ=255)                      gtTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATc                                                                                                                              >  8:27683/1‑149 (MQ=255)                       ttACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGAGCCGCGTAAGCGGGTCGCCGATCACCGGCATc                                                                                                                              <  1:27667/148‑1 (MQ=255)                          ccAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAg                                                                                                                           <  6:178099/148‑1 (MQ=255)                             ggCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGa                                                                                                                       <  6:85391/149‑1 (MQ=255)                              gCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGag                                                                                                                      <  7:27683/149‑1 (MQ=255)                                gTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCggga                                                                                                                                       >  1:33924/1‑127 (MQ=255)                                gTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCggga                                                                                                                                       <  2:33924/130‑4 (MQ=255)                                    aaCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGagcagcag                                                                                                                >  6:232867/1‑149 (MQ=255)                                      ccAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTGTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGAATAAACCCGCCGATCAGCGTAAGCGGGTCGCCGATAACCGGCATCCAGCTGagcagcagc                                                                                                               >  1:248553/1‑149 (MQ=255)                                       cAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TTGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGc                                                                                                                          >  8:212594/1‑136 (MQ=255)                                       cAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGc                                                                                                                          <  7:212594/136‑1 (MQ=255)                                             aGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTg                                                                                                       <  8:6389/149‑1 (MQ=255)                                             aGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTg                                                                                                       <  8:1226/149‑1 (MQ=255)                                                          cAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTgg                                                                                         <  2:205527/149‑1 (MQ=255)                                                          cAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTgg                                                                                         >  8:183993/1‑149 (MQ=255)                                                           aaGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTgg                                                                                         >  6:148215/1‑148 (MQ=255)                                                                ccGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAgcgc                                                                                   >  1:156855/1‑149 (MQ=255)                                                                   cATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCgctgc                                                                                >  3:149145/1‑149 (MQ=255)                                                                     tAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGcc‑cc                                                                                                  >  7:248326/1‑129 (MQ=255)                                                                     tAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGcc‑cc                                                                                                  <  8:248326/129‑1 (MQ=255)                                                                        cGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTgcg                                                                            >  3:155692/1‑148 (MQ=255)                                                                          gCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCAc                                                                         >  5:122843/1‑149 (MQ=255)                                                                              gCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACgcgc                                                                     <  6:122843/149‑1 (MQ=255)                                                                              gCCT‑TGGACAGGGTCACCAGCAACAGGAAGCGCCATAATGGCTCACTCATTATCCCGGCTATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACgcgc                                                                     <  4:149145/149‑1 (MQ=255)                                                                                     gCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCg                                                                                                                                   <  4:218183/82‑1 (MQ=255)                                                                                     gCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCg                                                                                                                                   >  3:218183/1‑82 (MQ=255)                                                                                      ccAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTgg                                                              <  5:232867/149‑1 (MQ=255)                                                                                       cAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTggg                                                             >  8:991/1‑149 (MQ=255)                                                                                        aGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTgg                                                                                         <  2:324714/120‑1 (MQ=255)                                                                                        aGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTgg                                                                                         >  1:324714/1‑120 (MQ=255)                                                                                            tCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAgcgc                                                                                   <  4:196764/122‑1 (MQ=255)                                                                                              aCCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCt                                                      >  7:216747/1‑149 (MQ=255)                                                                                               ccAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCAcgc                                                                       <  7:111446/131‑1 (MQ=255)                                                                                               ccAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCAcgc                                                                       >  8:111446/1‑131 (MQ=255)                                                                                                  gCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGaaa                                                  <  8:147771/149‑1 (MQ=255)                                                                                                  gCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGaaa                                                  <  2:106624/149‑1 (MQ=255)                                                                                                    aaCAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAgg                                                >  4:58042/1‑149 (MQ=255)                                                                                                     aCAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAgg                                                >  2:275232/1‑148 (MQ=255)                                                                                                      cAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAggg                                               >  4:100335/1‑148 (MQ=255)                                                                                                        ggAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTggg                                                             <  4:153462/132‑1 (MQ=255)                                                                                                        ggAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTggg                                                             >  3:153462/1‑132 (MQ=255)                                                                                                        ggAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAAc                                            <  4:155692/149‑1 (MQ=255)                                                                                                        ggAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCACCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAAc                                            >  2:28183/1‑149 (MQ=255)                                                                                                         gAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCgctgct                                                                               <  1:64933/113‑1 (MQ=255)                                                                                                         gAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCgctgct                                                                               >  2:64933/1‑113 (MQ=255)                                                                                                          aaGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCACCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACCCGCCAGTTGGGCAGCGCTGAAAGGGAACCa                                          <  1:28183/149‑1 (MQ=255)                                                                                                             cgcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgcg                                       <  7:183993/149‑1 (MQ=255)                                                                                                              gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgcgc                                      <  2:156855/149‑1 (MQ=255)                                                                                                              gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgcgc                                      <  3:100335/149‑1 (MQ=255)                                                                                                              gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgcgc                                      >  4:257347/1‑149 (MQ=255)                                                                                                                 aGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTg                                   >  7:148504/1‑149 (MQ=255)                                                                                                                   aaTGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTc                                 >  3:260352/1‑149 (MQ=255)                                                                                                                   aaTGGCTCACGCATTATCACGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑ACCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTc                                 >  5:191289/1‑149 (MQ=255)                                                                                                                     tGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTcgc                               <  6:122123/149‑1 (MQ=255)                                                                                                                       gCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgcgc                                      >  7:310798/1‑140 (MQ=255)                                                                                                                       gCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgcgc                                      <  8:310798/140‑1 (MQ=255)                                                                                                                       gCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTcgcg                              >  4:180919/1‑148 (MQ=255)                                                                                                                        cTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTAGGGCAGCGCTGAAAGGGAACCAGCGCTTGTcgcg                              >  2:78795/1‑147 (MQ=255)                                                                                                                         tCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAg                           >  3:274436/1‑149 (MQ=255)                                                                                                                               aTTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAAc                                            >  1:22637/1‑126 (MQ=255)                                                                                                                               aTTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAAc                                            <  2:22637/126‑1 (MQ=255)                                                                                                                               aTTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTc                     >  2:71874/1‑149 (MQ=255)                                                                                                                                 tATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGa                   >  7:319929/1‑149 (MQ=255)                                                                                                                                 tATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGa                   >  5:133101/1‑149 (MQ=255)                                                                                                                                  tgatcGGCGACCCGCTTACGCTGATCGCCGATCAACGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGAt                  >  4:252655/5‑149 (MQ=255)                                                                                                                                   tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAAc                                            <  6:76114/122‑1 (MQ=255)                                                                                                                                   tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGAt                  <  3:180919/148‑1 (MQ=255)                                                                                                                                   tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGAt                  <  4:274436/148‑1 (MQ=255)                                                                                                                                   tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATg                 >  5:159454/1‑149 (MQ=255)                                                                                                                                   tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATg                 <  5:148215/149‑1 (MQ=255)                                                                                                                                    cccGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATgg                >  2:207549/1‑149 (MQ=255)                                                                                                                                      cGGTGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcc              >  8:299408/1‑149 (MQ=255)                                                                                                                                      cGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcc              <  1:71874/149‑1 (MQ=255)                                                                                                                                      cGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcc              >  4:99882/1‑149 (MQ=255)                                                                                                                                       ggCTATCAGCTTACGGGGGGCGCCGATCACCGGCACCCCGCTGCGCAGCAGCGACCTTGGCC‑CCCTGCGTTGGTAGCGCGGCTGCGCACGCGCCAGTTGGGCCGCGCTGAAAGGCACCCAGCGCTTGGGGCGCAGGTGCTCGATggccc             <  6:159454/149‑1 (MQ=255)                                                                                                                                         cGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCGGGTGCTCGATggcccgg           <  1:275232/149‑1 (MQ=255)                                                                                                                                         cGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcccgg           <  3:99882/149‑1 (MQ=255)                                                                                                                                          gATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTCCTCGATggcccggc          <  7:299408/149‑1 (MQ=255)                                                                                                                                               gCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCcagc     <  2:241507/149‑1 (MQ=255)                                                                                                                                                    aGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCGCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCagcagc  >  1:257356/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACCACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCT‑TGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCC‑CCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCAGCAGC  >  NC_002947/2508227‑2508519

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