Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,333,363	T→C	8.7%	Q148R (CAG→CGG)	PP_2931 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,333,363	0	T	C	8.7%	191.6 / 7.4	69	Q148R (CAG→CGG)	PP_2931	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/29);  new base C (2/4);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

ATGACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAATCGGGCAGCA  >  NC_002947/3333219‑3333510                                                                                                                                                 |                                                                                                                                                    aTGACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGctgctg                                                                                                                                                 <  8:268889/149‑1 (MQ=255)   gACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGctgctg                                                                                                                                                 >  8:215266/1‑147 (MQ=255)    aCGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTgggg                                                                                                                                              <  1:17572/149‑3 (MQ=255)    aCGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTggg                                                                                                                                               <  8:69700/148‑2 (MQ=255)      gcgcACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTgggg                                                                                                                                              >  4:241616/1‑145 (MQ=255)      gcgcACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCAGCTGCAGCCGCTgggg                                                                                                                                              <  3:241616/147‑3 (MQ=255)               cTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATgg                                                                                                                                   >  8:15915/1‑149 (MQ=255)                  gACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGgca                                                                                                                                 <  7:215266/148‑1 (MQ=255)                     cTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCAt                                                                                                                                     <  8:233193/141‑1 (MQ=255)                     cTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCAt                                                                                                                                     >  7:233193/1‑141 (MQ=255)                     cTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCAt                                                                                                                             >  5:79232/1‑149 (MQ=255)                        ggCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATc                                                                                                                            <  5:199481/147‑1 (MQ=255)                         gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAgctgc                                                                                                                                                   <  6:282734/123‑1 (MQ=255)                         gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAg                                                                                                                         <  2:239282/149‑1 (MQ=255)                         gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAg                                                                                                                         <  2:157432/149‑1 (MQ=255)                         gCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAg                                                                                                                         <  2:136584/149‑1 (MQ=255)                                     aTCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGt                                                                                                              >  1:220027/1‑148 (MQ=255)                                     aTCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGt                                                                                                              <  2:220027/148‑1 (MQ=255)                                      tCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCg                                                                                                            >  3:59303/1‑149 (MQ=255)                                      tCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGACGGTGacgacg                                                                                                            >  8:345983/1‑149 (MQ=255)                                        cGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGgt                                                                                                          <  7:16336/149‑1 (MQ=255)                                         gCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGt                                                                                                                   <  5:354918/139‑1 (MQ=255)                                         gCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGt                                                                                                                   >  6:354918/1‑139 (MQ=255)                                            ccTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGtt                                                                                                       >  3:120578/1‑148 (MQ=255)                                              tGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGAc                                                                                                    >  3:342435/1‑149 (MQ=255)                                                  cgctgcgctgCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGcac                                                                                                <  7:137916/149‑1 (MQ=255)                                                       cgctgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCt                                                                                           >  6:38254/1‑149 (MQ=255)                                                          tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGAc                                                                                                    <  1:313884/137‑1 (MQ=255)                                                          tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGAc                                                                                                    >  2:313884/1‑137 (MQ=255)                                                          tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGc                                                                                        <  4:2483/149‑1 (MQ=255)                                                          tgcAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGc                                                                                        >  7:189113/1‑149 (MQ=255)                                                               ggcggcCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACg                                                                                   >  3:222450/1‑149 (MQ=255)                                                                  ggcCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGt                                                                                >  3:181884/1‑149 (MQ=255)                                                                   gcCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTg                                                                               <  4:342435/149‑1 (MQ=255)                                                                    ccTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTgg                                                                              >  5:242272/1‑149 (MQ=255)                                                                     cTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGgtg                                                                                                         >  2:354966/1‑121 (MQ=255)                                                                     cTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGgtg                                                                                                         <  1:354966/121‑1 (MQ=255)                                                                          tgatCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCgctg                                                                        <  7:15915/149‑1 (MQ=255)                                                                             tctcGCTCATCAGCAGGTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGctgctg                                                                     >  3:317320/1‑149 (MQ=255)                                                                             tctcGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGctgctg                                                                     >  6:288754/1‑149 (MQ=255)                                                                               tcGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGt                                                                    >  8:296435/1‑148 (MQ=255)                                                                                 gCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCAc                                                                 <  4:222450/149‑1 (MQ=255)                                                                                         gcagcTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGAtt                                                         >  6:378011/1‑149 (MQ=255)                                                                                             cTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCg                                                     >  5:55941/1‑149 (MQ=255)                                                                                              tCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGt                                                    >  1:328646/1‑149 (MQ=255)                                                                                               ccATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTg                                                   <  4:59303/149‑1 (MQ=255)                                                                                               ccATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTg                                                   <  7:296435/149‑1 (MQ=255)                                                                                               ccATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGAAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGACGCTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCg                                                                           >  2:316816/1‑125 (MQ=255)                                                                                                cATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTg                                                   <  4:120578/148‑1 (MQ=255)                                                                                                       gCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCg                                           >  4:358069/1‑149 (MQ=255)                                                                                                       gCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCg                                           >  6:333799/1‑149 (MQ=255)                                                                                                           gTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGc                                       >  5:144889/1‑149 (MQ=255)                                                                                                           gTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGc                                       >  3:48348/1‑149 (MQ=255)                                                                                                                ttGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCg                                  >  7:314561/1‑149 (MQ=255)                                                                                                                   aCCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGc                                <  6:55941/148‑1 (MQ=255)                                                                                                                     cGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCtt                             >  7:182588/1‑149 (MQ=255)                                                                                                                        aGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGa                          <  8:189113/149‑1 (MQ=255)                                                                                                                        aGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGa                          <  4:317320/149‑1 (MQ=255)                                                                                                                          gTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGc                        >  7:140531/1‑149 (MQ=255)                                                                                                                                cccccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTg                  <  4:181884/149‑1 (MQ=255)                                                                                                                                cccccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTg                  >  2:194380/1‑149 (MQ=255)                                                                                                                                  cccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTggg                <  2:328646/149‑1 (MQ=255)                                                                                                                                  cccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTggg                <  5:288754/149‑1 (MQ=255)                                                                                                                                  cccAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTggg                <  6:41754/149‑1 (MQ=255)                                                                                                                                   ccAGCAGCTGCAGCCGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGc                                >  4:28061/1‑132 (MQ=255)                                                                                                                                   ccAGCAGCTGCAGCCGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGc                                <  3:28061/132‑1 (MQ=255)                                                                                                                                      gCGGCTGCAGCTGCTGGAGAAAGTTTATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCGCTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTaa            >  8:103567/1‑149 (MQ=255)                                                                                                                                           tGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAATCggg       >  5:179700/1‑149 (MQ=255)                                                                                                                                           tGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAATCggg       <  8:314561/149‑1 (MQ=255)                                                                                                                                            gCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTgg                 >  3:19540/1‑138 (MQ=255)                                                                                                                                            gCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTgg                 <  4:19540/138‑1 (MQ=255)                                                                                                                                             cAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAATCGGgca     <  7:345983/149‑1 (MQ=255)                                                                                                                                               gctgctGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAATCGGgcagc   <  6:242272/149‑1 (MQ=255)                                                                                                                                                ctgctgGAGAAAGTTTATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCGCTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAATCGGgcagca  <  7:103567/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    ATGACGCGCACGGCCTCGACCTTGGCCTGTTTGTTCATCCGCCCCTGGCCGCTGCGCTGCAGGGCGGCCTGGATGATCTCGCTCATCAGCAGCTCCATGTCGGCAGGTTCGTTGACCGCAGGTATGGCCCCCAGCGGCTGCAGCTGCTGGAGAAAGTTCATGGCAGCATCCAGGCCGGTGACGTCGGTGTTGACGCACAGGCTGGCGAACGGGTGGCCGCTGCTGTCACGGAAGATTGCCGTGGCACTTCGTAGCGTGCGGCCTTTGAGCGTAGTGGGGTAATCGGGCAGCA  >  NC_002947/3333219‑3333510

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