Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	614,316	Δ1 bp	6.4%	intergenic (‑132/+46)	dxs ← / ← ispA	1‑deoxyxylulose‑5‑phosphate synthase/farnesyl diphosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	614,316	0	T	.	6.4%	180.8 / 7.3	63	intergenic (‑132/+46)	dxs/ispA	1‑deoxyxylulose‑5‑phosphate synthase/farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base . (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

TGAAACGTCGTGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTG  >  NC_002947/614168‑614465                                                                                                                                                      |                                                                                                                                                      tGAAACGTCGTGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATc                                                                                                                                                      <  1:1976032/150‑1 (MQ=255) tGAAACGTCGTGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATc                                                                                                                                                      <  1:1769982/150‑1 (MQ=255)         cgtGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTcgc                                                                                                                                              <  2:1943476/150‑1 (MQ=255)           tGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCgg                                                                                                                                            <  2:954933/150‑1 (MQ=255)           tGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCgg                                                                                                                                            <  1:1476689/150‑1 (MQ=255)           tGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCgg                                                                                                                                            <  2:2447284/150‑1 (MQ=255)           ccggcctacaacgacgccggccaggacaatgcccgttcgacaccgggcgaaggatttggcgtgtatgccggggtttcctacagcttccgctgacgggtgggcgCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGA‑CGCATTCGCgg                                                                                                                                            >  2:1147372‑M1/104‑150 (MQ=255)                    ggCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTcgc                                                                                                                                              >  1:628405/1‑139 (MQ=255)                    ggCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTcgc                                                                                                                                              <  2:628405/139‑1 (MQ=255)                    aacgacgccggccaggacaatgcccgttcgacaccgggcgaaggatttggcgtgtatgccggggtttcctatagcttccgctgacgggtgggcgCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGA‑CGCATTcgc                                                                                                                                              <  1:2416850‑M1/45‑1 (MQ=255)                    aacgacgccggccaggacaatgcccgttcgacaccgggcgaaggatttggcgtgtatgccggggtttcctatagcttccgctgacgggtgggcgCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑TGACGCATTcgc                                                                                                                                              >  2:2416850‑M1/95‑139 (MQ=255)                            ggtataGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCgg                                                                                                                                            <  1:405537/131‑1 (MQ=255)                            gTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCgg                                                                                                                                            >  2:405537/1‑133 (MQ=255)                                           ccgttcgacaccgggcgaaggatttggcgtgtatgccggggtttcctacagcttccgctgacgggtgggcgCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑TGACGCATTCGCggg                                                                                                                                           >  2:1544882‑M1/72‑119 (MQ=255)                                           ccgttcgacaccgggcgaaggatttggcgtgtatgccggggtttcctacagcttccgctgacgggtgggcgCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑TGACGCATTCGCggg                                                                                                                                           <  1:1544882‑M1/48‑1 (MQ=255)                                           aGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCggg                                                                                                                                           <  2:564044/119‑1 (MQ=255)                                           aGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCggg                                                                                                                                           >  1:564044/1‑119 (MQ=255)                                            cgttcgacaccgggcgaaggatttggcgtgtatgccggggtttcctacagcttccgctgacgggtgggagCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGA‑CGCATTCGCgg                                                                                                                                            >  2:272417‑M1/71‑117 (MQ=255)                                            cgttcgacaccgggcgaaggatttggcgtgtatgccggggtttcctacagcttccgctgacgggtgggagCGGGTTTACCCGCGAATGCGTCGGCAAATTCAATGA‑CGCATTCGCgg                                                                                                                                            <  1:272417‑M1/47‑1 (MQ=255)                                                                                                                aTCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCt                                       >  2:74604/1‑150 (MQ=255)                                                                                                                aTCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCt                                       >  1:1689615/1‑150 (MQ=255)                                                                                                                 tCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAg                                                                                                                    <  1:1499571/72‑1 (MQ=255)                                                                                                                 tCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAg                                                                                                                    >  2:1499571/1‑72 (MQ=255)                                                                                                                 tCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCtt                                      >  1:2746288/1‑150 (MQ=255)                                                                                                                  CGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTg                                     >  2:2746288‑M1/1‑150 (MQ=255)                                                                                                                          CCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAg                                                                                                                    <  2:2087316‑M1/63‑1 (MQ=255)                                                                                                                          CCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAg                                                                                                                    >  1:2087316‑M1/1‑63 (MQ=255)                                                                                                                          CCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCgg                             <  2:2840305‑M1/150‑1 (MQ=255)                                                                                                                           CCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGtt                                                                                                    <  2:1352859‑M1/78‑1 (MQ=255)                                                                                                                           CCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGtt                                                                                                    >  1:1352859‑M1/1‑78 (MQ=255)                                                                                                                           CCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACACTATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGc                            <  1:1574034‑M1/150‑1 (MQ=255)                                                                                                                           CCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAAccc                                        >  1:1814098‑M1/1‑138 (MQ=255)                                                                                                                           CCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAAccc                                        <  2:1814098‑M1/138‑1 (MQ=255)                                                                                                                           CCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGc                            <  1:567670‑M1/150‑1 (MQ=255)                                                                                                                           CCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGc                            <  1:1960389‑M1/150‑1 (MQ=255)                                                                                                                             GCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACCATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCa                          <  2:416017‑M1/150‑1 (MQ=255)                                                                                                                             GCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCa                          <  2:1073823‑M1/150‑1 (MQ=255)                                                                                                                             GCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCa                          <  2:2479615‑M1/150‑1 (MQ=255)                                                                                                                                 ATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAgcgc                      <  2:139087‑M1/150‑1 (MQ=255)                                                                                                                                     GTCGGCAAATTCAA‑GATCGCAGTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTgg                  <  2:128469‑M1/150‑1 (MQ=255)                                                                                                                                      TCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGt                 <  2:1342842‑M1/150‑1 (MQ=255)                                                                                                                                       CGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTc                >  1:128469‑M1/1‑150 (MQ=255)                                                                                                                                        ggCAAATTCAA‑GATCGCATGCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTcg               >  2:2972774/1‑150 (MQ=255)                                                                                                                                         gCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAat                                                                                     <  1:2871529/79‑1 (MQ=255)                                                                                                                                         gCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAat                                                                                     >  2:2871529/1‑79 (MQ=255)                                                                                                                                           aaaTTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCgcgc            >  1:44193/1‑150 (MQ=255)                                                                                                                                           aaaTTCAA‑GATCGCATTCGCGGGGAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCgcgc            >  1:2587392/1‑150 (MQ=255)                                                                                                                                             aTTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCt                                                     <  2:2267423/107‑1 (MQ=255)                                                                                                                                             aTTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCt                                                     >  1:2267423/1‑107 (MQ=255)                                                                                                                                             aTTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAgc                        <  2:141474/136‑1 (MQ=255)                                                                                                                                             aTTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAgc                        >  1:141474/1‑136 (MQ=255)                                                                                                                                              ttCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTTTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGt         >  1:2521932/1‑150 (MQ=255)                                                                                                                                              ttCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGt         >  1:2700049/1‑150 (MQ=255)                                                                                                                                               tCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGtt        >  2:822504/1‑150 (MQ=255)                                                                                                                                               tCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGtt        >  2:720232/1‑150 (MQ=255)                                                                                                                                                 aa‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCtttt                                                  <  1:1873761/106‑1 (MQ=255)                                                                                                                                                 aa‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCtttt                                                  >  2:1873761/1‑106 (MQ=255)                                                                                                                                                 aa‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCg      >  1:527790/1‑150 (MQ=255)                                                                                                                                                 aa‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAAAGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCAGAACCCTTGCAGTGCGGCCAGCGCATGGTCGCGCCGTTCg      >  2:1195701/1‑150 (MQ=255)                                                                                                                                                    gATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGAt    >  1:59626/1‑150 (MQ=255)                                                                                                                                                    gATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGACGTCCGTCCAGCGCCTGGTCGCGCAGTTCGAt    >  1:1542719/1‑150 (MQ=255)                                                                                                                                                     aTCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGAtt   <  1:74604/150‑1 (MQ=255)                                                                                                                                                     aTCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGAtt   >  2:1958007/1‑150 (MQ=255)                                                                                                                                                      tCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGAt    <  2:409791/148‑1 (MQ=255)                                                                                                                                                      tCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGAt    >  1:409791/1‑148 (MQ=255)                                                                                                                                                      tCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTg  >  1:2222204/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGAAACGTCGTGGGCATCAGGCGAGTCGTTATAGGTGTGTGAAGACGCGGCAGTTTACCCCATTGTGGGAAAAGTGCCCATTCGACCACAGGTTACGCGGCCCCCTGTAGGATCGGGTTTACCCGCGAATGCGTCGGCAAATTCAA‑GATCGCATTCGCGGGTAAACCCGATCCTACAGGGATAGGTGAAGCTGTTTAGTTACGGCGTTCGACAATATAACGCGCCAAGGCCCGCAGCGGCTCGGCCTTTTCGCCGAACCCTTGCAGTGCGGCCAGCGCCTGGTCGCGCAGTTCGATTG  >  NC_002947/614168‑614465

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