Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,951,594	G→A	8.3%	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	8.3%	242.7 / 8.6	84	L455L (CTG→TTG)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (37/40);  new base A (5/2);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCGCAC  >  NC_002947/5951446‑5951741                                                                                                                                                     |                                                                                                                                                     aGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAg                                                                                                                                                      <  4:352566/149‑1 (MQ=255)   gcgggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggc                                                                                                                                                    >  4:39781/1‑149 (MQ=255)     gggcgTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                                  <  6:312064/149‑1 (MQ=255)            gCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGcc                                                                                                                                            <  8:121922/148‑1 (MQ=255)            gCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggtgc                                                                                                                                           >  4:347003/1‑146 (MQ=255)             ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAt                                                                                                                                          <  8:248382/149‑1 (MQ=255)             ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCAt                                                                                                                                          >  8:235290/1‑149 (MQ=255)             ccGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCCTGCAAGCGGCTggtgcg                                                                                                                                          <  3:347003/149‑5 (MQ=255)               ggCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATca                                                                                                                                        <  1:373194/149‑2 (MQ=255)               ggCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCt                                                                                                                                        >  6:244913/1‑149 (MQ=255)                gCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                                  <  7:116195/138‑1 (MQ=255)                gCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggcgg                                                                                                                                                  >  8:116195/1‑138 (MQ=255)                 cGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCATAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGc                                                                                                                                      <  2:151642/149‑1 (MQ=255)                  gTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTg                                                                                                                                       >  4:115660/1‑147 (MQ=255)                     gTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCACCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGt                                                                                                                                  <  3:371012/149‑1 (MQ=255)                      tAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGtt                                                                                                                                 <  2:101212/149‑1 (MQ=255)                         ttGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGa                                                                                                                              <  3:149884/149‑1 (MQ=255)                               gTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                        <  6:117153/149‑1 (MQ=255)                               gTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTc                                                                                                                        <  7:118797/149‑1 (MQ=255)                                  cccAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCATCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGt                                                                                                                      <  1:269272/148‑1 (MQ=255)                                    cAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCggaa                                                                                                                    >  1:141846/1‑145 (MQ=255)                                      gATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCaaaa                                                                                                                 <  8:123511/149‑1 (MQ=255)                                       aTGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATAATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAc                                                                                                                >  5:178635/1‑149 (MQ=255)                                          cgcgTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAAcgc                                                                                                              <  3:212291/148‑1 (MQ=255)                                               aaCGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcacg                                                                                                         >  2:366736/1‑148 (MQ=255)                                                    aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggc                                                                                                                                                    <  5:229918/100‑1 (MQ=255)                                                    aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAggc                                                                                                                                                    >  6:229918/1‑100 (MQ=255)                                                    aCGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTc                                                                                                   <  4:323303/149‑1 (MQ=255)                                                            gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGAt                                                                                                <  7:77678/144‑1 (MQ=255)                                                            gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                           <  1:366736/149‑1 (MQ=255)                                                            gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                           >  8:231006/1‑149 (MQ=255)                                                            gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                           >  2:215195/1‑149 (MQ=255)                                                            gCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                           <  3:39781/149‑1 (MQ=255)                                                             cccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCg                                                                                           >  8:66277/1‑148 (MQ=255)                                                              ccTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCa                                                                                         <  5:244913/149‑1 (MQ=255)                                                                 gCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                       >  5:326520/1‑148 (MQ=255)                                                                  cTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGCAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                       <  1:204305/147‑1 (MQ=255)                                                                  cTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                       <  6:112931/147‑1 (MQ=255)                                                                  cTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGt                                                                                       <  1:132296/147‑1 (MQ=255)                                                                   tCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAc                                                                                    >  1:318346/1‑149 (MQ=255)                                                                     gcggcgCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcc                                                                                   <  3:212717/148‑1 (MQ=255)                                                                        gcgcAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctcct                                                                               >  3:232128/1‑149 (MQ=255)                                                                           cAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGc                                                                            >  8:336414/1‑149 (MQ=255)                                                                            aGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCCCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCCCCGCAGTGACCTCCTTGcc                                                                           >  4:29679/1‑149 (MQ=255)                                                                               aTCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGtt                                                                        <  4:14128/149‑1 (MQ=255)                                                                               aTCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGtt                                                                        <  5:123771/149‑1 (MQ=255)                                                                               aTCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGt                                                                         >  2:338157/1‑148 (MQ=255)                                                                                tCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTa                                                                       >  2:326305/1‑149 (MQ=255)                                                                                     atcatcTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCATTGACCTCCTTGCCGTTACgcgc                                                                  <  7:235290/149‑1 (MQ=255)                                                                                              gCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑ccac                                                        >  8:127722/1‑149 (MQ=255)                                                                                              gCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCGTGCCGTTACGCGCCTCAT‑ccac                                                        >  7:343772/1‑149 (MQ=255)                                                                                                gTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCA‑CAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAg                                                     >  1:251640/1‑149 (MQ=255)                                                                                                 tGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAg                                                     >  5:240390/1‑149 (MQ=255)                                                                                                 tGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAg                                                     >  8:264768/1‑149 (MQ=255)                                                                                                         gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctc                                                                                 >  2:260244/1‑114 (MQ=255)                                                                                                         gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGAcctc                                                                                 <  1:260244/114‑1 (MQ=255)                                                                                                         gAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTc                                               <  7:66277/147‑1 (MQ=255)                                                                                                             ccAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGAt                                         <  3:111289/149‑1 (MQ=255)                                                                                                              cAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGt                                                <  1:36959/141‑1 (MQ=255)                                                                                                              cAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGt                                                >  2:36959/1‑141 (MQ=255)                                                                                                                 ccACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTcc                                     <  6:326520/149‑1 (MQ=255)                                                                                                                   acacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGa                                   <  1:326305/149‑1 (MQ=255)                                                                                                                    cacCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAg                                  <  2:318346/149‑1 (MQ=255)                                                                                                                      ccGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCATCCCACCAGGGCGTCGACGATTTCCGAGt                                 <  1:215195/149‑1 (MQ=255)                                                                                                                          aGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcac                                                                                                          >  8:247732/1‑72 (MQ=255)                                                                                                                          aGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGcgcac                                                                                                          <  7:247732/72‑1 (MQ=255)                                                                                                                             tCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGccccc                         >  6:367669/1‑149 (MQ=255)                                                                                                                             tCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGccccc                         >  5:61138/1‑149 (MQ=255)                                                                                                                               aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGa                       <  7:127722/149‑1 (MQ=255)                                                                                                                               aCTACCGCACCAGCCGCTTGCAGGCGGCTGGCAATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGccccc                         >  6:129228/1‑147 (MQ=255)                                                                                                                                 tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGtt                                <  7:18181/138‑1 (MQ=255)                                                                                                                                 tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGtt                                >  8:18181/1‑138 (MQ=255)                                                                                                                                 tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc                     >  7:60767/1‑149 (MQ=255)                                                                                                                                 tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGAGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGc                     >  1:66197/1‑149 (MQ=255)                                                                                                                                 tACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTACGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGACGACGAATTCCGAGTTAGCCACCGAGc                     >  3:72710/1‑149 (MQ=255)                                                                                                                                  ggatggccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCg                    <  6:251626/143‑1 (MQ=255)                                                                                                                                   ccGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCAGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCagcggagcgg                   >  8:140428/1‑149 (MQ=255)                                                                                                                                     tggccAGCCGCCTGCAAGCGGCTGGTCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTa                 >  7:101465/4‑149 (MQ=255)                                                                                                                                      cACCAGCCGCTTGCAGGCGGCTGGCCAGCTGCAGGTTGGAGTCTGCGTCCAAACGCGCCCGCAACTCGATGACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAG‑GCGCCAGGGCGTCGCCGAGTTCCGAGGTAGCCGCCGAGCGGTag                >  4:222394/1‑148 (MQ=255)                                                                                                                                       aCCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTa                 >  4:245433/1‑147 (MQ=255)                                                                                                                                       aCCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACa               <  3:115660/149‑1 (MQ=255)                                                                                                                                        ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAg              >  7:309005/1‑149 (MQ=255)                                                                                                                                        ccAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAg              >  3:168474/1‑149 (MQ=255)                                                                                                                                         cAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg             <  5:367669/149‑1 (MQ=255)                                                                                                                                         cAGCCGCCTGCAAGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg             >  1:225462/1‑149 (MQ=255)                                                                                                                                          aGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAggg            >  3:127221/1‑149 (MQ=255)                                                                                                                                           gCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg             >  7:29176/1‑147 (MQ=255)                                                                                                                                           gCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAgg             <  8:29176/147‑1 (MQ=255)                                                                                                                                              gCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTg        <  1:338157/149‑1 (MQ=255)                                                                                                                                               cTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTAC‑GGGTCTGgc      <  2:251640/149‑1 (MQ=255)                                                                                                                                                    aGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCGCAc  >  6:245664/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGGCGGGCGTTGCCGGCGTGGTAGTTGCCGGTGCCCAGATGCGCGTAACGCACGATCTCGCCCTGCTCGCGGCGCAGAATCAGCATCATCTTGGCGTGGGTCTTGAAGCCAACCACACCGTAGATCACTACCGCACCAGCCGCTTGCAGGCGGCTGGCCATCTGCAGGTTGGACTCTTCGTCAAAACGCGCACGCAACTCGATCACCGCAGTGACCTCCTTGCCGTTACGCGCCGCAT‑CCACCAGGGCGTCGACGATTTCCGAGTTAGCCCCCGAGCGGTACAGGGTCTGGCGCAC  >  NC_002947/5951446‑5951741

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