Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,085,002	T→C	6.4%	L276S (TTG→TCG) ‡	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,002	0	T	C	6.4%	233.9 / 3.3	78	L276S (TTG→TCG) ‡	PP_0939	carbon‑nitrogen hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (35/38);  new base C (2/3);  total (37/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 = 1.00e+00

GCGATCGTCGACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCATCGATGCCGCCGAC  >  NC_002947/1084854‑1085146                                                                                                                                                       |                                                                                                                                                   gcgATCGTCGACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAGGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGtt                                                                                                                                                    <  2:57358/149‑1 (MQ=255) gcgATCGTCGACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGtt                                                                                                                                                    >  2:2973/1‑149 (MQ=255)   gATCGTCGACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGc                                                                                                                                                  >  5:370318/1‑149 (MQ=255)        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                                                                                                                      <  1:2973/148‑1 (MQ=255)                               gCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTAt                                                                                                                      >  8:238373/1‑148 (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                                                                                                            >  2:284610/1‑149 (MQ=255)                                        gcgcAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            >  6:376998/1‑149 (MQ=255)                                        gcgcAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            <  6:301340/149‑1 (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                                                                              >  2:207670/1‑149 (MQ=255)                                                                     gcgcGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGc                                                                              <  7:50240/149‑1 (MQ=255)                                                                         gACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            >  6:39770/1‑116 (MQ=255)                                                                         gACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATg                                                                                                            <  5:39770/116‑1 (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                                                                                                        >  2:43338/3‑58 (MQ=18)                                                                                                                                        tgcGCAGGC‑CGCAACGCGCCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTAt                                                                                                        <  1:43338/56‑1 (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)                                                                                                                                                       |                                                                                                                                                   GCGATCGTCGACCCTTGGGGGCGCATTCTCGCAGAACAGGCGCAAGGTGAAGCGGTACTGCTTGCCGAGCGCGACATTGAAGAACAAGCGTCCATTCGGACGCGCATGCCGGTGGT‑ATCGCACCGGCGCTTTTTTTCGCAGGG‑CGCGTTGCGGCCTGCGCATA‑CCTCGGAGTGACTATGAGCCAGATGTTATCCACCG‑TCAGTGAGCAGCTCCTGGCCCCAGGCGGCTTGACCCTGGACAGCCTGCAGAGCGTGCTGGGTGAGTTGGCTGGCCCCGGCATCGATGCCGCCGAC  >  NC_002947/1084854‑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: ‑