Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,070,263	G→C	10.7%	intergenic (+63/+79)	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,263	0	G	C	10.7%	164.6 / 10.3	56	intergenic (+63/+79)	aroP‑I/PP_0928	aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base C (4/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

ATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACA  >  NC_002947/1070117‑1070408                                                                                                                                                        |                                                                                                                                                  aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCgcgg                                                                                                                                                   <  4:144131/149‑1 (MQ=255) aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCgcgg                                                                                                                                                   >  6:302145/1‑149 (MQ=255)    tcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCgcgggcg                                                                                                                                                >  2:17964/1‑149 (MQ=255)    tcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCAAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGCG‑‑AT‑CGCGCGCCGCgcgggcg                                                                                                                                                >  4:39623/1‑149 (MQ=255)      ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGggcggc                                                                                                                                              <  4:330414/149‑1 (MQ=255)          aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCt                                                                                                                                           <  1:17964/148‑1 (MQ=255)             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                        <  8:266451/148‑1 (MQ=255)             cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGCATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                        <  1:287040/148‑1 (MQ=255)                 tCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTca                                                                                                                                   <  8:110692/149‑1 (MQ=255)                    cGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCAGATCGGGTTTGTTCGTAAGGGGGATGCCGTTAGGTTTTCGGCGCCTGT‑‑GACATACGAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                                 <  1:44230/149‑1 (MQ=38)                                ctgtgctttgGGTGGCTTTCAAGCTGAAGAATGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                     <  3:4283/132‑1 (MQ=38)                                ctgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCATGTGAG‑‑AT‑CGAGCGCCGCGCGGGCTGAGCTCGAtct                                                                                                                                     >  4:4283/1‑132 (MQ=38)                                 tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGcc                                                                                                                    >  8:101714/1‑148 (MQ=255)                                         gggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTACGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                           <  1:105612/149‑1 (MQ=255)                                          ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                       <  6:25221/136‑1 (MQ=255)                                          ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                       >  5:25221/1‑136 (MQ=255)                                          ggtggCCTTCAAGCCGAACAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTCAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                           <  3:39623/148‑1 (MQ=38)                                                tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                       >  5:70784/1‑130 (MQ=255)                                                tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGaaa                                                                                                                       <  6:70784/130‑1 (MQ=255)                                                tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                                                    >  3:22049/1‑149 (MQ=255)                                                tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                                                    >  7:84649/1‑149 (MQ=255)                                                     aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                                               >  4:276152/1‑149 (MQ=255)                                                      gCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGCCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGcc                                                                                              >  4:275598/1‑149 (MQ=255)                                                           agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGa                                                                                          <  8:84649/148‑1 (MQ=255)                                                           agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                         >  2:192835/1‑149 (MQ=255)                                                            gaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAcg                                                                                        >  5:273832/1‑149 (MQ=255)                                                                 cTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTc                                                                                   <  7:150094/149‑1 (MQ=255)                                                                   cGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                                                                 <  3:276152/149‑1 (MQ=255)                                                                    gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                                                                                >  7:186360/1‑149 (MQ=255)                                                                    gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcc                                                                                  <  1:177625/147‑1 (MQ=255)                                                                     ccAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCt                                                                               <  6:273832/149‑1 (MQ=255)                                                                       aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                                                             <  7:101714/149‑1 (MQ=255)                                                                       aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                                                             >  4:306097/1‑149 (MQ=255)                                                                            aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTcc                                                                            <  1:350548/145‑1 (MQ=255)                                                                            aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGAGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                                                        >  4:331286/1‑149 (MQ=255)                                                                                gCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCt                                                                    <  5:33098/149‑1 (MQ=255)                                                                                    aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAt                                                                >  6:70782/1‑149 (MQ=255)                                                                                           ttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCt                                                         >  8:325732/1‑149 (MQ=255)                                                                                           ttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCt                                                         >  1:201869/1‑149 (MQ=255)                                                                                              ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                                                      >  6:68341/1‑149 (MQ=255)                                                                                                 gTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                    <  7:325732/148‑1 (MQ=255)                                                                                                    aGGGGTATGCCGTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                                                 >  4:249538/1‑148 (MQ=255)                                                                                                      acgcaTGCCTTTAGGCTTTCATCGCCTGTGAG‑‑AT‑CGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                  >  1:219317/5‑145 (MQ=17)                                                                                                      acgcaTGCCTTTAGGCTTTCATCGCCTGTGAG‑‑AT‑CGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                  <  2:219317/141‑1 (MQ=17)                                                                                                              ctttAGGCTTTCATCGCCTGTGAG‑‑AT‑CGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                      >  7:242821/3‑149 (MQ=25)                                                                                                               gTTAGGTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtgg                                     >  7:195572/1‑149 (MQ=255)                                                                                                                   ggTTTTCGGCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGCCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                                 <  3:275598/149‑1 (MQ=255)                                                                                                                          ggCGCCTGTGAG‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaa                           <  4:22049/148‑1 (MQ=255)                                                                                                                                tgtgag‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCg                      >  3:174011/6‑147 (MQ=255)                                                                                                                                  tgag‑‑AT‑CGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                                                          >  8:339009/4‑109 (MQ=255)                                                                                                                                  tgag‑‑AT‑CGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                                                          <  7:339009/106‑1 (MQ=255)                                                                                                                                        aT‑CGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACAACGTAGCGAg               >  6:187510/1‑148 (MQ=255)                                                                                                                                         t‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCa             >  1:134324/1‑149 (MQ=255)                                                                                                                                         t‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCa             <  5:302145/149‑1 (MQ=255)                                                                                                                                            gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAcc           >  7:300839/1‑149 (MQ=255)                                                                                                                                             aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                                                      <  3:248187/89‑1 (MQ=255)                                                                                                                                             aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGccc                                                                      >  4:248187/1‑89 (MQ=255)                                                                                                                                             aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCg          <  2:134324/149‑1 (MQ=255)                                                                                                                                              gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCgg         <  8:186360/149‑1 (MQ=255)                                                                                                                                                 ccgcccgcGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAg      >  1:156585/1‑149 (MQ=255)                                                                                                                                                 ccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAg      >  8:16910/1‑149 (MQ=255)                                                                                                                                                  cgcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGt     <  5:68341/149‑1 (MQ=255)                                                                                                                                                     gcgggcggCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCTTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACa  <  5:70782/149‑1 (MQ=255)                                                                                                                                                     gcgggcggCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACa  <  8:195572/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  ATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCTGT‑‑G‑‑AT‑CGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTACA  >  NC_002947/1070117‑1070408

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