Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,149,790	C→G	7.3%	H333D (CAT→GAT) ‡	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,790	0	C	G	7.3%	197.4 / 4.1	65	H333D (CAT→GAT) ‡	PP_2765	sulfonate monooxygenase MsuD
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base G (2/3);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAAGTCGACCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTG  >  NC_002947/3149651‑3149938                                                                                                                                            |                                                                                                                                                      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                                                                                                                                      <  7:172057/91‑1 (MQ=255)                                                                   gCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTc                                                                                                                                      >  8:172057/1‑91 (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:62771/6‑109 (MQ=255)                                                                                                                                gacctcTTCCTGGAAATGCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                      >  6:37317/6‑109 (MQ=255)                                                                                                                                gacctcTTCCTGGAAATGCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                      <  5:62771/104‑1 (MQ=255)                                                                                                                                gacctcTTCCTGGAAATGCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAg                                                      <  5:37317/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                                                                                      >  2:136281/1‑67 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTggagca                                                                                      <  1:136281/67‑1 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGgc                                 >  3:216362/1‑120 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGgc                                 <  4:216362/120‑1 (MQ=255)                                                                                                                                          tGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGc    >  2:240385/1‑149 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCTTGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg  <  6:344013/149‑1 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg  <  6:349165/149‑1 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg  <  6:143507/149‑1 (MQ=255)                                                                                                                                            cATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTg  >  4:258373/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                      GGCCAAGTCGACCTTCGAAGACCTGGTGCAGTACAACGACGGCTTCAAGACCAACCTGATCGGTACGCCACGGCAGATCGCCGAGCGCATCGTCGCGCTCAAGGCGGTCGGCGTCGACCTGATCCTCTCCGGCTTCCTGCATTTCCAGGAAGAGGTCGAGTACTTC‑GGCAAGCATGTGTTGCCGCTGGTGCGCGAACTGGAGCAGGAGCAGCATGCCTCCCTCGCTGTGGCCTGAGGTGCGCGGCCATGCATGAGGCGCTCGCGAACAACCTGCCGGCCAGCTTGCTG  >  NC_002947/3149651‑3149938

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