Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,042,886	A→C	6.6%	M275L (ATG→CTG)	PP_5295 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,042,886	0	A	C	6.6%	219.3 / 4.6	74	M275L (ATG→CTG)	PP_5295	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base C (2/3);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCGCCTGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAGGCCACCCCGACCATGAA  >  NC_002947/6042744‑6043032                                                                                                                                               |                                                                                                                                                      acaGCGCCTGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCTTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAgg                                                                                                                                                 <  1:269305/149‑1 (MQ=255)       ccTGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGCATCAGGAAGTc                                                                                                                                            >  5:232900/1‑148 (MQ=255)        cTGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGCATCAGGAAGTcca                                                                                                                                          <  8:332396/149‑3 (MQ=255)         tGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGa                                                                                                                                                       <  4:5579/135‑1 (MQ=255)         tGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGa                                                                                                                                                       >  3:5579/1‑135 (MQ=255)         tGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAAc                                                                                                                                         >  2:241251/1‑149 (MQ=255)           gAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGCATCAGGAAGTccag                                                                                                                                         <  6:232900/147‑4 (MQ=255)           gAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGCATCAGGAAGTccag                                                                                                                                         <  7:250160/147‑4 (MQ=255)           gAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACcg                                                                                                                                       >  6:207582/1‑149 (MQ=255)            aGCCGGAGATGGGACTGCTGGCCCAGAAAAGCGACGGGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACcgc                                                                                                                                      <  4:304989/149‑1 (MQ=255)               aGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAg                                                                                                                                   <  4:170825/149‑1 (MQ=255)                ggAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAgg                                                                                                                                  <  3:305928/149‑1 (MQ=255)                 gagaTGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGc                                                                                                                                 <  4:86991/149‑1 (MQ=255)                 gagaTGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGc                                                                                                                                 >  4:50674/1‑149 (MQ=255)                 gagaTGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGc                                                                                                                                 <  8:106146/149‑1 (MQ=255)                             ctgGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTc                                                                                                                       <  4:229214/147‑1 (MQ=255)                             ctgGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCgg                                                                                                                     <  6:47221/149‑1 (MQ=255)                                  ccAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAAcc                                                                                                                                        <  4:311966/125‑1 (MQ=255)                                  ccAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAAcc                                                                                                                                        >  3:311966/1‑125 (MQ=255)                                  ccAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAacac                                                                                                                           <  7:65944/138‑1 (MQ=255)                                  ccAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAacac                                                                                                                           >  8:65944/1‑138 (MQ=255)                                  ccAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACATGCGGGTGCAGGAACTCAACCGCGAGGGCCACACCCTCGGCGcc                                                                                                                 >  2:330610/1‑148 (MQ=255)                                   cAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTcc                                                                                                                 >  3:227512/1‑147 (MQ=255)                                    aGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCCCGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAg                                                                                                              <  1:330610/149‑1 (MQ=255)                                    aGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAg                                                                                                              <  4:227512/149‑1 (MQ=255)                                    aGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGATGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAg                                                                                                              >  2:320760/1‑149 (MQ=255)                                            gACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGa                                                                                                      <  1:56014/149‑1 (MQ=255)                                               gTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGAccc                                                                                                   >  2:238967/1‑149 (MQ=255)                                                  gCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGc                                                                                                                                 <  3:1840/116‑1 (MQ=255)                                                  gCAGAGGAACTCGACCGCCTCACCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGc                                                                                                                                 >  4:1840/1‑116 (MQ=255)                                                     gagGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAAc                                                                                                                                              >  7:212239/1‑100 (MQ=255)                                                     gagGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAAc                                                                                                                                              <  8:212239/100‑1 (MQ=255)                                                         aaCTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACg                                                                                         >  8:236591/1‑149 (MQ=255)                                                           cTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCt                                                                                                          <  3:308261/130‑1 (MQ=255)                                                           cTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCt                                                                                                          >  4:308261/1‑130 (MQ=255)                                                            tCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGc                                                                                        >  1:205353/1‑147 (MQ=255)                                                             cGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACg                                                                                                >  1:306711/1‑138 (MQ=255)                                                             cGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACg                                                                                                <  2:306711/138‑1 (MQ=255)                                                                      cAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGTGGCG‑GT‑‑c                                                                          >  6:170971/1‑148 (MQ=255)                                                                            ccACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTg                                                                    >  4:131709/1‑148 (MQ=255)                                                                              aCGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAg                                                                 <  1:238967/149‑1 (MQ=255)                                                                              aCGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAg                                                                 <  3:50674/149‑1 (MQ=255)                                                                              aCGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAg                                                                 <  1:241251/149‑1 (MQ=255)                                                                                 tCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTa                                                              <  4:123137/149‑1 (MQ=255)                                                                                  cACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTAt                                                             <  5:207582/149‑1 (MQ=255)                                                                                   aCCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTAtt                                                            >  3:342649/1‑149 (MQ=255)                                                                                   aCCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAACACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTAtt                                                            >  7:164536/1‑149 (MQ=255)                                                                                     cGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCCCGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGGATTGa                                                          >  5:200001/1‑149 (MQ=255)                                                                                     cGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCGAGTTGCAACCTGAAGGTAt                                                             >  1:303806/1‑149 (MQ=255)                                                                                     cGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGa                                                          <  1:357877/149‑1 (MQ=255)                                                                                     cGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACAACCTCGGCTCCAAGGCCTTCGACCCCCGCAGCACGCAAGCGGCG‑GT‑‑CAACAAGAAGGTATTg                                                           >  4:290343/1‑148 (MQ=255)                                                                                      gAAGTGCGCCGCGTGCTCAAGGATGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGAt                                                         <  1:320760/149‑1 (MQ=255)                                                                                       aaGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATc                                                        <  3:290343/149‑1 (MQ=255)                                                                                         gTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGa                                                      <  7:236591/149‑1 (MQ=255)                                                                                          tGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGaa                                                     <  3:131709/149‑1 (MQ=255)                                                                                            cgccgcGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAAc                                                    <  1:285721/148‑1 (MQ=255)                                                                                                 cgTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGTGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGc                                              <  5:170971/149‑1 (MQ=255)                                                                                                     cTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTg                                           >  2:267426/1‑148 (MQ=255)                                                                                                       cAAGGGTGGCGGCCCCGCTGGCCGGCGCCTGGACATCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCATCGACCCACCCACCCCGCCAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGaa                                         >  4:304705/1‑148 (MQ=255)                                                                                                        aaGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCCTGaa                                         <  4:342649/147‑1 (MQ=255)                                                                                                          gggTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAg                                     <  3:304705/149‑1 (MQ=255)                                                                                                                 ggcgCCGCTGGCCGGCGCCTGGACTTCCTGCTGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGCCCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGa                                                 >  6:85267/1‑130 (MQ=255)                                                                                                                 ggcgCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGa                                                 <  5:85267/130‑1 (MQ=255)                                                                                                                  gcgcCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGCCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGaa                              >  7:142983/1‑148 (MQ=255)                                                                                                                     ccgcTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAACACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATc                          <  8:164536/149‑1 (MQ=255)                                                                                                                          ggCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGAt                                                <  1:148432/122‑1 (MQ=255)                                                                                                                             cGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCATCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAgg                  >  6:159889/1‑149 (MQ=255)                                                                                                                             cGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAgg                  >  4:32642/1‑149 (MQ=255)                                                                                                                                  ccTGGACTTCCTGATGCAGGAACTCAACCGCGAGCCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAGGCCAcc             <  8:142983/149‑1 (MQ=255)                                                                                                                                     ggACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTa                                   <  6:270780/124‑1 (MQ=255)                                                                                                                                     ggACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTa                                   >  5:270780/1‑124 (MQ=255)                                                                                                                                     ggACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAGGCCAcccc           >  2:361033/1‑148 (MQ=255)                                                                                                                                       aCTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAg                                                  <  6:358429/107‑1 (MQ=255)                                                                                                                                       aCTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAg                                                  >  5:358429/1‑107 (MQ=255)                                                                                                                                         ttCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAGGCCACCCCGACCa      >  2:217136/1‑149 (MQ=255)                                                                                                                                            cTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAGGCCACCCCGACCATg    >  6:295223/1‑148 (MQ=255)                                                                                                                                               aTGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAGGCCACCCCGACCATGaa  <  5:307015/147‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                      ACAGCGCCTGGAGCAGGAGATGGTACTGCTGGCCCAGAAAAGCGACGTGGCAGAGGAACTCGACCGCCTCAGCACCCACGTCACCGAAGTGCGCCGCGTGCTCAAGGGTGGCGGCGCCGCTGGCCGGCGCCTGGACTTCCTGATGCAGGAACTCAACCGCGAGGCCAACACCCTCGGCTCCAAGGCCTTCGACCCACGCAGCACGCAAGCGGCG‑GT‑‑CAACCTGAAGGTATTGATCGAACAGATGCGTGAACAAGTACAGAACATCGAGTAAGGCCACCCCGACCATGAA  >  NC_002947/6042744‑6043032

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