Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,953,402	Δ1 bp	6.6%	coding (173/645 nt)	PP_5218 →	DedA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,953,402	0	G	.	6.6%	193.2 / 6.8	61	coding (173/645 nt)	PP_5218	DedA family protein
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base . (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGCTACTTCGCCCTTTTTCTAGGCACCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTCTATGGTTTACGCACGGT  >  NC_002947/5953255‑5953549                                                                                                                                                     |                                                                                                                                                     tCGGCTACTTCGCCCTTTTTCTAGGCACCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTtgtg                                                                                                                                                      <  8:4410/148‑1 (MQ=255)         ttCGCCCTTTTTCTAGGCACCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCAt                                                                                                                                             >  8:61723/1‑149 (MQ=255)            gCCCTTTTTCTAGGCACCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATggg                                                                                                                                          >  8:155750/1‑149 (MQ=255)                        ggCACCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACgg                                                                                                                             >  6:301091/1‑149 (MQ=255)                           aCCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCAt                                                                                                                                             >  5:55465/1‑131 (MQ=255)                           aCCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCAt                                                                                                                                             <  6:55465/131‑1 (MQ=255)                           aCCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCATGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGT‑GTACTTCATGGGCCG‑ACGCCACGGgcg                                                                                                                          >  6:175155/1‑148 (MQ=255)                              ttcttTGAAGGCGAGAACATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTACGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACAAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAg                                                                                                                       >  6:277289/1‑149 (MQ=255)                                cttTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGa                                                                                                                     >  8:197433/1‑149 (MQ=255)                                   tGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATc                                                                                                                   <  2:69275/148‑1 (MQ=255)                                     aaGGCGAGACCATTCTGGTGCTTGCGGGATTTCATTGCGTTCGGCGGTTACATGGATATCAATCTGTTATGCATGGTGGCGTTATTCGGAAGCTACGCGGGCGACCAGTTGTGTTAATTCATGGGc                                                                                                                                         >  1:187823/1‑126 (MQ=255)                                     aaGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGc                                                                                                                                         <  2:187823/125‑1 (MQ=255)                                     aaGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCt                                                                                                                 <  1:237996/148‑1 (MQ=255)                                       ggCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgc                                                                                                              <  5:301091/149‑1 (MQ=255)                                             aCCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATgg                                                                                                                                           <  3:302371/115‑1 (MQ=255)                                             aCCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATgg                                                                                                                                           >  4:302371/1‑115 (MQ=255)                                             aCCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgcgcgc                                                                                                          >  5:296202/1‑147 (MQ=255)                                                aTTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGcc                                                                                                     >  7:81594/1‑149 (MQ=255)                                                   cTGGGGCTTGCGGGA‑TTCCTTTCGTTCCGCGGTTACATGGATTTCAATCTTGTATGGCTGGTGGCCTTCTTCGGCAGCTTCGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAATCcgc                                                                                                   <  4:284445/148‑1 (MQ=255)                                                      gTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCgcgc                                                                                                 <  7:61723/147‑1 (MQ=255)                                                        gCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTgg                                                                                              >  4:117306/1‑148 (MQ=255)                                                          ttGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGCCCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTggcag                                                                                           >  5:115758/1‑149 (MQ=255)                                                          ttGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgcgc                                                                                                            <  4:250390/132‑1 (MQ=255)                                                          ttGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgcgc                                                                                                            >  3:250390/1‑132 (MQ=255)                                                              gggA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgc                                                                                                              <  6:247790/126‑1 (MQ=255)                                                              gggA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgc                                                                                                              >  5:247790/1‑126 (MQ=255)                                                              gggA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCaa                                                                                       <  8:81594/149‑1 (MQ=255)                                                                 a‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATgg                                                                                    <  7:155750/149‑1 (MQ=255)                                                                      cTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGT‑GTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGAc                                                                               >  5:44559/1‑149 (MQ=255)                                                                         gCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGAACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTgg                                                                                              <  7:203267/132‑1 (MQ=255)                                                                         gCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGAACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTgg                                                                                              >  8:203267/1‑132 (MQ=255)                                                                           gTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTAGGCGGTCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgcgc                                                                                                            <  5:81745/116‑1 (MQ=255)                                                                           gTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgcgc                                                                                                            >  7:17522/1‑116 (MQ=255)                                                                           gTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgcgc                                                                                                            <  8:17522/116‑1 (MQ=255)                                                                           gTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTAAGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGgcgc                                                                                                            >  6:81745/1‑116 (MQ=255)                                                                            ttCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGgc                                                                          <  6:115758/149‑1 (MQ=255)                                                                              ccGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGT‑GTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGgcgc                                                                        <  6:44559/148‑1 (MQ=255)                                                                               cgcgGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGgcgc                                                                        <  3:117306/148‑1 (MQ=255)                                                                               cgcgGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCt                                                                       >  4:61683/1‑149 (MQ=255)                                                                               cgcgGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGT‑GTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCt                                                                       <  5:175155/148‑1 (MQ=255)                                                                                        aTGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACAt                                                              >  2:24186/1‑149 (MQ=255)                                                                                             tatCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATccgcc                                                         >  8:104548/1‑149 (MQ=255)                                                                                              atCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATAGGTGACCGGGCGCTGGACCACATccgcc                                                         <  2:149501/148‑1 (MQ=255)                                                                                                       gTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCAccc                                                  <  6:296202/146‑1 (MQ=255)                                                                                                       gTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCAcc                                                   <  1:24186/145‑1 (MQ=255)                                                                                                           gCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATc                                           >  5:64967/1‑149 (MQ=255)                                                                                                                   gCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCt                                   <  7:197433/149‑1 (MQ=255)                                                                                                                   cccTTCATCCGCAGCTACGCGGGCGACCAGAAGTGGAACAACATGCGCCG‑ACGCCACGCGCGCAGGATCCTGGCGCGCAAGCCGCGCACGCAGCCAATGGGTCACCCGGCGCTCGACCACATCCGCCTCCACCCCGATATCTGGGTCCt                                   >  6:94158/2‑149 (MQ=255)                                                                                                                      ttcttcGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAg                                >  6:99959/1‑149 (MQ=255)                                                                                                                        cttcCGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGc                               <  5:277677/148‑1 (MQ=255)                                                                                                                         ttcGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCtt                             >  8:199320/1‑149 (MQ=255)                                                                                                                              cAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCt                        <  7:104548/149‑1 (MQ=255)                                                                                                                                  tACGCGGGCGACCAGTTGTGGTACTTCATGGGCCGAACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTg                     >  1:205751/1‑149 (MQ=255)                                                                                                                                  tACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGt                    >  1:10853/1‑149 (MQ=255)                                                                                                                                   aCGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTc                   >  5:57581/1‑149 (MQ=255)                                                                                                                                    cgcgGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTc                   >  1:301090/1‑148 (MQ=255)                                                                                                                                     gcgGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTg                                  <  7:109131/132‑1 (MQ=255)                                                                                                                                     gcgGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTg                                  >  8:109131/1‑132 (MQ=255)                                                                                                                                      cgGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTCTAt                <  6:198771/149‑1 (MQ=255)                                                                                                                                          cGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTCTATGGtt            <  5:277289/149‑1 (MQ=255)                                                                                                                                            aCCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGTGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTCTATGGTTTa          <  2:301090/149‑1 (MQ=255)                                                                                                                                             ccAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTCTATGGTTTAc         >  4:121770/1‑149 (MQ=255)                                                                                                                                                    tgGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTCTATGGTTTACGCACGGt  >  2:68431/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TCGGCTACTTCGCCCTTTTTCTAGGCACCTTCTTTGAAGGCGAGACCATTCTGGTGCTTGCGGGA‑TTCCTTGCGTTCCGCGGTTACATGGATATCAATCTGGTATGCCTGGTGGCCTTCTTCGGCAGCTACGCGGGCGACCAGTTGTGGTACTTCATGGGCCG‑ACGCCACGGGCGCAGGATCCTGGCGCGCAAGCCGCGCTGGCAGGCAATGGGTGACCGGGCGCTGGACCACATCCGCCGCCACCCCGATATCTGGGTGCTGAGCTTCCGCTTTGTCTATGGTTTACGCACGGT  >  NC_002947/5953255‑5953549

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