Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,875,947	A→G	6.1%	S379G (AGC→GGC)	PP_4284 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,875,947	0	A	G	6.1%	188.6 / 2.6	65	S379G (AGC→GGC)	PP_4284	transporter
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base G (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAG  >  NC_002947/4875800‑4876093                                                                                                                                                    |                                                                                                                                                   ggtGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGgcag                                                                                                                                                   <  4:197620/147‑1 (MQ=255)   tGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGATGATGATGGgcagca                                                                                                                                                 >  6:391803/1‑149 (MQ=255)      cGCAGGTGGCCGAACTGGCCTGAGGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATgg                                                                                                                                              >  8:235550/1‑149 (MQ=255)       gCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCTCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGgc                                                                                                                                             >  6:235358/1‑149 (MQ=255)           gTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGca                                                                                                                                          >  2:387521/1‑148 (MQ=255)                cGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCa                                                                                                                                    >  5:341566/1‑149 (MQ=255)                 gAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCAt                                                                                                                                   >  3:201496/1‑149 (MQ=255)                  aaCTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCAtt                                                                                                                                  >  8:361117/1‑149 (MQ=255)                         cTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGa                                                                                                                           <  5:68282/149‑1 (MQ=255)                         cTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGa                                                                                                                           <  5:62224/149‑1 (MQ=255)                            aCGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAg                                                                                                                          <  1:387521/147‑1 (MQ=255)                             cGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGTGCGCCATGCTGCCCATcatca                                                                                                                                  >  2:123690/1‑134 (MQ=255)                             cGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGTGCGCCATGCTGCCCATcatca                                                                                                                                  <  1:123690/138‑5 (MQ=255)                              ggCTGTGGTGGTAGGCCTGGTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATTCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGcc                                                                                                                       <  4:106534/148‑1 (MQ=255)                              ggCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCa                                                                                                                      >  4:280085/1‑149 (MQ=255)                               gCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCAc                                                                                                                     <  6:341566/149‑1 (MQ=255)                                     gtggtAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGTGCGCCATGCTGCCCATc                                                                                                                                      >  4:242711/1‑126 (MQ=255)                                     gtggtAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGTGCGCCATGCTGCCCATc                                                                                                                                      <  3:242711/126‑1 (MQ=255)                                        gtAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCAt                                                                                                            >  1:285828/1‑149 (MQ=255)                                             ccTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCgg                                                                                                       >  3:287088/1‑149 (MQ=255)                                                 ctgTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGAt                                                                                                                             >  1:269151/1‑123 (MQ=255)                                                 ctgTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGAt                                                                                                                             <  2:269151/123‑1 (MQ=255)                                                   gTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCg                                                                                                 >  6:137065/1‑149 (MQ=255)                                                    ttCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATc                                                                                                  <  2:19017/147‑1 (MQ=255)                                                      cATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTca                                                                                              >  3:111775/1‑149 (MQ=255)                                                       aTCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTca                                                                                              <  7:361117/148‑1 (MQ=255)                                                       aTCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATCATTCCCGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTca                                                                                              >  6:277880/1‑148 (MQ=255)                                                        tCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCAtt                                                                                                                                  <  2:325019/111‑1 (MQ=255)                                                        tCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCTTTCCGGCGATCGTca                                                                                              >  8:358435/1‑147 (MQ=255)                                                        tCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACAACGAtt                                                                                                                                  >  1:325019/1‑111 (MQ=255)                                                           ccgcGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGt                                                                                          >  8:267280/1‑148 (MQ=255)                                                           ccgcGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGt                                                                                          >  8:242136/1‑148 (MQ=255)                                                            cgcgATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGa                                                                                        <  7:267280/149‑1 (MQ=255)                                                            cgcgATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGa                                                                                        >  2:193450/1‑149 (MQ=255)                                                             gcgATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGAt                                                                                       <  4:287088/149‑1 (MQ=255)                                                                     cttcGCTCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTTCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCatgatg                                                                                <  1:382373/148‑1 (MQ=255)                                                                     cttcGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCatgatg                                                                                <  4:201496/148‑1 (MQ=255)                                                                     cttcGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCatgatg                                                                                <  7:235550/148‑1 (MQ=255)                                                                      ttcGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGcc                                                                              >  6:148238/1‑149 (MQ=255)                                                                        cgcgCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCg                                                                             <  5:391803/148‑1 (MQ=255)                                                                        cgcgCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGc                                                                            >  4:1774/1‑149 (MQ=255)                                                                         gcgcCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCt                                                                           >  3:361370/1‑149 (MQ=255)                                                                          cgccgcTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGACACAGTGATCATGATGCCGCTg                                                                          <  7:242136/149‑1 (MQ=255)                                                                                    cGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCgg                                                                >  4:140186/1‑149 (MQ=255)                                                                                     gggATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGt                                                               >  1:333102/1‑149 (MQ=255)                                                                                      ggATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGAGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATTCCGCTGACGTTCTCGGTc                                                              <  5:235358/149‑1 (MQ=255)                                                                                           atTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGa                                                         >  6:186781/1‑149 (MQ=255)                                                                                             tCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATg                                                       >  5:300460/1‑149 (MQ=255)                                                                                              ccGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATg                                                       <  3:1774/148‑1 (MQ=255)                                                                                                   cTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTg                                                 <  5:148238/149‑1 (MQ=255)                                                                                                        cGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTg                                            <  3:280085/149‑1 (MQ=255)                                                                                                               agGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGTGATGAGGCCACCGACAGCATTCCGGAGATCGTCACAGTTATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGTCTtca                                     <  5:277880/149‑1 (MQ=255)                                                                                                                ggTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTtcat                                    <  3:140186/149‑1 (MQ=255)                                                                                                                ggTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCGGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTtcat                                    <  2:285828/149‑1 (MQ=255)                                                                                                                    cACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGc                                >  7:17141/1‑149 (MQ=255)                                                                                                                         aTCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATAATGATGCCGCTGACGTTATCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTAtgt                           >  4:309449/1‑149 (MQ=255)                                                                                                                            tACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTgg                         >  5:79772/1‑148 (MQ=255)                                                                                                                              cGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGc                        <  5:137065/147‑1 (MQ=255)                                                                                                                                        cTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCgg             <  4:111775/148‑1 (MQ=255)                                                                                                                                        cTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCgg             <  2:333102/148‑1 (MQ=255)                                                                                                                                        cTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCgg             <  8:57176/148‑1 (MQ=255)                                                                                                                                         tgatgaTGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTtc                                      <  6:261284/122‑1 (MQ=255)                                                                                                                                         tgatgaTGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTtc                                      >  5:261284/1‑122 (MQ=255)                                                                                                                                          gatgatGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCAGCGGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTa          >  7:249740/1‑149 (MQ=255)                                                                                                                                              atgGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCg                                                             <  4:378351/94‑1 (MQ=255)                                                                                                                                              atgGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCg                                                             >  3:378351/1‑94 (MQ=255)                                                                                                                                                gggCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGc     >  1:277565/1‑148 (MQ=255)                                                                                                                                                  gcagcaTGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAg  >  7:190619/1‑149 (MQ=255)                                                                                                                                                  gcagcaTGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAg  >  4:341665/1‑149 (MQ=255)                                                                                                                                                   cagcaTGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAg  >  4:206307/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GTTGCCGCAGGTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAG  >  NC_002947/4875800‑4876093

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