Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,376,329	A→T	10.7%	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,329	0	A	T	10.7%	218.9 / 15.2	77	intergenic (+32/+17)	kup/PP_1201	potassium transport system protein Kup/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base T (4/4);  total (39/38)
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.92e-01

TCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTG  >  NC_002947/1376181‑1376473                                                                                                                                                     |                                                                                                                                                 tCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCagag                                                                                                                                                 >  2:1034652/1‑150 (MQ=255) tCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCagag                                                                                                                                                 <  2:1562524/150‑1 (MQ=255)  cgcgCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGc                                                                                                                                                <  2:1169167/150‑1 (MQ=255)   gcgcTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCg                                                                                                                                               <  2:1640915/150‑1 (MQ=255)     gcTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGc                                                                                                                                             <  1:268033/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        <  1:708210/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        <  1:2101404/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        <  1:2101222/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  1:1034652/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  1:467978/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  1:465055/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  2:1650698/150‑1 (MQ=255)          cGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  2:1834681/150‑1 (MQ=255)              gCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                          <  2:1693690/144‑1 (MQ=255)              gCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                          >  1:1693690/1‑144 (MQ=255)                   cTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCa                                                                                                                               >  2:134728/1‑150 (MQ=255)                            ttCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        >  1:1896920/1‑132 (MQ=255)                            ttCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  2:1896920/132‑1 (MQ=255)                                       gaaTGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTgg                                                                                                           >  2:2377818/1‑150 (MQ=255)                                        aaTGCTAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGa                                                                                                          <  2:2307114/150‑1 (MQ=255)                                             cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        >  1:1924603/1‑115 (MQ=255)                                             cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  2:1924603/115‑1 (MQ=255)                                             cAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGAGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGCGCGGGCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGc                                                                                                     >  1:2427710/1‑150 (MQ=255)                                               aCGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTGTGGTGTGCTCAGCTACGCTGGACTTGCCt                                                                                                   >  1:1818918/1‑150 (MQ=255)                                                cGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCtt                                                                                                  <  1:134728/150‑1 (MQ=255)                                                             ttcttcAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAt                                                                                     >  2:1448945/1‑150 (MQ=255)                                                                       tGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        >  2:420867/1‑89 (MQ=255)                                                                       tGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGa                                                                                                                                        <  1:420867/89‑1 (MQ=255)                                                                          cGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTGGCCTTGGCGCGTTTCGATCTCGCCAATCAAc                                                                        <  2:2105022/150‑1 (MQ=255)                                                                           gTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTg                                                                                                            <  2:832457/113‑1 (MQ=255)                                                                           gTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTg                                                                                                            >  1:832457/1‑113 (MQ=255)                                                                              gAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACtt                                                                                                       <  2:915449/115‑1 (MQ=255)                                                                              gAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACtt                                                                                                       >  1:915449/1‑115 (MQ=255)                                                                                         tCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGcc                                                         >  1:379866/1‑150 (MQ=255)                                                                                           gAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                       >  2:321086/1‑150 (MQ=255)                                                                                            aGCTAGGGACCCAGGTCGAGATCTGACTGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCggg                                                      >  2:712051/1‑150 (MQ=255)                                                                                              cTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        <  2:57959/66‑1 (MQ=25)                                                                                              cTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                        >  1:57959/1‑66 (MQ=25)                                                                                                      ccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGTGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgt                                            >  1:88224/1‑150 (MQ=255)                                                                                                      ccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTgg                                                                                                <  1:2169070/98‑1 (MQ=255)                                                                                                      ccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTgg                                                                                                >  2:2169070/1‑98 (MQ=255)                                                                                                             gagaTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTGCTCGTCGAAgtg                                     <  2:2436198/150‑1 (MQ=255)                                                                                                              agaTCTGAATGGGGTTATACAGCAGCGTTCGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtgg                                    >  1:1243635/1‑150 (MQ=255)                                                                                                                agcTGAATGGGGTTATACAGCAGGCTCAGGGCCGACGAGTTGTATGGAGTTTTTTTAGTTCTGCTCAGCTCCACAGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtggtg                                  <  1:2377818/148‑1 (MQ=255)                                                                                                                aTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                              <  1:1596299/138‑1 (MQ=255)                                                                                                                aTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTtctc                                              >  2:1596299/1‑138 (MQ=255)                                                                                                                 tCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGc                                 >  1:2075017/1‑150 (MQ=255)                                                                                                                     aaTGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCggg                                                      >  1:1949857/1‑125 (MQ=255)                                                                                                                     aaTGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCggg                                                      <  2:1949857/125‑1 (MQ=255)                                                                                                                                aCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                  <  1:1448945/150‑1 (MQ=255)                                                                                                                                aCAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTc                  >  1:1897536/1‑150 (MQ=255)                                                                                                                                  agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGtt                    >  1:1550063/1‑146 (MQ=255)                                                                                                                                  agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGtt                    <  2:1550063/146‑1 (MQ=255)                                                                                                                                  agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg                <  1:321086/150‑1 (MQ=255)                                                                                                                                  agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg                <  2:2075017/150‑1 (MQ=255)                                                                                                                                  agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg                >  1:2129851/1‑150 (MQ=255)                                                                                                                                     gggCTCAGGGCGGGGGGGCGGCCCTGGGCCTGGTCAGTTCTGCTCAGCTACGCGGGACTTGCCTGGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGGGCC‑GGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgcgc            <  2:379866/149‑1 (MQ=255)                                                                                                                                     agGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAAACGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGgcg             >  2:916725/1‑150 (MQ=255)                                                                                                                                       gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGttt                                                                                         <  1:1417124/72‑1 (MQ=255)                                                                                                                                       gCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGttt                                                                                         >  2:1417124/1‑72 (MQ=255)                                                                                                                                        cTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTc          <  1:712051/150‑1 (MQ=255)                                                                                                                                        cTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAACTAATGGGTTCTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTc          >  2:1532946/1‑150 (MQ=255)                                                                                                                                         tCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTcc         >  1:215446/1‑149 (MQ=255)                                                                                                                                         tCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa         >  2:1029781/1‑150 (MQ=255)                                                                                                                                         tCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa         >  1:2421983/1‑150 (MQ=255)                                                                                                                                         tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCaa                                                                             <  2:2128110/82‑1 (MQ=255)                                                                                                                                         tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCaa                                                                             >  1:2128110/1‑82 (MQ=255)                                                                                                                                         tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCGGGCGCTCa         >  2:1309870/1‑150 (MQ=255)                                                                                                                                         tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa         >  1:1913082/1‑150 (MQ=255)                                                                                                                                         tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa         >  2:798435/1‑150 (MQ=255)                                                                                                                                         tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa         <  1:2164061/150‑1 (MQ=255)                                                                                                                                           aGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                    >  1:751685/1‑105 (MQ=255)                                                                                                                                           aGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTa                                                    <  2:751685/105‑1 (MQ=255)                                                                                                                                            gggCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtg                                     >  2:1494435/1‑119 (MQ=255)                                                                                                                                            gggCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAgtg                                     <  1:1494435/119‑1 (MQ=255)                                                                                                                                             ggCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCAc        <  2:2341318/147‑1 (MQ=255)                                                                                                                                             ggCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCAc        >  1:2341318/1‑147 (MQ=255)                                                                                                                                                cGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTg  >  2:2247043/1‑150 (MQ=255)                                                                                                                                                cGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTg  >  2:1742839/1‑150 (MQ=255)                                                                                                                                                    gagCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCt                                                                                                   >  1:1873133/1‑49 (MQ=255)                                                                                                                                                    gagCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCt                                                                                                   <  2:1873133/49‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTACTG  >  NC_002947/1376181‑1376473

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