Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,442,673	C→G	12.1%	intergenic (‑99/+39)	nicP‑II ← / ← nicT	porin‑like protein/metabolite transport protein NicT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,442,673	0	C	G	12.1%	160.3 / 14.6	58	intergenic (‑99/+39)	nicP‑II/nicT	porin‑like protein/metabolite transport protein NicT
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base G (4/3);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

CGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTG  >  NC_002947/4442526‑4442820                                                                                                                                                    |                                                                                                                                                     cgGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGACACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCAGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCgg                                                                                                                                                    <  1:293870/150‑1 (MQ=255)  ggCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTCTTCGACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGc                                                                                                                                                   <  1:118423/150‑1 (MQ=255)      aGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGc                                                                                                                                               <  2:201754/150‑1 (MQ=255)       ggCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcc                                                                                                                                              <  1:1043619/150‑1 (MQ=255)            gCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGGCTGCATGCGGGCCCCGGCCCGcctgccc                                                                                                                                         >  2:1753337/1‑150 (MQ=255)             cAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGcctgccct                                                                                                                                        >  1:1568405/1‑150 (MQ=255)                   aagaaAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGtt                                                                                                                                  <  1:229093/150‑1 (MQ=255)                       aaaaCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGt                                                                                                                             >  2:1821908/1‑150 (MQ=255)                          aCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGaa                                                                                                                          <  1:764720/150‑1 (MQ=255)                          aCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGaa                                                                                                                          <  2:1737140/150‑1 (MQ=255)                          aCGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGaa                                                                                                                          <  1:1360972/150‑1 (MQ=255)                                     tGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGAGATGCATCGAGAGAGGTCGGCAGGCGGG‑CCCGGCAGCACTGGACTCCCGTTGGCGTG‑AGCGCTGGGCCGc                                                                                                              >  1:1468073/1‑150 (MQ=255)                                        aaCCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTa                                                                                                            >  1:540931/1‑150 (MQ=255)                                         aCCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGACGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCAGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGATaa                                                                                                           >  2:432726/1‑149 (MQ=255)                                           cTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCg                                                                                                         <  2:1236982/150‑1 (MQ=255)                                            ttGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCgg                                                                                                        <  2:540931/150‑1 (MQ=255)                                            gtgcgtggggCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGGGAAGCGCTGGGCCGCTAGCgg                                                                                                        >  2:968087/8‑150 (MQ=255)                                                aTAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGa                                                                                                    >  2:639978/1‑150 (MQ=255)                                                 tAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGAc                                                                                                   <  1:280626/150‑1 (MQ=255)                                                            cTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGCGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCg                                                                                        >  2:103832/1‑150 (MQ=255)                                                            cTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCg                                                                                        >  1:629573/1‑150 (MQ=255)                                                                       cAACAGCAACCCGGCGCGGCGGCACTGCATGGCTACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATg                                                                             >  2:159944/1‑150 (MQ=255)                                                                           aGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCa                                                                         >  2:83036/1‑150 (MQ=255)                                                                           aGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCa                                                                         >  1:606925/1‑150 (MQ=255)                                                                           aGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCa                                                                         <  2:1568405/150‑1 (MQ=255)                                                                                 ccGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAagc                                                                   >  1:506015/1‑150 (MQ=255)                                                                                       cggcggCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                             <  1:1753337/150‑1 (MQ=255)                                                                                       cggcggCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                             <  1:1420458/150‑1 (MQ=255)                                                                                              aCTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAg                                                      >  1:46191/1‑150 (MQ=255)                                                                                                       ccACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGc                                             >  1:1684866/1‑150 (MQ=255)                                                                                                                 gCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCcac                                   <  1:159944/150‑1 (MQ=255)                                                                                                                     gTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGc                                                                                          <  1:1673232/91‑1 (MQ=255)                                                                                                                     gTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGc                                                                                          >  2:1673232/1‑91 (MQ=255)                                                                                                                      ttGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                              >  2:741627/1‑102 (MQ=255)                                                                                                                      ttGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGAt                                                                              <  1:741627/102‑1 (MQ=255)                                                                                                                       tGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgca                             <  2:554866/150‑1 (MQ=255)                                                                                                                          aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgcagcg                          >  1:1760710/1‑150 (MQ=255)                                                                                                                          aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgca                             <  2:370934/147‑1 (MQ=255)                                                                                                                          aTCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCAcagcgca                             >  1:370934/1‑147 (MQ=255)                                                                                                                           tCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAg                                                      <  1:746201/121‑1 (MQ=255)                                                                                                                           tCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAg                                                      >  2:746201/1‑121 (MQ=255)                                                                                                                            cGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTCATGCACTTGGCCGGCAAGGCGATGACCAGGAAGCAGGTGAGGATCATGCAGCCGGCAAACACCCACAGCGCAGCGTg                        >  1:1562063/1‑150 (MQ=255)                                                                                                                                 gagGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAgg                   <  2:222695/150‑1 (MQ=255)                                                                                                                                    gTCGGCAGGCGGGCCCCGGCCCGCCAGTTGTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGc                >  1:1204519/1‑150 (MQ=255)                                                                                                                                     tCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCACGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGcc               >  1:960212/1‑150 (MQ=255)                                                                                                                                       ggCAGGCGGGCCGGGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGaa                                                                     <  2:944478/94‑1 (MQ=255)                                                                                                                                       ggCAGGCGGGCCGGGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGaa                                                                     >  1:944478/1‑94 (MQ=255)                                                                                                                                       ggCAGGCGGGCCGGGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGaa                                                                     >  1:945381/1‑94 (MQ=255)                                                                                                                                       ggCAGGCGGGCCGGGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                             <  1:954534/102‑1 (MQ=255)                                                                                                                                       ggCAGGCGGGCCGGGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGa                                                             >  2:954534/1‑102 (MQ=255)                                                                                                                                       ggCAGGCGGGCCGGGGCCCGACTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGaa                                                                     <  2:945381/94‑1 (MQ=255)                                                                                                                                          aGGCGGGCCGGGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa          >  1:6397/1‑150 (MQ=255)                                                                                                                                          aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGTCAAGGCGAt                                                                              >  2:474271/1‑82 (MQ=255)                                                                                                                                          aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGTCAAGGCGAt                                                                              <  1:474271/82‑1 (MQ=255)                                                                                                                                          aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGCAGCACCCGAGGATCAGGCAGCAGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa          >  1:1454364/1‑150 (MQ=255)                                                                                                                                          aGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGa          >  2:895847/1‑150 (MQ=255)                                                                                                                                            gCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAgg        <  1:1821908/150‑1 (MQ=255)                                                                                                                                                gCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCgg    <  2:629573/150‑1 (MQ=255)                                                                                                                                                 ccccGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGt   <  1:639978/150‑1 (MQ=255)                                                                                                                                                  cccGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCgg    <  2:542829/148‑1 (MQ=255)                                                                                                                                                  cccGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCgg    >  1:542829/1‑148 (MQ=255)                                                                                                                                                  cccGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTg  >  1:123043/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGGCCAGGCGTGCAGACGAAGAAAAACGCGGCGCAGTGAAACCTTGCATAGGGCTACTCCTGGATTTGGGCAACAGCAACCCGGCGCGGCGGCACTGCATGGCCACCGGCGTGCGGGTTGCATCGAGAGAGGTCGGCAGGCGGGCCCCGGCCCGCCTGCCCTGCCGTT‑GCGTGAAGCGCTGGGCCGCTAGCGGTTGACCAGCTTGGCCGGCAAGGCGATGACCAGGAAGCAGCTGAGGATCAGGCAGCCGGCAAACACCCACAGCGCAGCGTGCGAGGTGCCGGTGAGGTCGGTG  >  NC_002947/4442526‑4442820

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