Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,209,547	C→G	6.7%	G52G (GGC→GGG)	ytnA →	amino acid permease YtnA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,209,547	0	C	G	6.7%	193.4 / 4.0	60	G52G (GGC→GGG)	ytnA	amino acid permease YtnA
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TCTGCCCACGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGCGGGCTTC  >  NC_002947/1209410‑1209691                                                                                                                                          |                                                                                                                                                 tctGCCCACGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCcatcat                                                                                                                                       >  7:36509/1‑149 (MQ=255)   tGCCCACGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGGCCGGCcatttc                                                                                                                                       <  6:186807/147‑4 (MQ=255)     cccACGGCCAGGCCACTGACGGCGGCCCGCGCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCAtgttg                                                                                                                                   <  8:138534/149‑1 (MQ=255)     cccACGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCAtgtt                                                                                                                                    <  4:150424/148‑1 (MQ=255)      ccACGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCAtgttgt                                                                                                                                  >  1:150934/1‑149 (MQ=255)      ccACGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCAtgttgt                                                                                                                                  >  8:15797/1‑149 (MQ=255)        aCGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTcc                                                                                                                                >  2:262937/1‑149 (MQ=255)          ggccaggccaCTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTa                                                                                                                              >  6:293609/1‑149 (MQ=255)                  caCTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATcg                                                                                                                      >  7:238634/1‑149 (MQ=255)                     tGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTACCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGCTCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATcggcg                                                                                                                   >  7:101871/1‑149 (MQ=255)                          gcggcCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTg                                                                                                              <  2:44529/149‑1 (MQ=255)                               ccGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGACCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCAt                                                                                                         >  8:183712/1‑149 (MQ=255)                                cGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCAtt                                                                                                        <  5:129503/149‑1 (MQ=255)                                 gCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGGTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTc                                                                                                       <  6:25582/149‑1 (MQ=255)                                     aaGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGt                                                                                                   <  4:128876/149‑1 (MQ=255)                                     aaGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTAGCGGTGTCGGTCTTTTCCTCGGTTCGGCCATGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCGGACATCATCGGCGGCCTGGCCATGCTGGt                                                                                                   >  7:207473/1‑149 (MQ=255)                                           gAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTgg                                                                                                    <  1:308178/142‑1 (MQ=255)                                           gAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTgg                                                                                                    >  2:308178/1‑142 (MQ=255)                                                 ggCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATgcgcgc                                                                                       <  4:4293/149‑1 (MQ=255)                                                     aaCGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCGTCATGTTGTCCTACATCATCGGCGGCCTGGc                                                                                                            >  8:295648/1‑124 (MQ=255)                                                      aCGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTg                                                                                  <  7:15797/149‑1 (MQ=255)                                                      aCGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTg                                                                                  <  5:293609/149‑1 (MQ=255)                                                       cGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTgg                                                                                 <  1:262937/149‑1 (MQ=255)                                                        ggCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGt                                                                                <  8:36509/149‑1 (MQ=255)                                                        ggCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGt                                                                                <  8:238634/149‑1 (MQ=255)                                                         gCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCAtt                                                                                                        >  2:186250/1‑124 (MQ=255)                                                         gCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCAtt                                                                                                        <  1:186250/124‑1 (MQ=255)                                                          cacaTCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTg                                                                               >  7:237348/1‑148 (MQ=255)                                                           acaTCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCGGGCAAGTCGGGTGATCATGCGCGCCCTGGGTGa                                                                              >  8:49808/1‑148 (MQ=255)                                                             aTCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCt                                                                                                                               >  2:9906/1‑97 (MQ=255)                                                             aTCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCt                                                                                                                               <  1:9906/97‑1 (MQ=255)                                                             aTCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAAt                                                                           <  1:168429/149‑1 (MQ=255)                                                             aTCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAAt                                                                           >  2:168429/1‑149 (MQ=255)                                                             aTCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAAt                                                                           <  8:289689/149‑1 (MQ=255)                                                                 gCCTGATGGCGCTGGGCGCCTGTAGCGGTGTCTGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGcc                                                                       <  6:228579/149‑1 (MQ=255)                                                                   cTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGt                                                                     <  8:237348/149‑1 (MQ=255)                                                                        ggCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGcac                                                                 >  5:210032/1‑148 (MQ=255)                                                                                gcgcCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATgc                                                                                           >  4:57295/1‑114 (MQ=255)                                                                                gcgcCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATgc                                                                                           <  3:57295/114‑1 (MQ=255)                                                                                            gtgtCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCt                                                           <  5:79934/134‑1 (MQ=255)                                                                                            gtgtCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCt                                                           >  6:79934/1‑134 (MQ=255)                                                                                                 ggTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGTTCGTTCAGCCg                                       >  1:281764/1‑149 (MQ=255)                                                                                                  gTCTTTTCCTCGGTTCGGCCAAGGCGATCGATATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGt                                      >  5:178106/1‑149 (MQ=255)                                                                                                  gTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGt                                      >  1:90576/1‑149 (MQ=255)                                                                                                   tCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGt                                                                                >  1:249328/1‑106 (MQ=255)                                                                                                   tCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGt                                                                                <  2:249328/106‑1 (MQ=255)                                                                                                        tCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGc                                 <  2:90576/148‑1 (MQ=255)                                                                                                        tCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGACCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGc                                 <  7:183712/148‑1 (MQ=255)                                                                                                         ccTCGGTTCGGCCAAGGCGATCGAGATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGcc                                <  6:178106/148‑1 (MQ=255)                                                                                                                   gCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCt                     >  6:266023/1‑149 (MQ=255)                                                                                                                     cAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAg                             >  5:239001/1‑139 (MQ=255)                                                                                                                     cAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAg                             <  6:239001/139‑1 (MQ=255)                                                                                                                         gCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGTTCGTTCAGCCGTTATGCCCAGGATTACCTCGGccc               <  2:281764/149‑1 (MQ=255)                                                                                                                            aTCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCAtt            >  8:39075/1‑149 (MQ=255)                                                                                                                               gAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGc         <  6:210032/149‑1 (MQ=255)                                                                                                                               atgaTGGCCGGGCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGc         >  6:232076/4‑149 (MQ=255)                                                                                                                                tgaTGGCCGGGCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCgg                  <  5:166591/137‑1 (MQ=255)                                                                                                                                tgaTGGCCGGGCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCgg                  >  6:166591/3‑139 (MQ=255)                                                                                                                                aaaTGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGCg        <  2:150934/149‑1 (MQ=255)                                                                                                                                  aTGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGCggg      <  5:266023/149‑1 (MQ=255)                                                                                                                                     gacggcccggccATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCTCCCTTGGAGAAATGGCCGTTCACAACCATGTCGCCGGCTCTTTCAGCCGTTATGACCAGGATTACATCGGACAATAGGCGGGCtt   >  7:206885/3‑149 (MQ=255)                                                                                                                                      ccgggccggcCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGCGGGCTTc  <  5:232076/149‑1 (MQ=255)                                                                                                                                      ccggcccggcCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGGCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGCGGGCTTc  <  8:101871/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                 TCTGCCCACGGCCAGGCCACTGAAGGCGGCCCGCTCAAGCGTGAACTGGGCGAACGGCACATCCGCCTGATGGCGCTGGGCGCCTGTATCGGTGTCGGTCTTTTCCTCGGTTCGGCCAAGGCGATCGAAATGGCCGGCCCGGCCATCATGTTGTCCTACATCATCGGCGGCCTGGCCATTCTGGTGATCATGCGCGCCCTTGGTGAAATGGCCGTGCACAACCCTGTCGCCGGCTCGTTCAGCCGTTATGCCCAGGATTACCTCGGCCCATTGGCGGGCTTC  >  NC_002947/1209410‑1209691

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