Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,644,951	A→G	8.1%	intergenic (‑46/+20)	PP_3211 ← / ← PP_3212	ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,644,951	0	A	G	8.1%	176.2 / 4.3	63	intergenic (‑46/+20)	PP_3211/PP_3212	ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base G (3/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAG  >  NC_002947/3644805‑3645097                                                                                                                                                   |                                                                                                                                                   gTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCAcc                                                                                                                                                  >  3:159666/1‑149 (MQ=255)     tttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTTCATCAACATAGCTCTCTCCCCGTGGGTCATTAGCGGTTGTGTATGCCCCGGCGCACCggg                                                                                                                                               <  2:253137/148‑1 (MQ=255)      ttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGc                                                                                                                                             >  1:116864/1‑149 (MQ=255)         cACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGtt                                                                                                                                          >  2:182751/1‑149 (MQ=255)            cTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTggg                                                                                                                                       <  3:249745/149‑1 (MQ=255)             tGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGc                                                                                                                                      <  4:159666/149‑1 (MQ=255)                    cGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTc                                                                                                                               <  1:252163/149‑1 (MQ=255)                     gAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGGCGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAgtgt                                                                                                                                >  2:178426/1‑147 (MQ=255)                     gAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGGCGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCa                                                                                                                              <  1:178426/149‑1 (MQ=255)                             gggTGTCGAGAGCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCTGCCTTGGATAGTGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                      <  2:116864/149‑1 (MQ=255)                             gggTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                      <  1:173643/149‑1 (MQ=255)                             gggTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                      <  1:182751/149‑1 (MQ=255)                                         cGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGc                                                                                                          >  8:163931/1‑149 (MQ=255)                                            tGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGctct                                                                                                       >  8:7551/1‑149 (MQ=255)                                              gtgtTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCaa                                                                                                                                    >  2:30005/1‑118 (MQ=255)                                              gtgtTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCaa                                                                                                                                    <  1:30005/118‑1 (MQ=255)                                                 ttCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCa                                                                                                  <  8:112391/149‑1 (MQ=255)                                                     gCAGGGCCTGCGGCTGGACATCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATggg                                                                                              <  2:47228/149‑1 (MQ=255)                                                              gCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAGCACTTCGCTCTCGCCATGGGCTTCAAg                                                                                       >  7:138824/1‑147 (MQ=255)                                                              gCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAg                                                                                       >  3:35901/1‑147 (MQ=255)                                                               cGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAg                                                                                       >  4:261594/1‑146 (MQ=255)                                                                  cTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTACGCTCTCGCCATGGGCTTCAAGCAGAcc                                                                                 <  2:36982/149‑1 (MQ=255)                                                                     gACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCag                                                                              >  2:174193/1‑149 (MQ=255)                                                                      aCGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCagt                                                                             <  7:163931/149‑1 (MQ=255)                                                                           ggCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGtt                                                                         <  7:7551/148‑1 (MQ=255)                                                                           ggCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTg                                                                        >  7:136777/1‑149 (MQ=255)                                                                             ccTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGcc                                                                      >  4:223695/1‑149 (MQ=255)                                                                                   aTAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGa                                                                <  4:35901/149‑1 (MQ=255)                                                                                      gggTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATaa                                                             <  8:136777/149‑1 (MQ=255)                                                                                        gTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATaac                                                            >  5:90800/1‑148 (MQ=255)                                                                                              aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACg                                                                  >  6:264077/1‑136 (MQ=255)                                                                                              aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACg                                                                  <  5:264077/136‑1 (MQ=255)                                                                                                cAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGAt                                                               <  8:88193/137‑1 (MQ=255)                                                                                                cAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGAt                                                               >  7:88193/1‑137 (MQ=255)                                                                                                   cATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGc                                                                       <  1:149811/126‑1 (MQ=255)                                                                                                   cATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGc                                                                       >  2:149811/1‑126 (MQ=255)                                                                                                        cgcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGc                                                                                                    >  4:265152/1‑92 (MQ=255)                                                                                                        cgcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGc                                                                                                    <  3:265152/92‑1 (MQ=255)                                                                                                        cgcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGt                                              <  3:185906/146‑1 (MQ=255)                                                                                                        cgcTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGt                                              >  4:185906/1‑146 (MQ=255)                                                                                                           tctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATg                                        >  4:97840/1‑149 (MQ=255)                                                                                                            ctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATgg                                       >  6:71096/1‑149 (MQ=255)                                                                                                               cccGTGGGGCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAt                                    <  6:84892/149‑1 (MQ=255)                                                                                                                ccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtg                                   >  8:161085/1‑149 (MQ=255)                                                                                                                                  ggtgTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTgcgg                 <  7:161085/149‑1 (MQ=255)                                                                                                                                     gtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtgtg                                 <  7:139988/130‑1 (MQ=255)                                                                                                                                     gtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtgtg                                 >  8:139988/1‑130 (MQ=255)                                                                                                                                     caacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGAc                                                                   >  4:176611/5‑96 (MQ=255)                                                                                                                                     caacgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGAc                                                                   <  3:176611/92‑1 (MQ=255)                                                                                                                                      tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAGCACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              <  8:138824/148‑1 (MQ=255)                                                                                                                                      tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGa                                          <  8:2624/120‑1 (MQ=255)                                                                                                                                      tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGa                                          >  7:2624/1‑120 (MQ=255)                                                                                                                                      tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              <  6:90800/148‑1 (MQ=255)                                                                                                                                      tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              <  3:97840/148‑1 (MQ=255)                                                                                                                                      tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              <  3:261594/148‑1 (MQ=255)                                                                                                                                       acgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgt            >  1:88364/3‑149 (MQ=255)                                                                                                                                        tGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgt            <  5:71096/148‑1 (MQ=255)                                                                                                                                         gCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCtt                                                <  3:17940/111‑1 (MQ=255)                                                                                                                                         gCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCtt                                                >  4:17940/1‑111 (MQ=255)                                                                                                                                          ccccGGCGCTCCGGGGCGTGGTGCAAGTGTCAGTGGCGGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCTGgtgt          <  8:187902/148‑1 (MQ=255)                                                                                                                                          ccccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc         <  7:196933/149‑1 (MQ=255)                                                                                                                                           cccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc         >  3:105571/1‑148 (MQ=255)                                                                                                                                                 gcACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAg  <  2:81915/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAG  >  NC_002947/3644805‑3645097

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