Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,085,006	G→C	5.2%	R277R (CGG→CGC)	PP_0939 →	carbon‑nitrogen hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,085,006	0	G	C	5.2%	255.6 / 2.9	77	R277R (CGG→CGC)	PP_0939	carbon‑nitrogen hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base C (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

TCGACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCATCGATGCCGCCGAC  >  NC_002947/1084861‑1085146                                                                                                                                                    |                                                                                                                                               tcgACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGc                                                                                                                                               >  1:306718/1‑147 (MQ=255) tcgACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGc                                                                                                                                               <  2:306718/147‑1 (MQ=255) tcgACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCt                                                                                                                                             <  6:290030/149‑1 (MQ=255)   gACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTgc                                                                                                                                           <  7:137726/149‑1 (MQ=255)   gACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGC‑CGCAACGCGCCCTgc                                                                                                                                           <  7:79654/149‑1 (MQ=255)       cTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCa                                                                                                                                        >  3:112611/1‑148 (MQ=255)       cTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCAt                                                                                                                                       <  8:103666/149‑1 (MQ=255)                 cATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGt                                                                                                                            <  5:350563/149‑1 (MQ=255)                  aTTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTTGCGAGCGCGACATTGCAGAACAAGCGTCCATTCGGACGCGCATGGCGGTGGTCTTCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGt                                                                                                                            <  1:103945/149‑1 (MQ=255)                     ctcGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGAc                                                                                                                         <  6:87098/148‑1 (MQ=255)                        gCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTAt                                                                                                                      >  8:238373/1‑148 (MQ=255)                        gCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTAt                                                                                                                      <  1:2973/148‑1 (MQ=255)                         cAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGa                                                                                                                    <  2:107637/149‑1 (MQ=255)                              cAGGCGCAAGGTGAAGCGGTACTGCTTGCAGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGAATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGGAGGACTTTCGGACAGCGCATA‑ACTACGAGTGACTATGAGCCa                                                                                                                >  1:315986/1‑149 (MQ=255)                                 gcgcAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            >  6:376998/1‑149 (MQ=255)                                 gcgcAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            <  6:301340/149‑1 (MQ=255)                                 gcgcAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            >  2:284610/1‑149 (MQ=255)                                   gcAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGGCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGtt                                                                                                          <  6:308201/149‑1 (MQ=255)                                     aaGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTAt                                                                                                        >  2:92492/1‑149 (MQ=255)                                     aaGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTAt                                                                                                        <  7:238373/149‑1 (MQ=255)                                          gAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCAcc                                                                                                   >  3:71152/1‑149 (MQ=255)                                            aGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑t                                                                                                <  6:370318/149‑1 (MQ=255)                                             gCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑t                                                                                                >  3:143731/1‑148 (MQ=255)                                              cGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCa                                                                                              <  4:143731/149‑1 (MQ=255)                                                gTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGt                                                                                            >  2:303380/1‑149 (MQ=255)                                                gTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGt                                                                                            >  7:351976/1‑149 (MQ=255)                                                  aCTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGt                                                                                            >  8:50240/1‑147 (MQ=255)                                                    tGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGag                                                                                         >  6:33177/1‑148 (MQ=255)                                                     gCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑TGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGag                                                                                         >  8:53036/1‑147 (MQ=255)                                                      cTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATACCCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGag                                                                                         <  1:264144/147‑1 (MQ=255)                                                       ttGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGagcagc                                                                                     >  8:269951/1‑149 (MQ=255)                                                            gAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTg                                                                                >  4:105682/1‑149 (MQ=255)                                                             aGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTgg                                                                               <  5:33177/149‑1 (MQ=255)                                                              gcgcGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGc                                                                              <  7:50240/149‑1 (MQ=255)                                                              gcgcGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGc                                                                              >  2:207670/1‑149 (MQ=255)                                                                  gACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            <  5:39770/116‑1 (MQ=255)                                                                  gACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            >  6:39770/1‑116 (MQ=255)                                                                        gaagaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAggcggc                                                                    >  7:184369/1‑149 (MQ=255)                                                                        gaagaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAggcggc                                                                    >  6:53060/1‑149 (MQ=255)                                                                        gaagaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAggcgg                                                                     <  3:105682/148‑1 (MQ=255)                                                                         aagaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑AGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TAAGTGAGCAGCTACTGGCCCCAggcggat                                                                   >  7:143446/1‑147 (MQ=255)                                                                          agaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCtt                                                                  <  8:351976/149‑1 (MQ=255)                                                                           gaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGAt                                                                                                             >  1:99789/1‑106 (MQ=255)                                                                           gaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGAt                                                                                                             <  2:99789/106‑1 (MQ=255)                                                                           gaaCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTg                                                                 >  7:43865/1‑149 (MQ=255)                                                                             aCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGAt                                                                                                             >  5:218947/1‑104 (MQ=255)                                                                             aCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGAt                                                                                                             <  6:218947/104‑1 (MQ=255)                                                                             aCAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGAc                                                               >  3:359162/1‑149 (MQ=255)                                                                               aaGCGTCCATTCTGACGCTCATGCCGGTGGT‑ATCGCACCGGAGCTTTTTATCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGAccc                                                             <  2:315986/149‑1 (MQ=255)                                                                                 gCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTg                                                           <  1:92492/149‑1 (MQ=255)                                                                                 gCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTg                                                           >  6:87105/1‑149 (MQ=255)                                                                                  cGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTgg                                                          <  5:366648/149‑1 (MQ=255)                                                                                    tCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGAc                                                        <  8:184369/149‑1 (MQ=255)                                                                                       aTTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGc                                                     >  7:159400/1‑149 (MQ=255)                                                                                           ggACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑TGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTg                                                  <  7:53036/148‑1 (MQ=255)                                                                                            gaagcgcATGCAGGTGGT‑ATCGCCCCGGCGATTTTTTTCGAAGTG‑CGCGTTGCTGCCTGAGCATA‑CCTAGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCa                                                <  8:143446/146‑1 (MQ=255)                                                                                            gACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCa                                                <  4:71152/149‑1 (MQ=255)                                                                                              cgcgCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCaga                                              <  5:376998/149‑1 (MQ=255)                                                                                               gcgcATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCGTTCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCaga                                              <  1:284610/149‑1 (MQ=255)                                                                                                cgcATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGGCCCTGGACAGCCTGCAGAGc                                            <  4:112611/149‑1 (MQ=255)                                                                                                       cGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCa                                                                                              >  8:39356/1‑92 (MQ=255)                                                                                                         gtggt‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGt                                   <  1:207670/149‑1 (MQ=255)                                                                                                               aTCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCa                                                                                              >  6:107429/1‑85 (MQ=255)                                                                                                                tCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTtg                             >  3:133380/1‑149 (MQ=255)                                                                                                                tCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTtg                             >  3:127864/1‑149 (MQ=255)                                                                                                                  gCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTtggc                           <  6:18678/149‑1 (MQ=255)                                                                                                                  gCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTtgg                            <  8:184631/148‑1 (MQ=255)                                                                                                                            tttttttCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CATCGGAGTGACTATGAGCCAGATGTTATCCTCGG‑TCAGTGAGCAGATCCTGGCACCAGGCGGCTTGACCCTGGACAGCCTGAAGAGAGTGCTGGGTGAGTtgtctggcacacg                  >  1:57907/1‑143 (MQ=255)                                                                                                                            agatatgcGCAGGG‑CGCCTCGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGc                 <  3:334844/142‑1 (MQ=255)                                                                                                                                 ttCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCATCGa            >  6:148041/1‑149 (MQ=255)                                                                                                                                 tgcGCAGGC‑CGCAACGCGCCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTAt                                                                                                        <  1:43338/56‑1 (MQ=18)                                                                                                                                 tgcGCAGGC‑CGCAACGCGCCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTAt                                                                                                        >  2:43338/3‑58 (MQ=18)                                                                                                                                   cGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCAt               >  7:213509/1‑144 (MQ=255)                                                                                                                                   cGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCAt               <  8:213509/144‑1 (MQ=255)                                                                                                                                   cGCAGGC‑CGCAACGCGCCCTGCGCATA‑AATCGGAGTGACTATGAGCCAGAGGTTATACAGAG‑TCAGTGAGCGGCTCCTGGCCCCCGGCGGCTTGACGCTGGACAGCCTGCAGAGCGTGCTGGGGGAGTTGGCGGGCCCCGGCaagaagg          >  4:334844/1‑143 (MQ=255)                                                                                                                                      aGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCATCGATgccg       <  7:269951/149‑1 (MQ=255)                                                                                                                                            gcgTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCATCGATGCCGCCGAc  <  5:87105/149‑1 (MQ=255)                                                                                                                                              aacgcgcCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCg                                           <  3:201831/103‑1 (MQ=255)                                                                                                                                              aacgcgcCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCg                                           >  4:201831/4‑106 (MQ=255)                                                                                                                                                   ggCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGt                                          <  2:105958/102‑1 (MQ=255)                                                                                                                                                   ggCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGt                                          >  1:105958/1‑102 (MQ=255)                                                                                                                                                    gCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAggcggc                                                                    <  4:205994/75‑1 (MQ=255)                                                                                                                                                    gCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAggcggc                                                                    >  3:205994/1‑75 (MQ=255)                                                                                                                                                    |                                                                                                                                               TCGACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCATCGATGCCGCCGAC  >  NC_002947/1084861‑1085146

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