Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	618,351	T→G	5.6%	Y73S (TAC→TCC)	PP_0532 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	618,351	0	T	G	5.6%	228.2 / 4.1	71	Y73S (TAC→TCC)	PP_0532	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base G (2/2);  total (33/38)
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.53e-01

GGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTC  >  NC_002947/618211‑618493                                                                                                                                              |                                                                                                                                               ggccggccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGAt                                                                                                                                        <  7:235217/149‑1 (MQ=255)  gccggccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATg                                                                                                                                       >  2:16741/1‑149 (MQ=255)     ggccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTg                                                                                                                                    <  1:281360/149‑1 (MQ=255)      gccCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTgg                                                                                                                                   >  7:253996/1‑149 (MQ=255)            tCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCCGGCGCTg                                                                                                                              <  4:78171/148‑1 (MQ=255)              gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGa                                                                                                                            <  8:12319/148‑1 (MQ=255)              gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGa                                                                                                                            <  1:16741/148‑1 (MQ=255)              gATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGAc                                                                                                                           >  4:237111/1‑149 (MQ=255)                      gTTCGGTCCGGTAACCCACCGCCCCCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGc                                                                                                                   <  5:163544/149‑1 (MQ=255)                       ttCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGa                                                                                                                         >  1:17797/1‑142 (MQ=255)                       ttCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGa                                                                                                                         <  2:17797/142‑1 (MQ=255)                                     ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgt                                                                                                                                   <  6:331356/118‑2 (MQ=255)                                     ccACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCATCGGGGTACCCGGTGCTgt                                                                                                                                   >  5:331356/1‑117 (MQ=255)                                     ccACCGCAACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAgg                                                                                                    >  4:161964/1‑149 (MQ=255)                                      cACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGc                                                                                                   >  1:109098/1‑149 (MQ=255)                                      cACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGAATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACAATTCGGCTTGACAACCCACAGGc                                                                                                   >  8:267417/1‑149 (MQ=255)                                       aCCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCg                                                                                                  >  6:192197/1‑149 (MQ=255)                                           ccaccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCAGGCTACCCCAGCGCAGCATTCCCATGCACCATCAACAGCACCGGGTACCCCGATGGTGGCGCGGGACGTTTCGGCTTGCCAACCCCCGGGCgg                                                                                                 >  1:250314/1‑147 (MQ=255)                                            caccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGc                                                                                                                   <  4:80619/127‑1 (MQ=255)                                            caccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGc                                                                                                                   >  3:80619/1‑127 (MQ=255)                                            caccaGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGt                                                                                             <  5:192197/149‑1 (MQ=255)                                               cagcagGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGAc                                                                                          <  3:195648/149‑1 (MQ=255)                                                              gCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCAGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGc                                                                             <  3:299331/147‑1 (MQ=255)                                                               cccACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGa                                                                          <  2:250314/149‑1 (MQ=255)                                                                 cACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGc                                                                                    >  3:20178/1‑137 (MQ=255)                                                                 cACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGc                                                                                    <  4:20178/137‑1 (MQ=255)                                                                      cGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTc                                                                       >  7:98904/1‑145 (MQ=255)                                                                       gAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGAAGCATTCACATCGAGCATCCACAGCACCGGCTACCACGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTc                                                                       >  5:212094/1‑144 (MQ=255)                                                                        aGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGAt                                                                  >  4:290011/1‑148 (MQ=255)                                                                         gCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATc                                                                 >  5:116195/1‑148 (MQ=255)                                                                          ccTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATc                                                                 >  7:265207/1‑147 (MQ=255)                                                                          ccTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCaa                                                               <  3:237111/149‑1 (MQ=255)                                                                           cTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATc                                                                 >  5:130505/1‑146 (MQ=255)                                                                           cTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCaa                                                               <  6:130505/148‑1 (MQ=255)                                                                             gTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCgg                                                                                                                    >  8:53729/1‑93 (MQ=255)                                                                             gTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCgg                                                                                                                    <  7:53729/93‑1 (MQ=255)                                                                             gTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCt                                                            >  7:137269/1‑149 (MQ=255)                                                                                 aGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGt                                                        >  8:48814/1‑149 (MQ=255)                                                                                       tgctgATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGCTCAAGCTGGGTGAAGt                                                   >  7:260652/1‑148 (MQ=255)                                                                                        gctgATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGa                                                 >  3:306034/1‑149 (MQ=255)                                                                                         ctgATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGAt                                                <  2:111516/149‑1 (MQ=255)                                                                                            aTCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAg                                              <  1:221843/148‑1 (MQ=255)                                                                                               cA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCa                                          <  4:306034/149‑1 (MQ=255)                                                                                                a‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGTGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCTGGCCATCGTCAGAGTCCAGATCAAGCTGTGTGAAGTGATAGGCCAg                                         <  6:212094/149‑1 (MQ=255)                                                                                                a‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGcc                                                                            <  8:71182/114‑1 (MQ=255)                                                                                                a‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGcc                                                                            >  7:71182/1‑114 (MQ=255)                                                                                                     tCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTcc                                     >  8:51649/1‑149 (MQ=255)                                                                                                        ccAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGcc                                                                            >  6:122794/1‑107 (MQ=255)                                                                                                        ccAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTgg                                  <  7:48814/149‑1 (MQ=255)                                                                                                        ccAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTgg                                  <  8:265207/149‑1 (MQ=255)                                                                                                            cgcAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCt                              <  3:161964/149‑1 (MQ=255)                                                                                                              cagcaTTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGt                            <  8:253996/149‑1 (MQ=255)                                                                                                                gcaTTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCg                                                                                                  >  6:144481/1‑77 (MQ=255)                                                                                                                gcaTTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCg                                                                                                  <  5:144481/77‑1 (MQ=255)                                                                                                                   ttCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAacac                       <  6:116195/149‑1 (MQ=255)                                                                                                                      ccATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTcc                                     >  3:7270/1‑132 (MQ=255)                                                                                                                      ccATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTcc                                     <  4:7270/132‑1 (MQ=255)                                                                                                                      ccATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGa                    >  7:23433/1‑149 (MQ=255)                                                                                                                      ccATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGa                    <  7:51649/149‑1 (MQ=255)                                                                                                                         tCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTcc                                     <  3:2349/129‑1 (MQ=255)                                                                                                                         tCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTcc                                     >  4:2349/1‑129 (MQ=255)                                                                                                                         tCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACgt                 <  8:137269/149‑1 (MQ=255)                                                                                                                               aTCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAt           <  3:223191/149‑1 (MQ=255)                                                                                                                                   ccaGCATCGGGGTACCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTg                             <  2:4173/126‑1 (MQ=255)                                                                                                                                   ccaGCATCGGGGTACCCGGTGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTg                             >  1:4173/2‑127 (MQ=255)                                                                                                                                   acaGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGCCGCAGGCCATCGGCCGACTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCAtcttc       >  1:79306/1‑149 (MQ=255)                                                                                                                                     aGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTg     <  3:290011/149‑1 (MQ=255)                                                                                                                                      gCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTgg    <  8:23433/149‑1 (MQ=255)                                                                                                                                        aCCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTc  <  2:109098/149‑1 (MQ=255)                                                                                                                                         ccGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                    >  1:172849/1‑98 (MQ=255)                                                                                                                                         ccGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAg                                                    <  2:172849/98‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               GGCCGGCCCTGTCGATGCGCAGTTCGGTCCGGTAACCCACCGCCACCAGCAGGGGCGCCTGGCCCACGGCGAGCCTGTCCAGTTGCTGCTGATCCA‑GCGTCCCCAGCGCAGCATTCCCATCGAGCATCCACAGCACCGGGTACCCCGATGCTGGCGCTGGACGATTCGGCTTGCCAACCCACAGGCGGTAGTGACGCAGGCCATCGGCCGAGTCCAGATCAAGCTGGGTGAAGTGATAGGCCAGGTCCTGGCGCTGTAACACCGACGTGTCCATCTTCTGGTC  >  NC_002947/618211‑618493

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