Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,594	G→A	7.9%	L455L (CTG→TTG)	ppk ←	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,951,594	0	G	A	7.9%	254.1 / 7.8	88	L455L (CTG→TTG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (42/39);  new base A (4/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

AGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTG  >  NC_002947/5951446‑5951735                                                                                                                                                     |                                                                                                                                               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                                                                                                                                        >  4:89954/1‑136 (MQ=255)                      tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTgg                                                                                                                                        <  3:89954/136‑1 (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                                                                                                       <  4:14532/148‑1 (MQ=255)                                          cgcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                       <  5:328172/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                                                                                            <  3:277503/145‑1 (MQ=255)                                                        tctcGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTc                                                                                            >  4:277503/1‑145 (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                                           >  5:391566/1‑149 (MQ=255)                                                                                                     tCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCg                                           >  2:275037/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:368461/1‑149 (MQ=255)                                                                                                                 ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                               >  7:356806/1‑149 (MQ=255)                                                                                                                  cacaCCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTcc                               >  4:338374/1‑148 (MQ=255)                                                                                                                   acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                             <  3:74519/149‑1 (MQ=255)                                                                                                                   acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGa                             <  2:81891/149‑1 (MQ=255)                                                                                                                     acCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATGCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGCGCGACGACGATTTCagag                            >  1:308489/1‑149 (MQ=255)                                                                                                                         tAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGc                       <  4:304826/149‑1 (MQ=255)                                                                                                                                 tACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGtt                          >  5:177418/1‑138 (MQ=255)                                                                                                                                 tACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGtt                          <  6:177418/138‑1 (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              <  5:176921/148‑1 (MQ=255)                                                                                                                                   ccGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCg              <  3:175894/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    >  3:142249/1‑149 (MQ=255)                                                                                                                                            ccGCCTGCAAGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTc    >  1:322111/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  >  1:281958/1‑149 (MQ=255)                                                                                                                                              gCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTg  >  4:45551/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                               AGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAT‑CTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTG  >  NC_002947/5951446‑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: ‑