Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,162,202	T→A	5.7%	W324R (TGG→AGG) ‡	PP_4544 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,162,202	0	T	A	5.7%	215.8 / 2.5	70	W324R (TGG→AGG) ‡	PP_4544	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base A (2/2);  total (35/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 = 6.06e-01

CATCCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162054‑5162347                                                                                                                                                     |                                                                                                                                                  cATCCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCct                                                                                                                                                   >  2:21381/1‑149 (MQ=255)  aTCCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctg                                                                                                                                                  >  3:9679/1‑149 (MQ=255)   tCCTGCAGTCGTTCGGCAACAACACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACCTGAAGAACACCGACAGCACAGACGGCGAAGATCGGGCctgg                                                                                                                                                 >  1:202368/1‑149 (MQ=255)     cTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggct                                                                                                                                               <  1:17090/149‑1 (MQ=255)        cAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCctggctgg                                                                                                                                             <  8:46396/148‑1 (MQ=255)                cGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTg                                                                                                                                     <  4:9679/148‑1 (MQ=255)                cGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCAGCCAGCCAGGCCTg                                                                                                                                     <  7:148153/148‑1 (MQ=255)                  gCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaa                                                                                                                                   >  8:36401/1‑148 (MQ=255)                   cAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaaaa                                                                                                                                 >  1:243121/1‑149 (MQ=255)                   cAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCAGCCAGCCAGGCCTGatct                                                                                                                                 <  7:145391/149‑4 (MQ=255)                    aaCACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAAc                                                                                                                                <  2:243121/149‑1 (MQ=255)                    aaCACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAAc                                                                                                                                <  7:30639/149‑1 (MQ=255)                       accaccAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCaa                                                                                                                             <  1:205260/149‑1 (MQ=255)                            caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATgg                                                                                                                        >  8:75180/1‑149 (MQ=255)                            caACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGATGGCCTGAAAACCAAAATgg                                                                                                                        >  4:250116/1‑149 (MQ=255)                                 aGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcc                                                                                                                   <  6:189605/149‑1 (MQ=255)                                  gTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcct                                                                                                                  >  3:25434/1‑149 (MQ=255)                                      aaTCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggag                                                                                                              >  3:142410/1‑149 (MQ=255)                                        tCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcc                                                                                                            >  3:120345/1‑149 (MQ=255)                                            gtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctgga                                                                                                        <  8:69898/149‑1 (MQ=255)                                            gtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                          <  8:92269/147‑1 (MQ=255)                                            gtggAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTgcctggagcctg                                                                                                          <  7:36401/147‑1 (MQ=255)                                               gAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACt                                                                                                      <  7:304392/148‑1 (MQ=255)                                                        ggCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcg                                                                                            >  8:221252/1‑149 (MQ=255)                                                         gCCAGATCGACTTCACTGGCTTGTCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCAACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACACCTACcgc                                                                                           >  6:27248/1‑149 (MQ=255)                                                         gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaaa                                                                                                                                  >  2:70200/1‑110 (MQ=255)                                                         gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGaaa                                                                                                                                  <  1:70200/110‑1 (MQ=255)                                                         gCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACcgc                                                                                           >  7:150356/1‑149 (MQ=255)                                                               tCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAtgtg                                                                                     <  7:75180/149‑1 (MQ=255)                                                                 gACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGAtgtg                                                                                     >  6:169376/1‑147 (MQ=255)                                                                 gACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCa                                                                                   <  4:25434/149‑1 (MQ=255)                                                                      cACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGt                                                                               <  5:169376/148‑1 (MQ=255)                                                                       aCTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTaa                                                                             >  8:247288/1‑149 (MQ=255)                                                                        cTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACAGACGGCGAAGATCACGCCTCGAAGGCAGGCCTCAAAACAAAAATGCCTGCCTGGATCCTGCACAACAACTACCGCCATGTGCAGCGAAAc                                                                            >  6:209656/1‑149 (MQ=255)                                                                           gCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAAcgcg                                                                         >  8:290774/1‑149 (MQ=255)                                                                             ttGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggt                                                                       >  2:259780/1‑149 (MQ=255)                                                                               ggCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtg                                                                      <  1:21381/148‑1 (MQ=255)                                                                               ggCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtg                                                                      <  3:250116/148‑1 (MQ=255)                                                                                  cATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGt                                                                               <  4:268198/136‑1 (MQ=255)                                                                                  cATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGt                                                                               >  3:268198/1‑136 (MQ=255)                                                                                    ttCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGc                                                                >  2:211833/1‑149 (MQ=255)                                                                                    ttCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCAACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGc                                                                <  5:27248/149‑1 (MQ=255)                                                                                      ccGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGcc                                                               <  4:90539/148‑1 (MQ=255)                                                                                            tCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCa                                                         <  1:259780/148‑1 (MQ=255)                                                                                                 cTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGCtt                                                   >  1:98731/1‑149 (MQ=255)                                                                                                      gCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCATCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGAc                                              >  3:194343/1‑149 (MQ=255)                                                                                                      gCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGAc                                              <  1:211833/149‑1 (MQ=255)                                                                                                       ccAGTACGTGAAGAACACCGACAGCACCGAGGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACt                                             >  1:27592/1‑149 (MQ=255)                                                                                                          gTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCg                                          <  1:123267/149‑1 (MQ=255)                                                                                                          gTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCg                                          <  7:247288/149‑1 (MQ=255)                                                                                                           tACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCgg                                         >  2:115953/1‑149 (MQ=255)                                                                                                            aCGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCgg                                         >  7:128260/1‑148 (MQ=255)                                                                                                               tGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAAc                                     >  5:114469/1‑149 (MQ=255)                                                                                                                gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtg                                                                      >  8:34143/1‑115 (MQ=255)                                                                                                                gaagaaCACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCggtg                                                                      <  7:34143/115‑1 (MQ=255)                                                                                                                      cacCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGAGGTGGTCAGCCTGTTCACCGATTCGTACTTCGGCAACGGCTaca                              <  2:98731/149‑1 (MQ=255)                                                                                                                        ccGACAGCAACGAAGGCGAAGATCAGGCCTGGCTGGCTGTCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGAAGCGTAACTCGGTGGTCAGCCTGTTAACATATTCGGACTTCGTCAAAGGCTACAcc                            >  4:92545/1‑149 (MQ=255)                                                                                                                           aCAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACg                                    >  4:157601/1‑138 (MQ=255)                                                                                                                           aCAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACg                                    <  3:157601/138‑1 (MQ=255)                                                                                                                               cACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCgg                                         >  1:308923/1‑129 (MQ=255)                                                                                                                               cACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCgg                                         <  2:308923/129‑1 (MQ=255)                                                                                                                                  cGACGGCGAAGATCAGGGCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTcgcgcg                  <  2:102812/149‑1 (MQ=255)                                                                                                                                  cGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTcgcgcg                  <  8:128260/149‑1 (MQ=255)                                                                                                                                       gCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTcgcgcg                  <  5:169399/144‑1 (MQ=255)                                                                                                                                       gCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTcgcgcg                  >  6:169399/1‑144 (MQ=255)                                                                                                                                         gAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTcac             >  1:9901/1‑147 (MQ=255)                                                                                                                                         gAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTcac             <  1:115953/147‑1 (MQ=255)                                                                                                                                         gAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACaa           <  7:221252/149‑1 (MQ=255)                                                                                                                                         gAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACaa           >  1:182808/1‑149 (MQ=255)                                                                                                                                           ttttCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGt         >  1:68310/4‑149 (MQ=255)                                                                                                                                           aGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGt         <  8:150356/149‑1 (MQ=255)                                                                                                                                              tCAGGCCAGCCAGCCAGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTACCGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTc       >  8:117839/1‑148 (MQ=255)                                                                                                                                               cAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGt         >  3:178882/1‑145 (MQ=255)                                                                                                                                               cAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGt         <  4:178882/145‑1 (MQ=255)                                                                                                                                                  gCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGt  >  8:172310/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CATCCTGCAGTCGTTCGGCAACACCACCAACGAGTTCAATCTGGTGGAAGGCTTTGGCCAGATCGACTTCACTGGCTTGGCCATTCCGTTGTCGGCCTATGGCCAGTACGTGAAGAACACCGACAGCACCGACGGCGAAGATCAGGCCTGGCTGGCTGGCCTGAAAACCAAAATGGGTGCCTGGAGCCTGGACTACAACTACCGCGATGTGCAGCGTAACGCGGTGGTCAGCCTGTTCACCGATTCGGACTTCGGCAACGGCTACACCGGTTCGCGCGGTCACAAGTTCAAGGT  >  NC_002947/5162054‑5162347

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