Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,149,795	C→G	6.5%	F334L (TTC→TTG)	PP_2765 →	sulfonate monooxygenase MsuD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,149,795	0	C	G	6.5%	218.5 / 2.1	69	F334L (TTC→TTG)	PP_2765	sulfonate monooxygenase MsuD
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (2/3);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAAGTCGACCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTGGCCCG  >  NC_002947/3149651‑3149943                                                                                                                                                 |                                                                                                                                                      ggcttaGTCGACCTTCGTAGAGCTTGTGCAGTACATCGACGGCTTCTAGACCTACCTGATCGGTACTCCACGGCAGATCGCCGAGCGCATCGTCGGTCTCAAGGCTGTCGGCGTCGACCTGATCCTCTGCGGCTTCCTGCATTTCCAgt                                                                                                                                                   <  4:47020/144‑2 (MQ=255) ggCCAAGTCGACCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAgg                                                                                                                                                   <  7:111996/149‑1 (MQ=255)     aaGTCGACCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGaa                                                                                                                                                 >  8:203670/1‑147 (MQ=255)          gACCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTc                                                                                                                                           <  8:139480/148‑1 (MQ=255)           aCCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGa                                                                                                                                         <  2:214602/149‑1 (MQ=255)             cTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAg                                                                                                                                        >  3:324859/1‑148 (MQ=255)              ttCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGACGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTa                                                                                                                                      >  3:323209/1‑149 (MQ=255)                cGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCAACGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACt                                                                                                                                    >  8:247075/1‑149 (MQ=255)                     aCCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGc                                                                                                                              <  7:203670/149‑1 (MQ=255)                             cAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCAtgt                                                                                                                      <  3:178577/149‑1 (MQ=255)                                tACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGtt                                                                                                                   >  5:105237/1‑149 (MQ=255)                                           cTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTgc                                                                                                        >  5:131831/1‑149 (MQ=255)                                               aaGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAaga                                                                                                                                               <  8:83692/107‑1 (MQ=255)                                               aaGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAaga                                                                                                                                               >  7:83692/1‑107 (MQ=255)                                               aaGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACtt                                                                                                                                   <  5:269149/119‑1 (MQ=255)                                               aaGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGcgcg                                                                                                     >  8:1164/1‑148 (MQ=255)                                               aaGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGAGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACtt                                                                                                                                   >  6:269149/1‑119 (MQ=255)                                                     aaCCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTg                                                                                                >  5:62993/1‑147 (MQ=255)                                                       ccTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagc                                                                                            <  6:131831/149‑1 (MQ=255)                                                         tGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagcag                                                                                          >  2:288694/1‑149 (MQ=255)                                                           aTCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagcagga                                                                                        <  1:288694/149‑1 (MQ=255)                                                           aTCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagcagg                                                                                         <  2:95231/148‑1 (MQ=255)                                                           aTCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagcag                                                                                          >  1:179089/1‑147 (MQ=255)                                                             cGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagcaggag                                                                                       >  5:359162/1‑148 (MQ=255)                                                                   gCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTc                                                                                                                                           >  8:172057/1‑91 (MQ=255)                                                                   gCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTc                                                                                                                                           <  7:172057/91‑1 (MQ=255)                                                                      aCGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGcct                                                                             <  3:31469/149‑1 (MQ=255)                                                                           aGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagcaggagc                                                                                      <  7:292088/135‑1 (MQ=255)                                                                           aGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagcaggagc                                                                                      >  8:292088/1‑135 (MQ=255)                                                                           aGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGcctccctc                                                                        >  6:914/1‑149 (MQ=255)                                                                                  cGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTgg                                                                 <  6:62993/149‑1 (MQ=255)                                                                                     gcgcATCGGCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCGTCCTGCATTTCCAGGAAGAGGTCGTGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCt                                                              <  5:80729/149‑1 (MQ=255)                                                                                        cATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTAGTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGTAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                           <  4:324859/149‑1 (MQ=255)                                                                                        cATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                           >  7:13120/1‑149 (MQ=255)                                                                                         aTCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAgg                                                          >  8:142634/1‑149 (MQ=255)                                                                                             tcgCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGACCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTgcg                                                      >  4:104480/1‑149 (MQ=255)                                                                                                cgcTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCgg                                                   <  8:13120/149‑1 (MQ=255)                                                                                                       ggCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCcatgca                                            >  7:28654/1‑149 (MQ=255)                                                                                                        gCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCcatgcat                                           <  6:178658/149‑1 (MQ=255)                                                                                                               gCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGgcgc                                    >  4:209198/1‑149 (MQ=255)                                                                                                                 gTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTCGGGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCt                                   >  1:94288/1‑149 (MQ=255)                                                                                                                  tCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGcc                                                 >  7:307733/1‑133 (MQ=255)                                                                                                                  tCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGcc                                                 <  8:307733/133‑1 (MQ=255)                                                                                                                     aCCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGa                              <  2:54475/149‑1 (MQ=255)                                                                                                                      ccTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGaa                             >  2:281072/1‑149 (MQ=255)                                                                                                                         gATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGaacaa                          >  5:143507/1‑149 (MQ=255)                                                                                                                          aTCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGaacaac                         <  3:209198/149‑1 (MQ=255)                                                                                                                           tCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAAcc                        <  7:247075/149‑1 (MQ=255)                                                                                                                             ctctCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTg                      >  5:47298/1‑149 (MQ=255)                                                                                                                             ctctCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTg                      >  5:235037/1‑149 (MQ=255)                                                                                                                                tcCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGc                     <  3:104480/147‑1 (MQ=255)                                                                                                                                gacctcTTCCTGGAAATGCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                           >  6:37317/6‑109 (MQ=255)                                                                                                                                gacctcTTCCTGGAAATGCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                           >  6:62771/6‑109 (MQ=255)                                                                                                                                gacctcTTCCTGGAAATGCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                           <  5:37317/104‑1 (MQ=255)                                                                                                                                gacctcTTCCTGGAAATGCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                           <  5:62771/104‑1 (MQ=255)                                                                                                                                  cGGCTTCCTGGATTTGCAGGAAGAGGTCGTGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGGGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGGTTGAGGCGCTCGCGAACAACCTGCCg                   <  4:323209/147‑1 (MQ=255)                                                                                                                                  cGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCgg                  <  8:28654/148‑1 (MQ=255)                                                                                                                                  cGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGc                 >  7:131721/1‑149 (MQ=255)                                                                                                                                        ccTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCcat                                               >  2:320088/1‑113 (MQ=255)                                                                                                                                        ccTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCcat                                               <  1:320088/113‑1 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagca                                                                                           <  1:136281/67‑1 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagca                                                                                           >  2:136281/1‑67 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGgc                                      <  4:216362/120‑1 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGgc                                      >  3:216362/1‑120 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGc         >  2:240385/1‑149 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCTTGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg       <  6:344013/149‑1 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg       <  6:143507/149‑1 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg       <  6:349165/149‑1 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg       >  4:258373/1‑149 (MQ=255)                                                                                                                                                 ccAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTGGCCCg  >  1:157134/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GGCCAAGTCGACCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTGGCCCG  >  NC_002947/3149651‑3149943

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