Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,486,711	T→A	9.1%	L14M (TTG→ATG) ‡	PP_4824 →	two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,486,711	0	T	A	9.1%	186.1 / 8.1	66	L14M (TTG→ATG) ‡	PP_4824	two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base A (3/3);  total (31/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.82e-01

CCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCACTGCCTCT  >  NC_002947/5486566‑5486850                                                                                                                                                  |                                                                                                                                            ccTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGc                                                                                                                                          <  1:309561/149‑1 (MQ=255)          tCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCAAGCTGAcgatca                                                                                                                                 <  1:242431/149‑6 (MQ=255)            cGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGct                                                                                                                                <  4:34180/148‑1 (MQ=255)                tcgtcgACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCGAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCa                                                                                                                           <  8:71824/149‑1 (MQ=255)                  gtggACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGctct                                                                                                                              <  6:388420/141‑1 (MQ=255)                  gtcgACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCt                                                                                                                                         >  6:246676/1‑133 (MQ=255)                  gtcgACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCt                                                                                                                                         <  5:246676/133‑1 (MQ=255)                  gtcgACAGGGCCCGGAATGCGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGctct                                                                                                                              >  5:388420/1‑144 (MQ=255)                     gACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGttt                                                                                                                      <  5:238975/149‑1 (MQ=255)                       cAGGGCCCGGATTGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTcc                                                                                                                    <  8:361716/149‑1 (MQ=255)                          ggCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCAAGCTGAcgat                                                                                                                                   >  4:56397/1‑128 (MQ=255)                          ggCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCAAGCTGAcgat                                                                                                                                   <  3:56397/131‑4 (MQ=255)                           gCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGa                                                                                                                                       >  5:167160/1‑126 (MQ=255)                           gCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGa                                                                                                                                       <  6:167160/126‑1 (MQ=255)                             ccGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGc                                                                                                              <  8:206441/149‑1 (MQ=255)                                        cttgccttCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGcc                                                                                                          >  4:222876/1‑142 (MQ=255)                                        cttgccttCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGcc                                                                                                          <  3:222876/142‑1 (MQ=255)                                          tgccttCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGAc                                                                                                 >  2:107318/1‑149 (MQ=255)                                              ttCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTg                                                                                             >  5:263209/1‑149 (MQ=255)                                                cGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTgg                                                                                            <  6:386776/148‑1 (MQ=255)                                                  aCTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCtt                                                                                                                                            >  1:360654/1‑98 (MQ=255)                                                  aCTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCtt                                                                                                                                            <  2:360654/98‑1 (MQ=255)                                                      ggCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGAc                                                                                      <  2:162566/148‑1 (MQ=255)                                                                 cATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTcga                                                                          <  1:107318/149‑1 (MQ=255)                                                                 cATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCAAGCTGAcgat                                                                                                                                   <  6:383766/92‑4 (MQ=255)                                                                 cATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCAAGCTGAcgat                                                                                                                                   >  5:383766/1‑89 (MQ=255)                                                                    aGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTc                                                                                 >  7:143679/1‑139 (MQ=255)                                                                    aGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTc                                                                                 <  8:143679/139‑1 (MQ=255)                                                                        aGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAg                                                                   <  4:162830/149‑1 (MQ=255)                                                                        aGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAg                                                                   >  3:162830/1‑149 (MQ=255)                                                                          tACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGc                                                                 <  6:272405/149‑1 (MQ=255)                                                                            cccGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGcc                                                                >  7:161504/1‑148 (MQ=255)                                                                              cGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAAc                                                             <  7:263718/149‑1 (MQ=255)                                                                                    aTCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCaa                                                       <  3:59934/149‑1 (MQ=255)                                                                                          cAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGc                                                 >  3:358866/1‑149 (MQ=255)                                                                                           aagaagGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGc                                                 >  5:384725/1‑148 (MQ=255)                                                                                             gaagGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAGCAGGCCaa                                                              >  8:280378/1‑133 (MQ=255)                                                                                             gaagGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCaa                                                              <  7:280378/133‑1 (MQ=255)                                                                                             gaagGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAg                                                  <  7:75400/145‑1 (MQ=255)                                                                                             gaagGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAg                                                  >  8:75400/1‑145 (MQ=255)                                                                                              aagGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGAtg                                             >  1:37631/1‑149 (MQ=255)                                                                                                gggATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTggttgcga                                             >  5:190910/1‑140 (MQ=255)                                                                                                gggATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGAtgtg                                           <  6:263209/149‑1 (MQ=255)                                                                                                  gATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTgcg                                         >  6:71759/1‑149 (MQ=255)                                                                                                       cGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTcc                                     <  3:316646/148‑1 (MQ=255)                                                                                                           gTGCGTCCGCTTCGGATTGCCTCACCTCGGATCGTCACCTTGCCGACCCTGCTCTGCATGTTTCCGGCTGCGGCCTAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAgc                                 <  4:212080/148‑1 (MQ=255)                                                                                                            tGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAgcgc                               >  3:274849/1‑149 (MQ=255)                                                                                                              cGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCt                              <  4:358866/148‑1 (MQ=255)                                                                                                                tCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCt                                       <  5:207465/137‑1 (MQ=255)                                                                                                                tCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGAGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCt                                       >  6:207465/1‑137 (MQ=255)                                                                                                                    cTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCaa                       >  4:206602/1‑149 (MQ=255)                                                                                                                     ttCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCac                       >  1:323539/1‑147 (MQ=255)                                                                                                                       cGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAAc                      <  2:323539/147‑1 (MQ=255)                                                                                                                         gATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAAc                                                      >  3:209394/1‑113 (MQ=255)                                                                                                                         gATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAAc                                                      <  4:209394/113‑1 (MQ=255)                                                                                                                              ccTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAg              <  8:161504/148‑1 (MQ=255)                                                                                                                               cTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGcc            >  2:145671/1‑149 (MQ=255)                                                                                                                               cTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGcc            >  1:250830/1‑149 (MQ=255)                                                                                                                               cTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGcc            >  8:42740/1‑149 (MQ=255)                                                                                                                                  aGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCAc         <  6:130685/149‑1 (MQ=255)                                                                                                                                   gCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCACt        <  3:206602/149‑1 (MQ=255)                                                                                                                                    ctctGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCACTg       >  3:108685/1‑149 (MQ=255)                                                                                                                                     tctGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGCCGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCACTGc      >  8:59227/1‑149 (MQ=255)                                                                                                                                     tctGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCACTGc      <  4:274849/149‑1 (MQ=255)                                                                                                                                        gATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCACTGCctc   >  5:242294/1‑149 (MQ=255)                                                                                                                                         agggTCAGCAAGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGCGCGATGTGCGCTCCGAGCGCTATCGCCACCAGTTCAGCCCACTGCcgct  >  8:45182/4‑149 (MQ=255)                                                                                                                                                 cTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGc                                                                              <  5:24541/65‑1 (MQ=255)                                                                                                                                                 cTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGc                                                                              >  6:24541/1‑65 (MQ=255)                                                                                                                                                  |                                                                                                                                            CCTGCGGCATCCGGCTCGTCGACAGGGCCCGGAATGGGCCTTGCCTTCGACTTGGCGCGCCGCGCATAGTGAGTACCCGGCACATCGGCCAAGAAGGGATTTCGTAGTGCGTCGGCTTCGGATTGCCTCAGCTCTGATCGTCAGCTTGCTGACCCTGCTCTGCATGTTTCCGGCTGCGGCCGAACATGACGGTGGCTGGACCGTTCTGCTCGACGAACAGGCCAACCTGCAACTGAGCGATGTGCGCTCCGAGCGCTATCGCAACCAGTTCAGCCCACTGCCTCT  >  NC_002947/5486566‑5486850

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