Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,376,323	A→T	6.8%	intergenic (+32/+17)	kup → / ← PP_1201	potassium transport system protein Kup/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,376,323	0	A	T	6.8%	167.1 / 2.5	59	intergenic (+32/+17)	kup/PP_1201	potassium transport system protein Kup/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/30);  new base T (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTG  >  NC_002947/1376175‑1376467                                                                                                                                                       |                                                                                                                                                 tctcgcTGGCGGGTCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCaga                                                                                                                                                  <  7:190264/146‑1 (MQ=255) tCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCaga                                                                                                                                                  <  4:102690/149‑1 (MQ=255)         gCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                          >  6:104983/1‑149 (MQ=255)           ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAAGCGGGTGATCGCGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        <  3:22694/149‑1 (MQ=255)           ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                          >  6:262901/1‑147 (MQ=255)           ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  8:216873/149‑1 (MQ=255)           ggggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  3:17371/149‑1 (MQ=255)             ggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGc                                                                                                                                      <  3:130566/149‑1 (MQ=255)                  ccTGTTCGCGTTCTTGCTGCAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGt                                                                                                                                 <  8:151668/149‑1 (MQ=255)                      ttCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  4:86020/138‑1 (MQ=255)                      ttCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        >  3:86020/1‑138 (MQ=255)                      ttCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGt                                                                                                                             >  7:54838/1‑149 (MQ=255)                      ttCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTAAGGGCCGCAGAGCGGCCCAGAGCCTTGTAAGt                                                                                                                             >  6:221281/1‑149 (MQ=255)                         gcgTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTc                                                                                                                           <  5:104983/148‑1 (MQ=255)                                   tGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTAc                                                                                                                >  8:92085/1‑149 (MQ=255)                                   tGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTAc                                                                                                                >  3:264605/1‑149 (MQ=255)                                       gaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTg                                                                                                            <  3:104943/149‑1 (MQ=255)                                         aTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGa                                                                                                          <  6:3552/149‑1 (MQ=255)                                             cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTg                                                                                                      >  2:250561/1‑149 (MQ=255)                                                 ggCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCtt                                                                                                  >  2:298321/1‑149 (MQ=255)                                                          cgcTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGCGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGttt                                                                                         >  2:56720/1‑149 (MQ=255)                                                           gcTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGttt                                                                                         >  1:117621/1‑148 (MQ=255)                                                           gcTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGttt                                                                                         <  2:117621/148‑1 (MQ=255)                                                           gcTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCAGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTc                                                                                        >  6:327075/1‑149 (MQ=255)                                                            cttcttCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCg                                                                                       <  5:262901/149‑1 (MQ=255)                                                             ttcttcAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCCCAGAGCGGCCCAGAGCCATGTGAGTTCTGCTAAGCTAAGCTGGACTTGCCTTTGCGCGtttag                                                                                       >  2:321353/1‑146 (MQ=255)                                                                ttcAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtct                                                                                   >  1:78523/1‑149 (MQ=255)                                                                  cAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtctc                                                                                  <  1:250561/148‑1 (MQ=255)                                                                    aCCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCGGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGc                                                                                >  2:113229/1‑148 (MQ=255)                                                                      cTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtctc                                                                                  <  8:21629/144‑1 (MQ=255)                                                                      cTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCgagat                                                                             >  7:21629/1‑145 (MQ=255)                                                                      cTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCaa                                                                             <  5:221281/149‑1 (MQ=255)                                                                         ccGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATc                                                                           <  7:189258/148‑1 (MQ=255)                                                                           gTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGgc                                                                                               <  5:185948/126‑1 (MQ=255)                                                                            ttGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAGATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAAc                                                                        >  1:62802/1‑149 (MQ=255)                                                                               aaCCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCt                                                                    >  4:95579/1‑149 (MQ=255)                                                                                 ccGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTAAGGGCCGCAGAGCGGCCCAGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCCCTTg                                                                  >  6:109097/1‑149 (MQ=255)                                                                                     gTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGa                                                              <  7:235684/149‑1 (MQ=255)                                                                                     gTGATCGAGCTAGGGACCCAGGGCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGa                                                              <  7:231204/149‑1 (MQ=255)                                                                                        aTCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTg                                                           >  6:147464/1‑149 (MQ=255)                                                                                           gAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGcc                                                         <  2:78523/148‑1 (MQ=255)                                                                                            aGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACTCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                       <  2:62802/149‑1 (MQ=255)                                                                                                 gggACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGATCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                    <  1:250790/147‑1 (MQ=255)                                                                                                     cccAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                              <  8:54838/149‑1 (MQ=255)                                                                                                     cccAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                              <  1:298321/149‑1 (MQ=255)                                                                                                      ccAGGTCGAGATCTGA‑ATGGGGTTATACAGCAAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                              >  1:198776/1‑149 (MQ=255)                                                                                                                     a‑aTGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGccc                               <  1:321353/148‑1 (MQ=255)                                                                                                                     a‑aTGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCg                              >  3:294331/1‑149 (MQ=255)                                                                                                                                   agc‑agGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCa                 <  3:95579/149‑1 (MQ=255)                                                                                                                                    gc‑agGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCTATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg                <  7:92085/149‑1 (MQ=255)                                                                                                                                    gc‑agGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg                <  4:264605/149‑1 (MQ=255)                                                                                                                                         gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCg                              >  1:106143/1‑131 (MQ=255)                                                                                                                                         gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCg                              <  2:106143/131‑1 (MQ=255)                                                                                                                                           tCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt           >  1:196525/1‑148 (MQ=255)                                                                                                                                           tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt           >  3:135341/1‑148 (MQ=255)                                                                                                                                             aGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt           <  2:196525/146‑1 (MQ=255)                                                                                                                                             aGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCAc        >  6:203633/1‑149 (MQ=255)                                                                                                                                             aGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa         <  2:198776/148‑1 (MQ=255)                                                                                                                                             aGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCAc        <  4:294331/149‑1 (MQ=255)                                                                                                                                                  cGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACt   >  3:25622/1‑149 (MQ=255)                                                                                                                                                   gCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTg  <  4:25622/149‑1 (MQ=255)                                                                                                                                                     agagCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCg                                                                                       >  4:66518/1‑62 (MQ=255)                                                                                                                                                     agagCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCg                                                                                       <  3:66518/62‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                 TCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGA‑ATGGGGTTATACAGC‑AGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTG  >  NC_002947/1376175‑1376467

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