Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	545,031	C→G	7.9%	P1091A (CCA→GCA)	rpoC →	DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	545,031	0	C	G	7.9%	194.0 / 6.2	63	P1091A (CCA→GCA)	rpoC	DNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (2/3);  total (31/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 = 9.57e-01

ACCATCAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGTGACATTACCGGTGGTCTGCCACG  >  NC_002947/544893‑545179                                                                                                                                           |                                                                                                                                                      gccATCAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGa                                                                                                                                             <  7:74057/148‑1 (MQ=255) aCCATCAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGa                                                                                                                                             <  7:248391/149‑1 (MQ=255) aCCATCAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGa                                                                                                                                             <  7:251716/149‑1 (MQ=255)   cATCAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCAATGCTGCAGGCAAGGATCTGTACCTGGCAGGTACAGatc                                                                                                                                           <  8:5385/149‑3 (MQ=255)      cAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAg                                                                                                                                                    >  6:107313/1‑137 (MQ=255)      cAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAg                                                                                                                                                    <  5:107313/137‑1 (MQ=255)          cGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGc                                                                                                                                    <  5:45003/149‑1 (MQ=255)           gTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGcc                                                                                                                                   <  7:73423/149‑1 (MQ=255)            tCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGcc                                                                                                                                   >  2:241780/1‑148 (MQ=255)               gACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAg                                                                                                                                >  3:59650/1‑148 (MQ=255)                aCCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTa                                                                                                                              <  4:169100/149‑1 (MQ=255)                        aCTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGc                                                                                                                      >  5:100467/1‑149 (MQ=255)                            aCCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTTCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCgg                                                                                                                   <  8:146260/148‑1 (MQ=255)                             ccGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTg                                                                                                                                              <  1:146211/120‑1 (MQ=255)                             ccGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTg                                                                                                                                              >  2:146211/1‑120 (MQ=255)                                gCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGACGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTAATTCCTGCCGGCCAAc                                                                                                              >  7:33202/1‑149 (MQ=255)                                       cAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGGCAGGTAca                                                                                                                                               >  5:174221/1‑108 (MQ=255)                                       cAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGGCAGGTAca                                                                                                                                               <  6:174221/109‑2 (MQ=255)                                       cAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTg                                                                                                       >  2:77745/1‑149 (MQ=255)                                       cAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTg                                                                                                       >  6:154486/1‑149 (MQ=255)                                        aaCATTGAAGTACTGGACGTCAAGGATCGCGCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGt                                                                                                      <  1:241780/149‑1 (MQ=255)                                        aaCATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGt                                                                                                      <  4:134071/149‑1 (MQ=255)                                         aCATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGccc                                                                                                          <  5:185855/144‑1 (MQ=255)                                         aCATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGccc                                                                                                          >  6:185855/1‑144 (MQ=255)                                               aaGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGg‑tgt                                                                                               <  1:77745/148‑1 (MQ=255)                                                  tACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCt                                                                                           >  7:225492/1‑149 (MQ=255)                                                       gACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACgg                                                                                      >  3:247038/1‑149 (MQ=255)                                                          gTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGc                                                                                   >  7:11176/1‑149 (MQ=255)                                                          gTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGc                                                                                   >  2:191629/1‑149 (MQ=255)                                                             aaGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGTCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCt                                                                                  >  8:100720/1‑147 (MQ=255)                                                             aaGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTc                                                                                 <  5:154486/148‑1 (MQ=255)                                                             aaGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCa                                                                                <  4:247038/149‑1 (MQ=255)                                                                gATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCt                                                                                  >  5:221170/1‑144 (MQ=255)                                                                    gCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCgg                                                                          <  6:221170/148‑1 (MQ=255)                                                                      ccTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGGCAg                                                                                                                                                    <  2:180039/73‑1 (MQ=255)                                                                      ccTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGGCAg                                                                                                                                                    >  1:180039/1‑73 (MQ=255)                                                                          ccGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTACCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGt                                                                   >  1:169006/1‑149 (MQ=255)                                                                             cAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGAc                                                                <  6:100467/149‑1 (MQ=255)                                                                                gCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGt                                                              >  4:205550/1‑148 (MQ=255)                                                                                  aaGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTAt                                                           >  8:170714/1‑149 (MQ=255)                                                                                    ggAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGc                                                                                   >  8:5821/1‑123 (MQ=255)                                                                                    ggAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGc                                                                                   <  7:5821/123‑1 (MQ=255)                                                                                    ggAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTg                                                                     >  1:182023/1‑137 (MQ=255)                                                                                    ggAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTg                                                                     <  2:182023/137‑1 (MQ=255)                                                                                       aaTCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGAGCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATcgcgc                                                      <  4:99230/149‑1 (MQ=255)                                                                                        aTCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGGCTGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATcgcgc                                                      <  1:191629/149‑1 (MQ=255)                                                                                              ccgGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATccc                                               <  8:225492/149‑1 (MQ=255)                                                                                                 gCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCATGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGTTCGGTGTTGGTGACGTTATCGCGCGTATCCCGc                                             >  6:87036/1‑148 (MQ=255)                                                                                                    aTCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCATGCCGGCCAACTCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGTTCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAaga                                         <  5:87036/149‑1 (MQ=255)                                                                                                              tCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCg                                                                                                                    >  3:60024/1‑65 (MQ=255)                                                                                                              tCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCg                                                                                                                    <  4:60024/65‑1 (MQ=255)                                                                                                                gATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGAc                             <  3:205550/149‑1 (MQ=255)                                                                                                                 aTGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGaa                                >  1:29826/1‑145 (MQ=255)                                                                                                                 aTGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGaa                                <  2:29826/145‑1 (MQ=255)                                                                                                                     tgcAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGt                         >  2:50192/1‑148 (MQ=255)                                                                                                                      gcAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGc                                             >  6:13928/1‑127 (MQ=255)                                                                                                                      gcAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGc                                             <  5:13928/127‑1 (MQ=255)                                                                                                                           cAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGTGACATTa                  >  8:147304/1‑149 (MQ=255)                                                                                                                            aaGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGCTCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGTGACATTAc                 >  3:95843/1‑149 (MQ=255)                                                                                                                                    gTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGCCCCGTGACATTACCGGTGGTc         >  7:156793/1‑149 (MQ=255)                                                                                                                                     tACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCTGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACGCGTGACATTACCGGTGGTc         <  8:33202/148‑1 (MQ=255)                                                                                                                                     tACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGTGACATTACCggt             <  7:170714/144‑1 (MQ=255)                                                                                                                                          gCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGTGACATTACCGGTGGTCTGCCAc   >  6:19700/1‑149 (MQ=255)                                                                                                                                           ccAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCg                          <  6:182039/125‑1 (MQ=255)                                                                                                                                           ccAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCg                          >  5:182039/1‑125 (MQ=255)                                                                                                                                           ccAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGTGACATTACCGGTGGTCTGCCAcg  <  8:11176/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                      ACCATCAAGCGTCAGACCGACGAACTGACCGGCTTGACCAACATTGAAGTACTGGACGTCAAGGATCGCCCTGCCGCAGGCAAGGAAATCCGTCCGGCAATCAAGATGGTCGATGCTGCAGGCAAGGATCTGTACCTGCCAGGTACTGACGTACCTGCCCAGTACTTCCTGCCGGCCAACGCCCTTGTCGG‑TGTGGCTGACGGTGCTCAGATCGGTGTTGGTGACGTTATCGCGCGTATCCCGCAAGAGACGTCGAAGACCCGTGACATTACCGGTGGTCTGCCACG  >  NC_002947/544893‑545179

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