Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,188,586	G→A	6.3%	T104M (ACG→ATG)	PP_4568 ←	metabolite efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,188,586	0	G	A	6.3%	186.6 / 2.2	63	T104M (ACG→ATG)	PP_4568	metabolite efflux pump
Reads supporting (aligned to +/- strand):  ref base G (26/33);  new base A (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CCGCCAGCGATACCCCGCCCAGGTTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAGC‑GCAGCCCTAGCAACGGGCCGATCA  >  NC_002947/5188440‑5188733                                                                                                                                                   |                                                                                                                                                        ccGCCAGCGATACCCCTCCCAGGGTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTg                                                                                                                                                       <  6:250544/149‑1 (MQ=255)    ccAGCGATACCCCGCCCAGGTTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGcc                                                                                                                                                    <  4:102684/149‑1 (MQ=255)     cAGCGATACCCCGCCCAGGTTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCa                                                                                                                                                   <  7:45490/149‑1 (MQ=255)        cGATACCCCGCCCATGTTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGa                                                                                                                                                <  5:310098/149‑1 (MQ=255)             ccccGCCCAGGTTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGGGCCCATGAC‑GCCg                                                                                                                                          >  5:25951/1‑149 (MQ=255)              cccGCCCAGGTTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCgg                                                                                                                                         <  3:72007/149‑1 (MQ=255)                          gATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATggg                                                                                                                            >  7:245348/1‑149 (MQ=255)                          gATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATgg                                                                                                                             <  6:308738/148‑1 (MQ=255)                          gATGCCCTCGTAGACCACCCTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATgg                                                                                                                             <  8:306066/148‑1 (MQ=255)                            tGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTa                                                                                                                          <  1:215926/149‑1 (MQ=255)                              ccatGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGTCTCCAGGAGCGCTGGAGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATAACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑t                                                                                                                        >  5:117124/4‑148 (MQ=255)                                cTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑Tc                                                                                                                       >  2:85938/1‑147 (MQ=255)                                        ccaccaTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAAc                                                                                                             <  6:203644/149‑1 (MQ=255)                                         caccaTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGc                                                                                                                                                     >  4:174840/1‑111 (MQ=255)                                         caccaTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGc                                                                                                                                                     <  3:174840/111‑1 (MQ=255)                                           ccaTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGTTCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCTACAGa                                                                                                          <  1:163886/149‑1 (MQ=255)                                           ccaTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCaa                                                                                                              >  5:11191/1‑145 (MQ=255)                                             aTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAg                                                                                                                    >  6:44224/1‑137 (MQ=255)                                             aTGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAg                                                                                                                    <  5:44224/137‑1 (MQ=255)                                                aTCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATagaa                                                                                                     <  3:144468/149‑1 (MQ=255)                                                 tCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACAt                                                                                                                               <  4:241277/123‑1 (MQ=255)                                                 tCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACAt                                                                                                                               >  3:241277/1‑123 (MQ=255)                                                 tCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATag                                                                                                       >  1:152035/1‑146 (MQ=255)                                                   caGCCGCCCAACCCCAAACCCAGGCCGAAAAAACGGGTGCAGGCGGGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATagaa                                                                                                     <  6:25951/145‑1 (MQ=255)                                                     gccgccCACCCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GGGACATGGGTA‑TCCAGCGCCAACAGATAGAAGATCt                                                                                                >  1:227306/1‑149 (MQ=255)                                                       cgccCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGa                                                                                              <  7:299450/149‑1 (MQ=255)                                                        gccCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGTTCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGAt                                                                                             <  7:209289/149‑1 (MQ=255)                                                         cccAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATa                                                                                            <  6:122095/149‑1 (MQ=255)                                                          ccAACCCCAAACCCAGGCCGAACAAACGGCTGCAGGAGCGCTGGCGTTCCATCAGCACCTCGGTGATGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GGCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACTGATAGAAGATCTGATAg                                                                                           <  6:117124/149‑1 (MQ=255)                                                                ccAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCAcc                                                                                      >  8:47181/1‑148 (MQ=255)                                                                ccAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCAcc                                                                                      >  2:264672/1‑148 (MQ=255)                                                                 cAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGACTGCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATag                                                                                                       <  1:59055/131‑1 (MQ=255)                                                                 cAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATag                                                                                                       >  2:59055/1‑130 (MQ=255)                                                                 cAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCAcc                                                                                      >  5:17608/1‑147 (MQ=255)                                                                      ccAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGagc                                                                                  <  8:144072/146‑1 (MQ=255)                                                                      ccAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGagc                                                                                  >  7:144072/1‑146 (MQ=255)                                                                          gCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGa                                                                                                                                                <  7:182910/83‑1 (MQ=255)                                                                          gCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGa                                                                                                                                                >  8:182910/1‑83 (MQ=255)                                                                              aaCAAACGGCTCCAGGAGCGCTGGCTTTACATCATCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAAcc                                                                       <  6:17608/149‑1 (MQ=255)                                                                                 aaaCGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGagcag                                                                                >  3:82720/1‑137 (MQ=255)                                                                                 aaaCGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGagcag                                                                                <  4:82720/137‑1 (MQ=255)                                                                                 aaaCGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCgg                                                                     >  6:36803/1‑148 (MQ=255)                                                                                   aCGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTAc                                                                  >  1:41683/1‑149 (MQ=255)                                                                                   aCGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCATGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTAc                                                                  >  5:265665/1‑149 (MQ=255)                                                                                   aCGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCATGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTAc                                                                  >  5:270159/1‑149 (MQ=255)                                                                                     ggCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCTATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTAccc                                                                <  1:83613/149‑1 (MQ=255)                                                                                     ggCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCTCCGAGCAGCACACAACCGGTAcc                                                                 >  7:63797/1‑148 (MQ=255)                                                                                            ggAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCACCAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGgcgc                                                         <  1:85938/149‑1 (MQ=255)                                                                                                 gcTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTAACCAGGGCGCGCAg                                                     >  8:198541/1‑148 (MQ=255)                                                                                                     gctttCCATCAGCACCAGGGTGTGGATGGTTTTTACCCCAATCACGGGGGCATTGAC‑GCGGGGGTGTCACTTCCTCATTCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGaa                                                   <  1:264672/145‑1 (MQ=255)                                                                                                           cATCAGCACCACGGTGATGAACGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGCTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCACCACACAACCGGTACCCAGCGCGCTCa                                                      >  4:70698/1‑137 (MQ=255)                                                                                                           cATAAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCa                                                      <  3:70698/137‑1 (MQ=255)                                                                                                             tCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGa                                        <  2:41683/149‑1 (MQ=255)                                                                                                                 caccacGGTGAGGATCGGTTGTACCCCCATCACCATGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTAc                                    <  6:270159/149‑1 (MQ=255)                                                                                                                 caccacGGTGAGGATCGGTTGTACCCCCATCACCATGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTAc                                    <  6:265665/149‑1 (MQ=255)                                                                                                                      cGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGc                               <  5:36803/149‑1 (MQ=255)                                                                                                                             gATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAgc‑g                         >  7:22002/1‑147 (MQ=255)                                                                                                                             gATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAgc‑g                         <  8:22002/147‑1 (MQ=255)                                                                                                                                  gTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑TTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTc                                              >  5:79896/1‑122 (MQ=255)                                                                                                                                  gTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCtg                                              <  6:79896/122‑2 (MQ=255)                                                                                                                                  gTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAGC‑GCAGccc                   >  6:78769/1‑148 (MQ=255)                                                                                                                                             aTCACCGTGGCCACGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAGC‑GCAGCCCTAGCAACGGGcc       >  4:215455/1‑149 (MQ=255)                                                                                                                                                aCCGTGGCCACGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAGC‑GCAGCCCTAGCAACGGGCCGAt    <  3:215455/149‑1 (MQ=255)                                                                                                                                                   gTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAGCAGCAGCCCTAGCAACGGGCCGATCa  <  7:47181/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                        CCGCCAGCGATACCCCGCCCAGGTTGATGCCCTGGTAGACCACCATGATCAAGCCGCCCAACCCCAAACCCAGGCCGAACAAACGGCTCCAGGAGCGCTGGCGTTCCATCAGCACCACGGTGAGGATCGGTTGTACCCCCATCACCGTGGCCATGAC‑GCCGGGG‑GTGACATGGGTA‑TCAAGCGCCAACAGATAGAAGATCTGATAGGCACCGAGCAGCACACAACCGGTACCCAGGGCGCGCAGAATGCTCCCGCGAGTACGCGGCCAGC‑GCAGCCCTAGCAACGGGCCGATCA  >  NC_002947/5188440‑5188733

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