Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,619,466	C→T	7.1%	T345T (ACG→ACA)	PP_4938 ←	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,619,466	0	C	T	7.1%	248.2 / 6.5	84	T345T (ACG→ACA)	PP_4938	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (41/37);  new base T (2/4);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

CCCTTGCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCA  >  NC_002947/5619320‑5619614                                                                                                                                                   |                                                                                                                                                      cccTTGCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCg                                                                                                                                                      >  6:110387/1‑148 (MQ=255)   cTTGCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTa                                                                                                                                                    >  1:79512/1‑148 (MQ=255)      gccgccTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAgc                                                                                                                                                >  7:218269/1‑149 (MQ=255)      gccgccTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAg                                                                                                                                                 <  5:60160/148‑1 (MQ=255)        cgccTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGgcgc                                                                                                                                              >  6:257272/1‑149 (MQ=255)           cTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTca                                                                                                                                           <  4:266467/149‑2 (MQ=255)              cGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                         <  5:257272/148‑4 (MQ=255)              cGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcag                                                                                                                                          <  1:210452/147‑3 (MQ=255)              cGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGa                                                                                                                                          <  5:311759/147‑1 (MQ=255)                  gCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGcc                                                                                                                                     >  6:128710/1‑148 (MQ=255)                  gCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCt                                                                                                                                    >  8:53250/1‑149 (MQ=255)                            tGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                         >  5:29217/1‑131 (MQ=255)                            tGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCTGTACGGCGCTcagt                                                                                                                                         <  6:29217/134‑4 (MQ=255)                             gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCTCCGTACAGCGCTCGATCGCCTCAGC‑CCAGtcg                                                                                                                        >  8:308952/1‑149 (MQ=255)                             gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATc                                                                                                                                        >  4:121201/1‑134 (MQ=255)                             gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATc                                                                                                                                        <  3:121201/134‑1 (MQ=255)                             gACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGgcg                                                                                                                        >  5:78462/1‑149 (MQ=255)                                 cGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCttt                                                                                                                    <  4:124196/149‑1 (MQ=255)                                   cTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTtcg                                                                                                                  >  6:122486/1‑149 (MQ=255)                                      aaGGTCAGCTGCGCAACAGGCGCTGCGAGCCGCTGCAGAACGCCGGCCAGGCCTGCTCCGGTTCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATTCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGAGCTttcttc                                                                                                                <  2:101025/148‑1 (MQ=255)                                       aGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCgg                                                                                                              <  5:48532/149‑1 (MQ=255)                                        tgTCAGCTGCGCATCAGGCTCTGCGAGGCGCTGCTGAACGGCTGCCAGGCCTGCTCCGGTGCCAGTGCCTTGCGCTGCTTAGCGGCGATTGTCGTTGGACTGAGCGCCGTAGAGCTGTCGATCGCCTCAGG‑CCAGGCGATTTCGTCgg                                                                                                              <  4:270599/147‑1 (MQ=255)                                             gCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCg                                                                                                        >  4:143438/1‑149 (MQ=255)                                              cTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCgg                                                                                                       <  5:122486/149‑1 (MQ=255)                                                 cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGc                                                                                                      <  5:110387/147‑1 (MQ=255)                                                 cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAg                                                                                                    <  4:265378/149‑1 (MQ=255)                                                 cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAg                                                                                                    <  4:260398/149‑1 (MQ=255)                                                 cgcAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAg                                                                                                    <  7:53250/149‑1 (MQ=255)                                                   cAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAgg                                                                                                   >  8:157247/1‑148 (MQ=255)                                                    aaCAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTa                                                                                                 <  2:79512/149‑1 (MQ=255)                                                         gCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAc                                                                                            <  5:128710/149‑1 (MQ=255)                                                         gCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAc                                                                                            <  8:218269/149‑1 (MQ=255)                                                         gCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAc                                                                                            <  4:69976/149‑1 (MQ=255)                                                          ctctGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAcc                                                                                           >  2:52049/1‑149 (MQ=255)                                                           tctGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCCGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAccc                                                                                          >  2:171709/1‑149 (MQ=255)                                                           tctGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAccc                                                                                          >  6:364601/1‑149 (MQ=255)                                                             tGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCgt                                                                                        >  6:13517/1‑149 (MQ=255)                                                               cGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCTACCTGCCGGTTGccgacgaaa                                                                                       <  6:78462/148‑8 (MQ=255)                                                               cGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCgt                                                                                        <  1:38049/147‑1 (MQ=255)                                                                gAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTc                                                                                     <  1:52049/149‑1 (MQ=255)                                                                  gCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTcgc                                                                                   <  8:187785/149‑1 (MQ=255)                                                                  gCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGAGAcgc                                                                                   >  4:61406/1‑149 (MQ=255)                                                                          aGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGc                                                                           >  6:1373/1‑149 (MQ=255)                                                                              cGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTcgc                                                                       <  6:169339/149‑1 (MQ=255)                                                                                ccTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTcgcg                                                                      >  2:215100/1‑148 (MQ=255)                                                                                    ccAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCAGCCGTGTCGCGCAACTGCTCGCGaaa                                                                   <  3:61406/147‑1 (MQ=255)                                                                                    ccAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAaca                                                                 >  1:352235/1‑149 (MQ=255)                                                                                        gCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTc                                                             >  6:206058/1‑149 (MQ=255)                                                                                          cTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTc                                                             >  2:147826/1‑147 (MQ=255)                                                                                          cTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCa                                                           <  5:364601/149‑1 (MQ=255)                                                                                           tGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAg                                                          >  2:79607/1‑149 (MQ=255)                                                                                              tCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTc                                                       >  7:2077/1‑149 (MQ=255)                                                                                                cGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTc                                                       <  1:147826/147‑1 (MQ=255)                                                                                                cGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGc                                                     >  2:61077/1‑149 (MQ=255)                                                                                                cGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACGGCTCGCGAAACACTTCCAGGTCGc                                                     >  4:279445/1‑149 (MQ=255)                                                                                                 ggTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCg                                                      <  8:223238/147‑1 (MQ=255)                                                                                                 ggTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCg                                                      >  7:223238/1‑147 (MQ=255)                                                                                                  gTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTg                                                   <  2:352235/149‑1 (MQ=255)                                                                                                       aGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCATGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACtt                                                              <  7:286929/133‑1 (MQ=255)                                                                                                       aGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCacac                                              >  6:223300/1‑149 (MQ=255)                                                                                                       aGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTAGCGCAACTGCTCGCGAAACACtt                                                              >  8:286929/1‑133 (MQ=255)                                                                                                           ccTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACaccccc                                          >  5:363743/1‑149 (MQ=255)                                                                                                               gcgcTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGt                                       <  5:206058/148‑1 (MQ=255)                                                                                                               gcgcTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGt                                       <  1:61077/148‑1 (MQ=255)                                                                                                                      ttAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCCAGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATAGGc                               >  5:89645/1‑149 (MQ=255)                                                                                                                          cggcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATc                           <  7:157247/149‑1 (MQ=255)                                                                                                                            gcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTg                                                   >  2:192324/1‑123 (MQ=255)                                                                                                                            gcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCAc                         >  5:324777/1‑149 (MQ=255)                                                                                                                            gcgATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCAc                         >  1:237577/1‑149 (MQ=255)                                                                                                                              gATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCg                       >  6:380191/1‑149 (MQ=255)                                                                                                                                ttGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGcc                     <  8:2077/149‑1 (MQ=255)                                                                                                                                   tCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCctct                  <  5:223300/149‑1 (MQ=255)                                                                                                                                   tCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCctct                  <  1:215100/149‑1 (MQ=255)                                                                                                                                     gATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCtt                 <  1:171709/148‑1 (MQ=255)                                                                                                                                     gATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCCGTACATCGGCCATCACCGCCTCTTg                >  1:168867/1‑149 (MQ=255)                                                                                                                                      aTGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTg                >  2:185393/1‑148 (MQ=255)                                                                                                                                          aCTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCtt                 >  6:21873/1‑143 (MQ=255)                                                                                                                                          aCTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCtt                 <  5:21873/143‑1 (MQ=255)                                                                                                                                           cTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGTCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGcaccac          <  5:380191/149‑1 (MQ=255)                                                                                                                                            tGAGCGCCGTACAGCGCTCGATCGCCTCAGCGCCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCCGCAcaccac                                           >  1:40626/1‑116 (MQ=255)                                                                                                                                            tGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGcaccac          >  7:97965/1‑148 (MQ=255)                                                                                                                                            tGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCCGCAcaccac                                           <  2:40626/115‑1 (MQ=255)                                                                                                                                             gAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGcaccac          >  2:64598/1‑147 (MQ=255)                                                                                                                                                cgcCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCCTCGGCAACCGGCAGGTAGCCACCCGTGTCTCGCAACTGCTCGCCAAACACTACCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGcc     >  2:283322/1‑149 (MQ=255)                                                                                                                                                   cGTACCGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCGCAGCACCGGTACATCGGCCAGCCCCGCCTCTTGGACCACCCGccc    >  4:130825/1‑147 (MQ=255)                                                                                                                                                   cGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGccc    >  3:65631/1‑147 (MQ=255)                                                                                                                                                   cGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCa  <  7:308952/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CCCTTGCCGCCTTCGAGGCAGCCTTGCTGACACGCTGAAGGTCAGCTGCGCAACAGGCTCTGCGAGCCGCTGCAGAACGCCTGCCAGGCCTGCTCCGGTGCCAGTGCCTGGCGCTGCTTAGCGGCGATTGTCGATGCACTGAGCGCCGTACAGCGCTCGATCGCCTCAGC‑CCAGGCGCTTTCGTCGGCAACCGGCAGGTAGCCACCCGTGTCGCGCAACTGCTCGCGAAACACTTCCAGGTCGCTGCACACCACCGGTACATCGGCCATCACCGCCTCTTGCACCACCAGCCCCA  >  NC_002947/5619320‑5619614

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