Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,611,237	C→G	5.7%	R77R (CGC→CGG)	PP_4931 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,611,237	0	C	G	5.7%	226.5 / 3.6	70	R77R (CGC→CGG)	PP_4931	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base G (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GGCTACTCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGAC  >  NC_002947/5611091‑5611384                                                                                                                                                   |                                                                                                                                                    ggCTACTCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTg                                                                                                                                                   >  7:267825/1‑149 (MQ=255)   cTACTCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGgc                                                                                                                                                 <  2:38392/149‑1 (MQ=255)      cTCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGgcgc                                                                                                                                               >  6:160561/1‑148 (MQ=255)      cTCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGgcgc                                                                                                                                               >  7:38565/1‑148 (MQ=255)       tCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGTCGCGCCAGCGGcgcgg                                                                                                                                             >  2:100253/1‑148 (MQ=255)       tCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGTCGCGCCAGCGGcgcgg                                                                                                                                             >  2:95392/1‑148 (MQ=255)       tCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGgcgc                                                                                                                                               >  5:285998/1‑147 (MQ=255)         gACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAg                                                                                                                                            <  5:164926/148‑1 (MQ=255)           cAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGcc                                                                                                                                         <  8:38565/149‑1 (MQ=255)            aGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGTCGCGCCAGCGGcgcggcatgt                                                                                                                                        <  1:100253/149‑7 (MQ=255)            aGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGTCGCGCCAGCGGcgcggcatgt                                                                                                                                        <  1:95392/149‑7 (MQ=255)              tGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTg                                                                                                                                      <  8:93428/149‑1 (MQ=255)              tGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTg                                                                                                                                      <  2:192675/149‑1 (MQ=255)                 tggtggAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGcc                                                                                                                                   >  3:366301/1‑149 (MQ=255)                   gtggAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGa                                                                                                                                             <  5:73165/137‑1 (MQ=255)                   gtggAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGa                                                                                                                                             >  6:73165/1‑137 (MQ=255)                       aGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAAc                                                                                                                             >  8:35793/1‑149 (MQ=255)                        gcgcgcCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACg                                                                                                                            >  8:80038/1‑149 (MQ=255)                        gcgcgcCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACg                                                                                                                            >  1:304895/1‑149 (MQ=255)                          gcgcCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGc                                                                                                                          <  3:320160/149‑1 (MQ=255)                           cgcCAGCCGTGGCGGCGAGCAGCCCCTCACGGCGCAGGGCCTCAGCCACGGCGGCACCAACTCCGCCCTGCACGGCGCCCTGCCCGGCGCCTCGGACGCCCTCGCCGCCATGCCGCGCCGCTGGCGCGAGGCCGTGGCCGGCAACGGCg                                                                                                                         <  3:293278/149‑1 (MQ=255)                            gcCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGCGGCCTTGGCCGGCAACGGCg                                                                                                                         >  4:191969/1‑148 (MQ=255)                             ccAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGaa                                                                                                                       <  4:19881/149‑1 (MQ=255)                                gCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCCCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTc                                                                                                                    <  5:298074/149‑1 (MQ=255)                                 ccTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTc                                                                                                                    >  3:84116/1‑148 (MQ=255)                                       cggcgAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCa                                                                                                             <  8:267825/149‑1 (MQ=255)                                          cgAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTc                                                                                                          <  3:126309/149‑1 (MQ=255)                                           gAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGACCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCAACGCCTACATGCCGCGCCGCTGGAGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCg                                                                                                         >  7:170452/1‑149 (MQ=255)                                                 aCCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCg                                                                                                   >  6:199278/1‑149 (MQ=255)                                                  ccATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGc                                                                                                  <  2:94802/149‑1 (MQ=255)                                                    aTCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCt                                                                                                <  2:304317/149‑1 (MQ=255)                                                      cAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCt                                                                                                <  1:194616/147‑1 (MQ=255)                                                      cAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCt                                                                                                >  2:194616/1‑147 (MQ=255)                                                          tCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCtctc                                                                                          >  7:262240/1‑149 (MQ=255)                                                            gCAGGGCATCATCCCCGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTcc                                                                                         >  8:265924/1‑148 (MQ=255)                                                              aGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGAt                                                                                      <  2:243146/149‑1 (MQ=255)                                                               gggCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATg                                                                                     >  5:337122/1‑149 (MQ=255)                                                                ggCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGc                                                                                    <  7:35793/149‑1 (MQ=255)                                                                 gCATCATCCACGGCGGCACCAAATACGCCCTGCACGGAGCAATGACCCGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGCCGAACTCGACATCAGTCGTACCCGCCTCCTCTCAGATGcc                                                                                   >  8:188082/1‑149 (MQ=255)                                                                  catcatCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGcc                                                                                   >  8:54201/1‑148 (MQ=255)                                                                      atcCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCAGATGCCCACTa                                                                              <  7:188082/149‑1 (MQ=255)                                                                          aCGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCtg                                                                          >  8:107179/1‑149 (MQ=255)                                                                          aCGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCtg                                                                          <  4:366301/149‑1 (MQ=255)                                                                                 cACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGc                                                                   <  7:54201/149‑1 (MQ=255)                                                                                  aCCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGcc                                                                  >  7:30153/1‑149 (MQ=255)                                                                                    cataTACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTGAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCgg                                                                <  8:170452/146‑1 (MQ=255)                                                                                    cAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGcc                                                                                                 >  5:71987/1‑116 (MQ=255)                                                                                    cAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGcc                                                                                                 <  6:71987/116‑1 (MQ=255)                                                                                       aTACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCa                                                             >  8:11597/1‑149 (MQ=255)                                                                                         aCGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCAcc                                                           >  7:177165/1‑149 (MQ=255)                                                                                          cGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCt                                                          >  1:284026/1‑149 (MQ=255)                                                                                             ccTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTgg                                                       >  5:276970/1‑149 (MQ=255)                                                                                                 cACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTg                                                        >  5:249588/1‑144 (MQ=255)                                                                                                 cACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTg                                                        <  6:249588/144‑1 (MQ=255)                                                                                                       gcCCTGACCGGCGCCGCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTTGCCGGCCAAGGTGCCCGGCGAcca                                             <  7:80038/149‑2 (MQ=255)                                                                                                         ccTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTca                                           >  6:76410/1‑149 (MQ=255)                                                                                                         ccTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTca                                           >  5:156347/1‑149 (MQ=255)                                                                                                          cTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTcac                                          <  8:262240/149‑1 (MQ=255)                                                                                                           tGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTcac                                          >  4:237162/1‑148 (MQ=255)                                                                                                             aCCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTcaca                                         <  6:337122/147‑1 (MQ=255)                                                                                                             aCCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAg                                       <  2:104741/149‑1 (MQ=255)                                                                                                              ccGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGc                                      >  8:38577/1‑149 (MQ=255)                                                                                                                 gcgcCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGcttc                                   >  6:80758/1‑149 (MQ=255)                                                                                                                   gcCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGcttctt                                 <  8:177165/149‑1 (MQ=255)                                                                                                                    ccTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCttcttc                                <  5:160561/149‑1 (MQ=255)                                                                                                                         gAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTAc                                                                             <  6:191444/99‑1 (MQ=255)                                                                                                                         gAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTAc                                                                             >  5:191444/1‑99 (MQ=255)                                                                                                                               aTCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGc                     >  5:104477/1‑149 (MQ=255)                                                                                                                               aTCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGc                     >  2:136679/1‑149 (MQ=255)                                                                                                                                       cATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGcgcg             <  6:156347/149‑1 (MQ=255)                                                                                                                                           ccgcgccgcTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCg         >  8:86602/1‑149 (MQ=255)                                                                                                                                             gcgccgcTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCggg       >  8:260918/1‑149 (MQ=255)                                                                                                                                              cgccgcTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGGCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGt      <  3:191969/149‑1 (MQ=255)                                                                                                                                                 cgcTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGa   >  3:122805/1‑149 (MQ=255)                                                                                                                                                  gcTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGAc  <  4:122805/149‑1 (MQ=255)                                                                                                                                                   cTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTcac                                          >  5:36683/1‑108 (MQ=255)                                                                                                                                                   cTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTcac                                          <  6:36683/108‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GGCTACTCGACAGTGCTGGTGGAGCGCGCCAGCCTTGGCGGCGAGCAGACCATCAAGTCGCAGGGCATCATCCACGGCGGCACCAAATACGCCCTGCACGGCGCCCTGACCGGCGCCTCGGAAGCCATCGCCGACATGCCGCGCCGCTGGCGCGAGGCCTTGGCCGGCAACGGCGAACTCGACCTCAGTCGTACCCGCCTGCTCTCCGATGCCCACTACCTGTGGTCGCCGGGCACCTTGGCCGGCAACCTCACAAGCTTCTTCGCCAGCAAGGCCGTGCGCGGCCGGGTCGAC  >  NC_002947/5611091‑5611384

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