Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,237,053	A→C	11.0%	intergenic (+105/+23)	uvrB → / ← PP_1975	uvrABC excinuclease protein B/EmrB/QacA family drug resistance transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,237,053	0	A	C	11.0%	146.7 / 12.1	63	intergenic (+105/+23)	uvrB/PP_1975	uvrABC excinuclease protein B/EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base C (4/3);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

GAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGAAGGCATCGTTGTAAGCCATGATCAGCGCCTG  >  NC_002947/2236910‑2237201                                                                                                                                                |                                                                                                                                                     gagaTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGAc                                                                                                                                                     >  2:1275616/1‑145 (MQ=255) gagaTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGAc                                                                                                                                                     <  1:1275616/145‑1 (MQ=255) gagaTTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCACGCCCGGGTCCAGAGGGACCGCat                                                                                                                                                <  2:315610/150‑1 (MQ=37)     ttGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCatctga                                                                                                                                            <  1:1246517/150‑1 (MQ=255)       gcgcAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCACGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCatctgatc                                                                                                                                          >  1:938104/1‑150 (MQ=255)        cgcAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCatctgatct                                                                                                                                         <  1:633486/150‑1 (MQ=255)        cgcAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCatctgatct                                                                                                                                         <  1:118436/150‑1 (MQ=255)           aGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTa                                                                                                                                      <  2:938104/150‑1 (MQ=255)                   gAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTTTTCGCGGGCATGCCCGCTCCCACAGCGACCGCATCTGATCTGTAATAGCCTc                                                                                                                              >  1:89237/1‑150 (MQ=255)                      cGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAAt                                                                                                                           <  2:1128159/150‑1 (MQ=255)                      cGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAAt                                                                                                                           <  1:840974/150‑1 (MQ=255)                         ttGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAAt                                                                                                                           >  2:25603/1‑147 (MQ=255)                         ttGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAAt                                                                                                                           <  1:25603/147‑1 (MQ=255)                              cagcagCTGAAGCCGGGCGCGGGCCCTGTGGGAGCCGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTcc                                                                                                                   <  1:1056687/150‑1 (MQ=255)                                gcagcTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGCCCGCATCTGATCTGTAATAGCCTCAATGGGCTTcccc                                                                                                                 >  1:1030982/1‑150 (MQ=255)                                       aaGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTc                                                                                                          <  2:1030982/150‑1 (MQ=255)                                               cgcgGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTc                                                                                                  <  2:354778/150‑1 (MQ=255)                                                 cgGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGTCCCTCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCg                                                                                                <  1:655839/150‑1 (MQ=255)                                                      ccTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAAc                                                                                           >  1:1150779/1‑150 (MQ=255)                                                       cTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAAt                                                                                                                           <  1:1224341/117‑1 (MQ=255)                                                       cTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAAt                                                                                                                           >  2:1224341/1‑117 (MQ=255)                                                          tgGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCt                                                                                                                               >  1:465163/1‑110 (MQ=255)                                                          tgGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCt                                                                                                                               <  2:465163/110‑1 (MQ=255)                                                            ggAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCgg                                                                                     >  2:661653/1‑150 (MQ=255)                                                              aGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGc                                                                                                              <  2:1016164/123‑1 (MQ=255)                                                              aGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGc                                                                                                              >  1:1016164/1‑123 (MQ=255)                                                              aGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCg                                                                                                >  1:1146441/1‑137 (MQ=255)                                                              aGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCg                                                                                                <  2:1146441/137‑1 (MQ=255)                                                               gcgggcatgcccgcgaacaccggcgcagccggtgccagccaccgcgttgaattcttcgcgggcacGCCCGCTCCCACAGGGCCCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTc                                                                                  >  1:813248‑M2/66‑150 (MQ=255)                                                               gCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTc                                                                                  <  1:661653/150‑1 (MQ=255)                                                                  ggcatgcccgcgaacaccggcgcagccggtgccagccaccgcgttgaattcttcgcgggcacGCCCGCTCCCACAGGGCCCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTcagc                                                                               >  2:1155185‑M2/63‑150 (MQ=255)                                                                   gcatgcccgcgaacaccggcgcagccggtgccagccaccgcgttgaattcttcgcgggcacGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCg                                                                                                >  2:1136081‑M2/62‑132 (MQ=255)                                                                   gcatgcccgcgaacaccggcgcagccggtgccagccaccgcgttgaattcttcgcgggcacGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCg                                                                                                <  1:1136081‑M2/71‑1 (MQ=255)                                                                                 aaTTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGGTCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATg                                                                >  2:590861/1‑150 (MQ=255)                                                                                   ttCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCa                                                                                                                <  1:1135560/100‑1 (MQ=255)                                                                                   ttCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCa                                                                                                                >  2:1135560/1‑100 (MQ=255)                                                                                   ttCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCt                                                              >  1:1218238/1‑150 (MQ=255)                                                                                           ggTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACTCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCa                                                      >  2:248620/1‑150 (MQ=255)                                                                                           ccggtgccagccaccgcgttgaattcttcgcgggcacGCCCGCTCCCACAGGGCCCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCa                                                      <  2:813248‑M2/113‑1 (MQ=255)                                                                                             tggctggcACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAAc                                                    >  1:285058/1‑150 (MQ=255)                                                                                             ggtgccagccaccgcgttgaattcttcgcgggcacGCCCGCTCCCACAGGGCCCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAAc                                                    <  1:1155185‑M2/115‑1 (MQ=255)                                                                                                      ccGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTccc                                                                                                                  >  1:810547/1‑79 (MQ=255)                                                                                                      ccGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTccc                                                                                                                  <  2:810547/79‑1 (MQ=255)                                                                                                      ccGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCg                                           <  1:270252/150‑1 (MQ=255)                                                                                                         gCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACa                                        >  1:162909/1‑150 (MQ=255)                                                                                                             cgcCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGt                                    <  2:89237/150‑1 (MQ=255)                                                                                                             cgcCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGt                                    <  1:512381/150‑1 (MQ=255)                                                                                                                 tgtgggagcgGGCATGCCCGCGAACACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGaa                                >  2:1195498/8‑150 (MQ=255)                                                                                                                 ggTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGaa                                <  2:1150779/150‑1 (MQ=255)                                                                                                                     ttCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCaa                                                                                            <  1:119445/86‑1 (MQ=255)                                                                                                                     ttCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCaa                                                                                            >  2:119445/1‑86 (MQ=255)                                                                                                                           ggCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCgg                                                                                     <  2:130625/87‑1 (MQ=255)                                                                                                                           ggCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCgg                                                                                     >  1:130625/1‑87 (MQ=255)                                                                                                                              aTGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACTCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGAAGGCATCGTTGTaa                   <  1:248620/150‑1 (MQ=255)                                                                                                                              aTGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGGGGAAGGCATCGTTGTaa                   <  2:285058/150‑1 (MQ=255)                                                                                                                                  ccGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGAAGGCATCGTTGTAAGCCa               >  1:284148/1‑150 (MQ=255)                                                                                                                                  ccGCTCCCACAGGGACCGCATCTGATCTATAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGAAGGCATCGTTGTAAGCCa               <  2:162909/150‑1 (MQ=255)                                                                                                                                        ccACAGGGCCCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGCTCCCGACAAAg                                     <  2:119920/122‑1 (MQ=255)                                                                                                                                        ccACAGGGCCCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAg                                     >  1:119920/1‑122 (MQ=255)                                                                                                                                          acaGGGACCGCATCTGATCTGTAATGGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGAAGGCATCGTTGTAAGCCATGATCAgc       <  2:564366/150‑1 (MQ=255)                                                                                                                                           caGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAcc                                                          <  2:266036/98‑1 (MQ=255)                                                                                                                                           caGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCAcc                                                          >  1:266036/1‑98 (MQ=255)                                                                                                                                               gACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGAAGGCATCGTTGTAAGCCATGATCAGCGCctg  >  2:308983/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                     GAGATTGCGCAGCTGCGGGAGCGGTTGATCAGCAGCTGAAGCCGGGCGCGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATTCAACGCGGTGGCTGGCACCGGCTGCGCCGGTGTTCGCGGGCATGCCCGCTCCCACAGGGACCGCATCTGATCTGTAATAGCCTCAATGGGCTTCCCCAGCCTTCAACCCTTCCGGCAACTTGCGGGTCAGCAGCACCGCCACCATGCTCACCGCCAACCCGATCCCGACAAAGTGGAAGGCATCGTTGTAAGCCATGATCAGCGCCTG  >  NC_002947/2236910‑2237201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑