Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,589	T→C	9.5%	Q456Q (CAA→CAG)	ppk ←	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,951,589	0	T	C	9.5%	254.2 / 14.2	95	Q456Q (CAA→CAG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (45/41);  new base C (4/5);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

TGTACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTG  >  NC_002947/5951441‑5951735                                                                                                                                                     |                                                                                                                                                    tgtACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCt                                                                                                                                                     >  2:14731/1‑149 (MQ=255)   tACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCtt                                                                                                                                                    <  6:383927/148‑1 (MQ=255)     cAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGc                                                                                                                                                  <  1:409180/148‑1 (MQ=255)      aGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAg                                                                                                                                                <  5:11303/149‑1 (MQ=255)      aGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAg                                                                                                                                                >  4:175894/1‑149 (MQ=255)        gcgggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCGCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggc                                                                                                                                              >  4:148743/1‑149 (MQ=255)          gggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAgg                                                                                                                                               <  2:72534/146‑1 (MQ=255)          gggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAgg                                                                                                                                               >  1:72534/1‑146 (MQ=255)           ggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAagcggc                                                                                                                                           <  5:173102/149‑1 (MQ=255)             cgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGGACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTg                                                                                                                                         <  5:13582/149‑1 (MQ=255)             cgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTg                                                                                                                                         >  3:399235/1‑149 (MQ=255)              gTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTgg                                                                                                                                        >  6:328172/1‑149 (MQ=255)              gTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTgg                                                                                                                                        <  7:8383/149‑1 (MQ=255)               ttGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGc                                                                                                                                       >  8:195908/1‑149 (MQ=255)               ttGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggt                                                                                                                                       <  1:222634/149‑2 (MQ=255)                tGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGc                                                                                                                                       >  2:140383/1‑148 (MQ=255)                  ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCGCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCa                                                                                                                                     <  3:148743/148‑1 (MQ=255)                  ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAt                                                                                                                                    <  3:12393/149‑1 (MQ=255)                          gTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAgg                                                                                                                            <  5:409923/148‑1 (MQ=255)                           tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTgg                                                                                                                                        <  3:89954/136‑1 (MQ=255)                           tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTgg                                                                                                                                        >  4:89954/1‑136 (MQ=255)                           tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGtt                                                                                                                          <  1:426511/149‑1 (MQ=255)                                gCCGGTGCCAACATACGCGTAACGCACGATCTCGCCCTGCTCGAGCCGCAGAATCAGCATCATCTTCTCGAGGGTCTTGAAGCCACACACACCGAAGAACCCTACCGCACCAGCCGCTTTCAGGCGGCTGGCCAT‑CTGCACGTAGGAc                                                                                                                      >  6:326038/1‑148 (MQ=255)                                     tGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCg                                                                                                                <  6:395500/149‑1 (MQ=255)                                     tGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCg                                                                                                                >  3:14532/1‑149 (MQ=255)                                      gCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACCCCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGt                                                                                                               >  3:340359/1‑149 (MQ=255)                                       cccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGt                                                                                                               <  1:14731/148‑1 (MQ=255)                                        ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGTAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCa                                                                                                             <  1:140383/149‑1 (MQ=255)                                        ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTc                                                                                                              >  4:74519/1‑148 (MQ=255)                                        ccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTc                                                                                                              >  5:346580/1‑148 (MQ=255)                                             tggtcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCAGCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAAc                                                                                                         <  8:352395/144‑1 (MQ=255)                                               cgcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                       <  5:328172/148‑1 (MQ=255)                                               cgcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                       <  4:14532/148‑1 (MQ=255)                                                   tAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGcgcacg                                                                                                  <  4:399235/149‑1 (MQ=255)                                                      cGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGGCAT‑CTGCAGGTTGGACTCGTCGTCAAAACGCGCACGCaa                                                                                               >  3:215617/1‑149 (MQ=255)                                                       gCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAAc                                                                                              >  3:79865/1‑149 (MQ=255)                                                        cACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACt                                                                                             <  4:408195/149‑1 (MQ=255)                                                        cACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACt                                                                                             >  1:81891/1‑149 (MQ=255)                                                         aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTc                                                                                            <  8:303499/149‑1 (MQ=255)                                                             tctcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTc                                                                                            >  4:277503/1‑145 (MQ=255)                                                             tctcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTc                                                                                            <  3:277503/145‑1 (MQ=255)                                                               tcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGTCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCAc                                                                                      >  8:277110/1‑149 (MQ=255)                                                               tcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGCTCAc                                                                                      >  3:350861/1‑149 (MQ=255)                                                                  cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGc                                                                                   <  6:346580/149‑1 (MQ=255)                                                                   ccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                  <  7:195908/149‑1 (MQ=255)                                                                       cTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGa                                                                              >  2:315632/1‑149 (MQ=255)                                                                          gcggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGTCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcc                                                                            <  7:277110/148‑1 (MQ=255)                                                                          gcggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAc                                                                             <  4:340359/147‑1 (MQ=255)                                                                               gcAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTg                                                                      <  1:64068/149‑1 (MQ=255)                                                                                         catcatCTTGGGGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAACCGCGCACGCAACTCGATCCCCGCAGTGACCTCCTGGCCGGTAcgcg                                                            >  4:269881/1‑149 (MQ=255)                                                                                           tcatcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGcgcc                                                          >  6:40238/1‑149 (MQ=255)                                                                                               cTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcct                                                                           <  6:74977/128‑1 (MQ=255)                                                                                                ttGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCCGTGACCTCCTTGCCGTTACGCGCCGCATc                                                     >  5:279936/1‑149 (MQ=255)                                                                                                  ggCGTGGGTCTTGAAGCCAACGACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATcc                                                    <  7:317812/148‑1 (MQ=255)                                                                                                  ggCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATc                                                     <  7:244940/147‑1 (MQ=255)                                                                                                     gTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATccac                                                  >  1:323946/1‑147 (MQ=255)                                                                                                          tCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCg                                           >  2:275037/1‑149 (MQ=255)                                                                                                          tCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCg                                           >  5:391566/1‑149 (MQ=255)                                                                                                             tGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTc                                         >  2:319180/1‑148 (MQ=255)                                                                                                              gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTAGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACAAGGGCGTcga                                       >  7:382784/1‑149 (MQ=255)                                                                                                               aaGCCAGCCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGAGTCTTCGTCAACACGCGCCAGCCACTCGATCACCGCAGTGACCTCCTTGGCGTTACGCGCCGCAGCCCGCTGGGGGTcggc                                      >  8:312956/1‑147 (MQ=255)                                                                                                               aaGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTcgac                                      >  6:15559/1‑149 (MQ=255)                                                                                                                 gCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGACAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTcgacga                                    >  3:304826/1‑149 (MQ=255)                                                                                                                  ccAACCACACCGTAGATTACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGAt                                   >  2:68698/1‑149 (MQ=255)                                                                                                                      ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                               >  7:356806/1‑149 (MQ=255)                                                                                                                      ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                               >  7:368461/1‑149 (MQ=255)                                                                                                                       cacaCCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                               >  4:338374/1‑148 (MQ=255)                                                                                                                        acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                             <  2:81891/149‑1 (MQ=255)                                                                                                                        acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                             <  3:74519/149‑1 (MQ=255)                                                                                                                          acCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATGCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGCGCGACGACGATTTCagag                            >  1:308489/1‑149 (MQ=255)                                                                                                                              tAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGc                       <  4:304826/149‑1 (MQ=255)                                                                                                                                      tACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGtt                          <  6:177418/138‑1 (MQ=255)                                                                                                                                      tACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGtt                          >  5:177418/1‑138 (MQ=255)                                                                                                                                      tACCGCACCACCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCATTGACCTCCTTGCCGTTACGAGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc               <  1:315632/149‑1 (MQ=255)                                                                                                                                       aCCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCg              >  8:294589/1‑149 (MQ=255)                                                                                                                                        gatggccAGCCGCCTGCAAGCGGCTGGTCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCgg             >  3:164150/6‑149 (MQ=255)                                                                                                                                        ccGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTa                         <  8:304466/137‑1 (MQ=255)                                                                                                                                        ccGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTa                         >  7:304466/1‑137 (MQ=255)                                                                                                                                        ccGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc               <  5:40238/147‑1 (MQ=255)                                                                                                                                        ccGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc               <  5:340085/147‑1 (MQ=255)                                                                                                                                        ccGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCg              <  3:175894/148‑1 (MQ=255)                                                                                                                                        ccGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCg              <  5:176921/148‑1 (MQ=255)                                                                                                                                         cGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTc                                <  2:200735/129‑1 (MQ=255)                                                                                                                                          tggccAGCCGCCTGCAAGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTa           >  2:411966/4‑149 (MQ=255)                                                                                                                                           cACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTa           >  1:25420/1‑148 (MQ=255)                                                                                                                                            gccAGCCGCCTGCAAGCGGCTGGTCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACa         <  4:164150/148‑1 (MQ=255)                                                                                                                                              cAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg       >  2:249495/1‑149 (MQ=255)                                                                                                                                              cAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg       <  1:275037/149‑1 (MQ=255)                                                                                                                                                gCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGt     >  4:159503/1‑149 (MQ=255)                                                                                                                                                 ccGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTc    <  3:269881/149‑1 (MQ=255)                                                                                                                                                 ccGCCTGCAAGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTc    >  1:322111/1‑149 (MQ=255)                                                                                                                                                 ccGCCTGCAAGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTc    >  3:142249/1‑149 (MQ=255)                                                                                                                                                  cGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCt   <  8:368461/149‑1 (MQ=255)                                                                                                                                                   gCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGt            <  7:43743/139‑1 (MQ=255)                                                                                                                                                   gCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGt            >  8:43743/1‑139 (MQ=255)                                                                                                                                                   gCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTg  >  4:45551/1‑149 (MQ=255)                                                                                                                                                   gCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTg  >  1:281958/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGTACAGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTG  >  NC_002947/5951441‑5951735

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