Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,875,957	A→C	6.0%	H382P (CAC→CCC)	PP_4284 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,875,957	0	A	C	6.0%	210.9 / 3.6	67	H382P (CAC→CCC)	PP_4284	transporter
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base C (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAGTACAAGGA  >  NC_002947/4875810‑4876101                                                                                                                                                    |                                                                                                                                                 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:62224/149‑1 (MQ=255)               cTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGa                                                                                                                                   <  5:68282/149‑1 (MQ=255)                  aCGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAg                                                                                                                                  <  1:387521/147‑1 (MQ=255)                   cGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGTGCGCCATGCTGCCCATcatca                                                                                                                                          <  1:123690/138‑5 (MQ=255)                   cGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGTGCGCCATGCTGCCCATcatca                                                                                                                                          >  2:123690/1‑134 (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                                                                                                >  2:193450/1‑149 (MQ=255)                                                  cgcgATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGa                                                                                                <  7:267280/149‑1 (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                                                                     >  3:378351/1‑94 (MQ=255)                                                                                                                                    atgGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCg                                                                     <  4:378351/94‑1 (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)                                                                                                                                              tGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAGTACAAg    <  5:186781/149‑1 (MQ=255)                                                                                                                                                 cgcACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAGTACAAGGa  >  4:281851/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GTGGCCGAACTGGCCTGACGGCTGTGGTGGTAGGCCTGCTGTTCATCGCCGCGATGTTCTTCGCGCCGCTGGCCGGGATGATTCCGGCCTATGCGACAGCAGGTGCACTGATCTACGTTGCGATGCTGATGATGGGCAGCATGGCGCACATCCATTGGGATGAGGCCACCGACAGCATTCCGGCGATCGTCACAGTGATCATGATGCCGCTGACGTTCTCGGTCGCCGATGGTATTGCCCTGGGCTTCATCAGCTATGTGGCATTGAAGGCGGGTACCGGCAAGTACAAGGA  >  NC_002947/4875810‑4876101

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