Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,070,286	A→C	6.9%	intergenic (+86/+56)	aroP‑I → / ← PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,070,286	0	A	C	6.9%	180.0 / 5.0	58	intergenic (+86/+56)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base C (2/2);  total (27/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 = 6.49e-01

TGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAA  >  NC_002947/1070141‑1070422                                                                                                                                                         |                                                                                                                                         tGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGAt                                                                                                                                         >  4:226367/1‑149 (MQ=255) tGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGAt                                                                                                                                         <  6:75657/149‑1 (MQ=255)   gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATg                                                                                                                                        >  1:201969/1‑148 (MQ=255)   gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATg                                                                                                                                        <  3:63455/148‑1 (MQ=255)    aTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGa                                                                                                                                       >  8:278510/1‑148 (MQ=255)     tCGCTGTGCTGGGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGaaa                                                                                                                                     <  3:260613/149‑1 (MQ=255)      cGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGaaa                                                                                                                                     <  6:165353/148‑1 (MQ=255)           tgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTaa                                                                                                                               <  1:23119/149‑1 (MQ=255)             ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAg                                                                                                                             <  7:11051/149‑1 (MQ=255)             cggtgggtggCTTTCAAGCTGAAGAAGGCGCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTACGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑AGCTCGCGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAg                                                                                                                             <  5:211259/147‑1 (MQ=38)                tgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCa                                                                                                                          >  1:301146/1‑149 (MQ=255)                 gggtggCTTTCACGCTGAAGAAGGCTCGCCAGGCCAAGGGCGGAGCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑GACA‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                                         <  1:312761/149‑1 (MQ=38)                  ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                                         <  3:174333/148‑1 (MQ=255)                   gtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                                                                       <  4:264018/149‑1 (MQ=255)                   gtggCTTTCAAGCTGAAGAAGGCTCGCCAAGACAAGGGCTGAACGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTATTCGGCGCCTGTGAG‑‑A‑‑ACGAGCGCCGCGCGGGCGGCGCAA‑GATCTCACAGGCGATGAAAACCTAAAgggaag                                                                                                                        >  7:198909/1‑144 (MQ=37)                   gtggCTTTCAACCTTAATAAGGCTCGCCAGGCCTTGGGCTTATCTGGTTTTTTCGTTTGTGCTTTGCCGTTAGCTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCCCCGCGCGGGGGGGGCTC‑GTTCTCACAGGCGATGAAAGCCTAAAGGCATg                                                                                                                        <  8:198909/148‑1 (MQ=255)                        tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTaa                                                                                                                               <  4:268694/136‑1 (MQ=255)                        tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTaa                                                                                                                               >  3:268694/1‑136 (MQ=255)                           cAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGGTTGTTCGTAAGGGGTAGGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                                                                <  4:212107/148‑1 (MQ=38)                            aaGCTGAAGCAGGCTCGCCAGGCCAAGGGCTGATCGGGGTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcgg                                                                                                              <  8:219464/149‑1 (MQ=38)                                tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                                                                    <  3:92997/139‑1 (MQ=255)                                tGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                                                                    >  4:92997/1‑139 (MQ=255)                                     aagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcgc                                                                                                     >  6:16325/1‑149 (MQ=255)                                            gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                                                                                              >  3:197250/1‑149 (MQ=255)                                            gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                                                                               <  2:201969/148‑1 (MQ=255)                                                ggCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                                                          <  3:226367/149‑1 (MQ=255)                                                 gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                                                                >  6:274002/1‑126 (MQ=255)                                                 gCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggc                                                                                                                <  5:274002/126‑1 (MQ=255)                                                    aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAg                                                                                                                             <  8:282139/110‑1 (MQ=255)                                                    aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAg                                                                                                                             >  7:282139/1‑110 (MQ=255)                                                                ggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                                           <  4:197250/148‑1 (MQ=255)                                                                         gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                                  >  8:151327/1‑148 (MQ=255)                                                                         gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                                  >  3:93400/1‑148 (MQ=255)                                                                         gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                                  >  3:88989/1‑148 (MQ=255)                                                                          tAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTTAG‑‑A‑‑TAGAGCGCCCCGCGGGCGGCGCTC‑GGTATCACAGGGGATGAAAGCCTAAAGGCATGCGTCGTGCGGCCCGGGCGCGGTGCCCGTGCAAGGCCCGCAATGCGGCCTCCGCCaa                                                                >  4:169396/1‑149 (MQ=255)                                                                              gggTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcac                                                            <  7:278510/149‑1 (MQ=255)                                                                                     cctttAGGCTTTCATCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgt                                                     >  3:91721/4‑149 (MQ=18)                                                                                      cGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                                    <  4:93400/149‑1 (MQ=255)                                                                                      cGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                                    <  4:88989/149‑1 (MQ=255)                                                                                      cGTTAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                                    <  1:264966/149‑1 (MQ=255)                                                                                        ttAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca                                                         >  4:233615/1‑142 (MQ=255)                                                                                        ttAGGTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca                                                         <  3:233615/142‑1 (MQ=255)                                                                                          aGGCTTTCATCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggtg                                                 <  4:91721/148‑1 (MQ=18)                                                                                           ggTTTTCGGCGCCTGTGAG‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                                               <  7:95705/149‑1 (MQ=255)                                                                                                       ctgtgag‑‑A‑‑TCGAGCGCCGCCCGCGCGGCGCTCGGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCg                                    >  1:38018/7‑149 (MQ=38)                                                                                                            ag‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGc                                >  2:189213/2‑147 (MQ=255)                                                                                                                a‑‑TCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAg                             >  1:156115/1‑148 (MQ=255)                                                                                                                   tCGAGCGCCGCCCGCGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCa                           >  6:181063/1‑149 (MQ=255)                                                                                                                     gAGCGCCGCGCGGGCGGCGCTT‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCCCGCCCAGTGGTGGACATCCACACCGTAGCGAGCAcc                         >  2:85815/1‑149 (MQ=255)                                                                                                                                gggcggCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACa                >  6:230042/1‑147 (MQ=255)                                                                                                                                 ggcggcGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAg             <  5:230042/149‑1 (MQ=255)                                                                                                                                  gcggcgCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCGTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAg             >  3:1699/1‑148 (MQ=255)                                                                                                                                       gcTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAg        >  3:39375/1‑148 (MQ=255)                                                                                                                                            gATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGc                                >  6:98781/1‑120 (MQ=255)                                                                                                                                            gATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGc                                <  5:98781/120‑1 (MQ=255)                                                                                                                                                  acaGGCGCTGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGa                              <  6:276883/116‑1 (MQ=255)                                                                                                                                                  acaGGCGCTGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGa                              >  5:276883/1‑116 (MQ=255)                                                                                                                                                  acaGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCaa  <  7:151327/144‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                         TGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑A‑‑TCGAGCGCCGCGCGGGCGGCGCTC‑GATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACATAGAACAGCATCAA  >  NC_002947/1070141‑1070422

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