Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,376,321	A→T	10.0%	intergenic (+30/+19)	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,321	0	A	T	10.0%	153.2 / 8.1	60	intergenic (+30/+19)	kup/PP_1201	potassium transport system protein Kup/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/31);  new base T (3/3);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

CTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGCC  >  NC_002947/1376180‑1376469                                                                                                                                                 |                                                                                                                                                     cTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGc                                                                                                                                               <  1:20659/149‑1 (MQ=255) cTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGc                                                                                                                                               <  4:210221/149‑1 (MQ=255)    gCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCTTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGTTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                            <  4:158256/149‑1 (MQ=255)       ggggAACATTTTCGCGTTCTTTCTTTATAATTGCAAAGGCAACTTTCGCTTCTTCAACCTTCCCTTTAAACGGGTGATCCTCAAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                            <  1:257800/149‑1 (MQ=255)        ggCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGGCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGAGc                                                                                                                                        <  7:6652/149‑1 (MQ=255)           aaCCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCtt                                                                                                                                     <  5:326751/149‑1 (MQ=255)              cTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCt                                                                                                                                            >  5:140564/1‑139 (MQ=255)              cTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCt                                                                                                                                            <  6:140564/139‑1 (MQ=255)                       gtCTTGCTGACGAATTCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGGGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCt                                                                                                                         <  7:245200/148‑1 (MQ=255)                        tCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGt                                                                                                                                   <  8:263490/138‑1 (MQ=255)                        tCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGt                                                                                                                                   >  7:263490/1‑138 (MQ=255)                             cTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCt                                                                                                                    <  2:71040/148‑1 (MQ=255)                               gaagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCa                                                                                                                                                      >  8:121772/1‑112 (MQ=255)                               gaagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCa                                                                                                                                                      <  7:121772/112‑1 (MQ=255)                                aagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTAc                                                                                                                  <  3:230236/147‑1 (MQ=255)                                aagaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTGATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGc                                                                                                                <  6:136098/149‑1 (MQ=255)                                 agaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCCCAGAGCGGCCCTGAGCATTGGACGTTATGATCAGCTaac                                                                                                                 >  3:155922/1‑145 (MQ=255)                                     tGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGCGCCGTGTCAGTTCTGCTCAGCTACGcgtgcc                                                                                                           >  8:223533/1‑144 (MQ=255)                                       ccAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGCACCGGGTGA‑‑‑TCGAGCTAGGGACCCCGGTCGAGAGCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACtt                                                                                                         <  7:223533/149‑1 (MQ=255)                                                  ttGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGcgcg                                                                                              >  1:171448/1‑149 (MQ=255)                                                    gcgcTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTg                                                                                                   <  7:69281/142‑1 (MQ=255)                                                    gcgcTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTg                                                                                                   >  8:69281/1‑142 (MQ=255)                                                      gcTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTc                                                                                          >  3:206290/1‑149 (MQ=255)                                                       cttcttCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTc                                                                                                                                  <  8:151790/108‑1 (MQ=255)                                                         tcttcAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCtt                                                                                                                                     <  3:1142/103‑1 (MQ=255)                                                            gcaagttGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtctc                                                                                    <  1:289737/143‑1 (MQ=255)                                                             cAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCg                                                                                   >  5:47975/1‑149 (MQ=255)                                                                 cTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCaa                                                                               <  1:145813/149‑1 (MQ=255)                                                                  tGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCGCAGGGCCGCCCAGCGGCCCAGAGCCCTGTGAGTTGTGCTCAGCTACCCGGGGCGTGCCATGGGGCGGTTGGATGTGGCCaa                                                                               >  6:222288/1‑148 (MQ=255)                                                                    ccGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCaa                                                                               <  5:222288/146‑1 (MQ=255)                                                                    ccGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATc                                                                             <  6:47975/148‑1 (MQ=255)                                                                                       tCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGc                                                            >  5:83140/1‑149 (MQ=255)                                                                                         gAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCt                                                                                                               >  1:128867/1‑96 (MQ=255)                                                                                         gAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCt                                                                                                               <  2:128867/96‑1 (MQ=255)                                                                                              aGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGGTCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAg                                                     >  4:66155/1‑149 (MQ=255)                                                                                               gggACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                      <  3:37660/147‑1 (MQ=255)                                                                                               gggACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGt                                                    >  6:110468/1‑149 (MQ=255)                                                                                                ggACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGtt                                                   >  1:230329/1‑149 (MQ=255)                                                                                                 gACCCAGGTCGAGATCTGAATGGGGTTATAAAGCAGGCTCAGGGCCGCAGAGCGGCCCTGATCCTTGTAAGTTCTGCTCAGCTTAGCAGGAATTTCCTTGTCGCGGTTCGATATAGCCCATCAACCGCTTGGCGAGTTCACGGGAGTtc                                                  >  4:30541/1‑149 (MQ=255)                                                                                                   cccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                      <  1:277847/143‑1 (MQ=255)                                                                                                   cccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                      >  2:277847/1‑143 (MQ=255)                                                                                                     cAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgt                                              >  5:320846/1‑149 (MQ=255)                                                                                                      aGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgtc                                             <  6:168455/149‑1 (MQ=255)                                                                                                         tCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGaa                                          >  1:235012/1‑149 (MQ=255)                                                                                                         tCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGaa                                          <  5:110468/149‑1 (MQ=255)                                                                                                                   aaTGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCg                                                          <  1:236533/123‑1 (MQ=255)                                                                                                                        gggtataCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCCGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAgg                           <  2:44065/146‑1 (MQ=255)                                                                                                                         gTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGCTCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGc                          <  3:66155/149‑1 (MQ=255)                                                                                                                         cttATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCAAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGc                          <  8:290501/148‑1 (MQ=255)                                                                                                                             taCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGt                       <  6:320846/148‑1 (MQ=255)                                                                                                                               cagcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGTCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                    <  4:155922/149‑1 (MQ=255)                                                                                                                                  cagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGaa                                          <  7:178103/124‑1 (MQ=255)                                                                                                                                  cagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGaa                                          >  8:178103/1‑124 (MQ=255)                                                                                                                                  cagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAgg                 >  5:42179/1‑149 (MQ=255)                                                                                                                                     gcgcAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgcg               <  7:293822/145‑1 (MQ=255)                                                                                                                                     gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTCGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAATGCCGGGTAGTTCACGTCGAAGTGGTGCCCGCAAGGAAGTTTCAGgcg               >  8:293822/1‑148 (MQ=255)                                                                                                                                      cTCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGc                               >  7:235546/1‑131 (MQ=255)                                                                                                                                      cTCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGc                               <  8:235546/131‑1 (MQ=255)                                                                                                                                       tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCt             >  1:174911/1‑148 (MQ=255)                                                                                                                                        cAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa           >  5:193170/1‑149 (MQ=255)                                                                                                                                            gCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTa       >  5:123969/1‑149 (MQ=255)                                                                                                                                              cGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACt     >  2:105775/1‑149 (MQ=255)                                                                                                                                                 agagCGGACCTGAGCCTTGTCAGTTCTGCTCAGCTACTCTTGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGGAGTTATCGTCGAAGTGGTGCACTCCAGGCAGTTTCAGTCGCTCACATCCTgac  >  2:179546/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGA‑‑‑TCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTGCC  >  NC_002947/1376180‑1376469

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑