Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	730,744	A→T	5.5%	D95V (GAC→GTC)	clpB →	chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	730,744	0	A	T	5.5%	262.6 / 2.5	83	D95V (GAC→GTC)	clpB	chaperone protein ClpB
Reads supporting (aligned to +/- strand):  ref base A (37/41);  new base T (2/3);  total (39/44)
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

CCGCTGCTGATGCAGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCATCAACAACCTGCGT  >  NC_002947/730602‑730892                                                                                                                                               |                                                                                                                                                      ccGCTGCTGATGCAGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCt                                                                                                                                                 <  4:107474/149‑1 (MQ=255)   gctgctGATGCAGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCGGTCTGCCTgg                                                                                                                                               <  1:241907/149‑1 (MQ=255)          aTGCAGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTgg                                                                                                                                               >  6:385256/1‑142 (MQ=255)          aTGCAGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTgg                                                                                                                                               <  5:385256/142‑1 (MQ=255)              aGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAg                                                                                                                                    >  1:300928/1‑149 (MQ=255)              aGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACCTGAACAAGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGTCACAGCAGAAg                                                                                                                                    >  5:398983/1‑149 (MQ=255)               ggTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAgg                                                                                                                                   <  8:198231/149‑1 (MQ=255)                    gCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGAc                                                                                                                              <  2:92949/149‑1 (MQ=255)                      ttCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGACACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGAc                                                                                                                              <  6:37213/147‑1 (MQ=255)                      ttCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGACACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGAc                                                                                                                              >  5:37213/1‑147 (MQ=255)                        cGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCa                                                                                                                            <  2:201578/147‑1 (MQ=255)                            agcactGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCAt                                                                                                                      <  6:108696/144‑1 (MQ=255)                               aaTGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGc                                                                                                                                 >  3:335083/1‑135 (MQ=255)                               aaTGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGc                                                                                                                                 <  4:335083/135‑1 (MQ=255)                               aaTGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTc                                                                                                                   >  4:319787/1‑149 (MQ=255)                               aaTGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTGTGCCAGGCGGTCTGCCTGGTTAAGCAGacgt                                                                                                                                   >  8:168404/1‑130 (MQ=37)                               aaTGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTGTGCCAGGCGGTCTGCCTGGTTAAGCAGacgt                                                                                                                                   <  7:168404/133‑4 (MQ=37)                                 tGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCa                                                                                                                 <  3:398374/149‑1 (MQ=255)                                     cTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGa                                                                                                             <  4:12635/149‑1 (MQ=255)                                         gcCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTg                                                                                                         <  6:266615/149‑1 (MQ=255)                                           cAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGt                                                                                                       >  2:45623/1‑149 (MQ=255)                                           cAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGt                                                                                                       >  3:266368/1‑149 (MQ=255)                                             ggggCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCa                                                                                                                 >  1:15945/1‑137 (MQ=255)                                             ggggCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCa                                                                                                                 <  2:15945/137‑1 (MQ=255)                                             ggggCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGc                                                                                                     <  3:319787/149‑1 (MQ=255)                                                   ggTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTgg                                                                                               >  5:323732/1‑149 (MQ=255)                                                    gTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCt                                                                                                    <  7:284347/143‑1 (MQ=255)                                                    gTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCt                                                                                                    >  8:284347/1‑143 (MQ=255)                                                    gTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTgg                                                                                               <  5:362611/148‑1 (MQ=255)                                                      gAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGgccg                                                                                            >  1:401015/1‑149 (MQ=255)                                                      gAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGgccg                                                                                            >  5:265153/1‑149 (MQ=255)                                                         agagCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCTACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGgc                                                                                              <  6:25634/144‑1 (MQ=255)                                                         agagCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGgc                                                                                              >  5:25634/1‑144 (MQ=255)                                                              tCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGAACAGCAGAAGGGCGACCAGTTAATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGAc                                                                                    >  3:222274/1‑149 (MQ=255)                                                                gACCAACTGCCGAAAATACAGAACCCTACCGGAGACGTGAACATGTCGCAGGACTTGGAACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACg                                                                                   <  4:105018/148‑1 (MQ=255)                                                                  ccAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGaa                                                                                 >  7:254890/1‑148 (MQ=255)                                                                       tGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAg                                                                            >  7:343388/1‑148 (MQ=255)                                                                         ccGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGc                                                                           >  7:231228/1‑147 (MQ=255)                                                                         ccGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGc                                                                           <  8:254890/147‑1 (MQ=255)                                                                         ccGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCaa                                                                         <  4:241481/149‑1 (MQ=255)                                                                           gAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑Gc                                                                      >  6:400338/1‑149 (MQ=255)                                                                             aaaTACAGAACCCTACCGGCGACGTGAACCTGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTc                                                                    <  3:115415/149‑1 (MQ=255)                                                                             aaaTACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTc                                                                    <  6:323732/149‑1 (MQ=255)                                                                              aaTACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCg                                                                   <  8:329682/149‑1 (MQ=255)                                                                                   aGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTcc                                                                                                                  >  2:347439/1‑98 (MQ=255)                                                                                   aGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTcc                                                                                                                  <  1:347439/98‑1 (MQ=255)                                                                                   aGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAg                                                              >  6:153664/1‑149 (MQ=255)                                                                                        ccTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAgctgct                                                         >  6:229853/1‑149 (MQ=255)                                                                                          tACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCg                                                                   <  3:116398/137‑1 (MQ=255)                                                                                          tACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCg                                                                   >  4:116398/1‑137 (MQ=255)                                                                                          tACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAgctgctgc                                                       >  7:252284/1‑149 (MQ=255)                                                                                             cGGCGACGTGAACATGTCGCAGGACTTTGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGa                                                    <  6:265153/149‑1 (MQ=255)                                                                                             cGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGAc                                                                                    >  8:304880/1‑118 (MQ=255)                                                                                             cGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGAc                                                                                    <  7:304880/118‑1 (MQ=255)                                                                                                     tGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGc                                            <  2:401015/149‑1 (MQ=255)                                                                                                        aCATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGt                                          <  8:343388/148‑1 (MQ=255)                                                                                                         cATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGa                                        <  8:90817/149‑1 (MQ=255)                                                                                                            gTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTgcgc                                      >  2:126579/1‑148 (MQ=255)                                                                                                            gTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGc                                      >  7:215472/1‑148 (MQ=255)                                                                                                             tCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAAc                                                                              >  7:110791/1‑108 (MQ=255)                                                                                                             tCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAAc                                                                              <  8:110791/108‑1 (MQ=255)                                                                                                             tCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑g                                                                       <  8:209340/114‑1 (MQ=255)                                                                                                             tCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑g                                                                       >  7:209340/1‑114 (MQ=255)                                                                                                               gCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAg                                               >  1:177619/1‑136 (MQ=255)                                                                                                               gCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAg                                               <  2:177619/136‑1 (MQ=255)                                                                                                               gCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGCGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGc                                      >  6:336866/1‑145 (MQ=255)                                                                                                               gCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGCGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGc                                      <  5:336866/145‑1 (MQ=255)                                                                                                                cAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGGTGAGCCCGGGCGTGAGCacg                                   >  8:314454/1‑145 (MQ=255)                                                                                                                cAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAgg                                              >  7:163636/1‑136 (MQ=255)                                                                                                                cAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAgg                                              <  8:163636/136‑1 (MQ=255)                                                                                                                        ggCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAACACTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCt                           <  1:290948/148‑1 (MQ=255)                                                                                                                        ggCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCt                           <  7:314454/147‑1 (MQ=255)                                                                                                                           aCGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGaa                       >  6:268285/1‑148 (MQ=255)                                                                                                                             gTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGt                                                                                                       <  4:69005/67‑1 (MQ=255)                                                                                                                             gTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGt                                                                                                       >  3:69005/1‑67 (MQ=255)                                                                                                                             gTCTGCTTAACCAGGCAGACCGCCTGACACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCCCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAAt                    >  1:219280/1‑149 (MQ=255)                                                                                                                                tgctgtgccAGGCGGTCTGCCTGGTTAAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGc                  >  2:348544/8‑148 (MQ=255)                                                                                                                                 gCTTAACCAGGCAGTCCGCCTGGCACAGCTGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGGTGTTGGCCGCAATGGACGAAAACAGCAA‑GCTAGGCAAGCTGCTGCTGAGCCAGGGCGTGTGCAAGAACGCGCTGGAAAATGcc                 <  4:222274/148‑1 (MQ=255)                                                                                                                                  cTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGc                  <  8:215472/146‑1 (MQ=255)                                                                                                                                  cTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCAt               <  5:229853/149‑1 (MQ=255)                                                                                                                                  cTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCAt               <  4:266368/149‑1 (MQ=255)                                                                                                                                     aaCCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCATc              <  5:268285/147‑1 (MQ=255)                                                                                                                                             aGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCATCAACAACCTGc    >  8:143501/1‑149 (MQ=255)                                                                                                                                               aCCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCATCAACAACCTGCGt  >  1:184885/1‑149 (MQ=255)                                                                                                                                               aCCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCATCAACAACCTGCGt  <  4:403381/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                      CCGCTGCTGATGCAGGTAGGCTTCGACATCAATGGCCTGCGCCAGGGGCTGGTGAAAGAGCTCGACCAACTGCCGAAAATACAGAACCCTACCGGCGACGTGAACATGTCGCAGGACTTGGCACGTCTGCTTAACCAGGCAGACCGCCTGGCACAGCAGAAGGGCGACCAGTTCATTTCCAGCGAGCTGGTGCTGTTGGCCGCCATGGACGAAAACAGCAA‑GCTCGGCAAGCTGCTGCTGAGCCAGGGCGTGAGCAAGAAAGCGCTGGAAAATGCCATCAACAACCTGCGT  >  NC_002947/730602‑730892

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