Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,682,935	A→G	8.0%	Q485R (CAG→CGG)	recJ →	ssDNA‑specific exonuclease RecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,682,935	0	A	G	8.0%	174.4 / 4.0	62	Q485R (CAG→CGG)	recJ	ssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (2/3);  total (29/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 = 1.00e+00

GCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGT  >  NC_002947/1682787‑1683079                                                                                                                                                      |                                                                                                                                                  gCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCa                                                                                                                                                   <  5:295825/149‑1 (MQ=255) gCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCa                                                                                                                                                   <  5:247290/149‑1 (MQ=255)  cGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAg                                                                                                                                                  <  5:296246/149‑1 (MQ=255)        ccGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCTGGCCCCATGGGCCGGCatgg                                                                                                                                             <  5:363916/148‑4 (MQ=255)        ccGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTc                                                                                                                                            >  5:91465/1‑149 (MQ=255)               gTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACTGGGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAg                                                                                                                                      >  1:86584/1‑149 (MQ=255)               gTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCCATGGGCCGGCat                                                                                                                                               >  6:106323/1‑138 (MQ=255)               gTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCCATGGGCCGGCat                                                                                                                                               <  5:106323/139‑2 (MQ=255)               gTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGc                                                                                                                                     <  5:228313/149‑1 (MQ=255)                    acgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTg                                                                                                                                >  1:230435/1‑149 (MQ=255)                    acgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTg                                                                                                                                <  8:237121/149‑1 (MQ=255)                    acgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTg                                                                                                                                >  6:78617/1‑149 (MQ=255)                    acgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTg                                                                                                                                >  6:73347/1‑149 (MQ=255)                     cgaAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGt                                                                                                                               <  4:7003/149‑1 (MQ=255)                           ggTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATg                                                                                                                        <  4:305895/149‑1 (MQ=255)                              gcggcgTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCg                                                                                                                     >  5:279728/1‑149 (MQ=255)                                 gcgTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCgtgt                                                                                                                  >  4:132550/1‑149 (MQ=255)                                       gCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGt                                                                                                            >  6:265224/1‑149 (MQ=255)                                        cggcgcGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGtt                                                                                                           <  8:262626/147‑1 (MQ=255)                                        cTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGtt                                                                                                           >  7:218220/1‑149 (MQ=255)                                          gcgcGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTtgg                                                                                                         <  8:303704/149‑1 (MQ=255)                                              gaagaaGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGa                                                                                                     >  4:195212/1‑149 (MQ=255)                                                  aagaCCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAg                                                                                                                                                  <  1:85336/101‑1 (MQ=255)                                                  aagaCCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAg                                                                                                                                                  >  2:85336/1‑101 (MQ=255)                                                     aCCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCgtgt                                                                                                                  <  8:61576/129‑1 (MQ=255)                                                     aCCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCgtgt                                                                                                                  >  7:61576/1‑129 (MQ=255)                                                     aCCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCgt                                                                                              <  2:230435/149‑1 (MQ=255)                                                      ccTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTCGGTGGAGCAGCgtg                                                                                             <  6:279728/149‑1 (MQ=255)                                                         gACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTtggtgg                                                                                                      <  4:326397/137‑1 (MQ=255)                                                         gACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTtggtgg                                                                                                      >  3:326397/1‑137 (MQ=255)                                                                  gccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGaa                                                                                  >  5:35178/1‑149 (MQ=255)                                                                   ccTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAAc                                                                                 <  6:91465/149‑1 (MQ=255)                                                                          tCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAg                                                                                                                                                  <  1:150450/78‑1 (MQ=255)                                                                          tCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGACATGCCCGCCCATGGGGCCAg                                                                                                                                                  >  2:150450/1‑78 (MQ=255)                                                                            cGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTg                                                                         <  4:170761/148‑1 (MQ=255)                                                                            cGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTA‑CATGGAGTGTTCCAGGAGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTg                                                                         >  6:25488/1‑148 (MQ=255)                                                                                 gCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAgtg                                                                   <  3:195212/149‑1 (MQ=255)                                                                                   aGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAgtggt                                                                 <  5:82437/149‑1 (MQ=255)                                                                                      cTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCt                                                              <  8:218220/149‑1 (MQ=255)                                                                                         gCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTgg                                                                                          <  7:336970/118‑1 (MQ=255)                                                                                         gCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTgg                                                                                          >  8:336970/1‑118 (MQ=255)                                                                                         gCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCaa                                                           <  5:78617/149‑1 (MQ=255)                                                                                         gCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCaa                                                           <  5:73347/149‑1 (MQ=255)                                                                                           ggTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAaga                                                         <  3:132550/149‑1 (MQ=255)                                                                                                  gagTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATgg                                                                                                                       >  3:1761/1‑80 (MQ=255)                                                                                                  gagTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATgg                                                                                                                       <  4:1761/80‑1 (MQ=255)                                                                                                                   ggCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTg                                  >  4:132819/1‑148 (MQ=255)                                                                                                                     ccAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGAt                               <  6:35178/149‑1 (MQ=255)                                                                                                                      cAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAaga                                                         <  8:79049/122‑1 (MQ=255)                                                                                                                      cAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAaga                                                         >  7:79049/1‑122 (MQ=255)                                                                                                                      cAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATg                              <  2:86584/149‑1 (MQ=255)                                                                                                                       aaGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGaa                                                    <  4:12467/126‑1 (MQ=255)                                                                                                                       aaGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGaa                                                    >  3:12467/1‑126 (MQ=255)                                                                                                                         ggCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTG‑CATGTCGTGTTCCAGTTGGTGGCGCAGCGTGTGGTTGGCGACCGGCACCTGAAAGTGGTGCTGAAGCGCGAATGTGGCTCGGTGCGGCTGGATGGCa                           >  5:258043/1‑149 (MQ=255)                                                                                                                            ccTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGGGAATGTGGCTCGGTGCGGCTGGATGGCATTg                        >  8:99067/1‑149 (MQ=255)                                                                                                                            ccTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTg                        <  5:265224/149‑1 (MQ=255)                                                                                                                                  ccATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTCGCATGGCGTGTGCCAGTTGGTGGGGCAGCGTGTGGTTGGCGACCGGCCCCTGAAAGTGTTGCTCAAGAGCGAATGTGGCTGGGTGCGGCTGGATGGCATTGCCTTc                   >  1:256632/1‑149 (MQ=255)                                                                                                                                      ctggccccAT‑GGGCCGGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAagag                                                        >  5:87323/6‑106 (MQ=255)                                                                                                                                      ctggccccAT‑GGGCCGGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAagag                                                        <  6:87323/101‑1 (MQ=255)                                                                                                                                        cGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCAtt             >  1:364467/1‑148 (MQ=255)                                                                                                                                        cGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCAtt             >  6:334457/1‑148 (MQ=255)                                                                                                                                          gCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGAc          >  7:225544/1‑149 (MQ=255)                                                                                                                                                 gggCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGaa    <  1:319792/148‑1 (MQ=255)                                                                                                                                                  ggCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAgt  <  2:256632/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GCGTTTGCCGAGGCGTTCGACGAAGAGGTGCGGCGTCAGCTGCGCGAAGAAGACCTGACT‑GGCCGCCTGTTGTCCGACGGCAGCCTGGCGGTAGAGGAGTTCCACCTCGACCTGGCCAAGGCCCTGCGCCATGCCGGCCCATGGGGCCAGCACTTCCCTGAGCCGCTGTTC‑CATGGCGTGTTCCAGTTGGTGGAGCAGCGTGTGGTTGGCGAACGGCACCTGAAAGTGGTGCTCAAGAGCGAATGTGGCTCGGTGCGGCTGGATGGCATTGCCTTCGGCATTGACCGCGAAGT  >  NC_002947/1682787‑1683079

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