Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,475,587	A→T	6.1%	intergenic (‑52/+12)	PP_3964 ← / ← PP_3965	transposase/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,475,587	0	A	T	6.1%	189.0 / 2.9	65	intergenic (‑52/+12)	PP_3964/PP_3965	transposase/transposase
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

CTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGT  >  NC_002947/4475446‑4475732                                                                                                                                              |                                                                                                                                                   cTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCATTTGGCCg                                                                                                                                             <  7:73456/149‑1 (MQ=255)    gACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCATTTGGCCggtc                                                                                                                                          <  6:23488/149‑3 (MQ=255)      ccTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGcc                                                                                                                                              <  7:110517/143‑1 (MQ=255)      ccTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGcc                                                                                                                                              >  8:110517/1‑143 (MQ=255)       cTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCg                                                                                                                                             >  7:240745/1‑143 (MQ=255)       cTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCg                                                                                                                                             <  8:240745/143‑1 (MQ=255)            attgcaattgcaACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTa                                                                                                                                     <  8:314038/145‑1 (MQ=255)            attgcaattgcaACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTa                                                                                                                                     >  7:314038/1‑145 (MQ=255)                 aattgcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGc                                                                                                                             >  4:8769/1‑148 (MQ=255)                     gcaaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAAt                                                                                                                         <  4:101544/148‑1 (MQ=255)                       aaCGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGAGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATca                                                                                                                       >  5:277014/1‑148 (MQ=255)                          gCCTGTGCCGCCAGGTCACGAAGCTGTTCCGGGTTCATCTTGTCGAGGTTGTGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGTGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATTCTAATCACAcc                                                                                                                   <  6:179263/149‑1 (MQ=17)                          gCCTGTGCCGCCAGGGCACGCAGCTTTACCGGGGTCATCAGGTCGAGGTTGGGAGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCCGGGCGACCACAAGTAGACCCGCCAAATCGCCGGCGT‑ATACAGCATGATAATCACAcc                                                                                                                   >  5:179263/1‑149 (MQ=17)                            cTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTTATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGc                                                                                                                 >  8:309100/1‑149 (MQ=255)                             tgtgCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCacac                                                                                                                    <  2:32315/145‑1 (MQ=255)                             tgtgCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCacac                                                                                                                    >  1:32315/1‑145 (MQ=255)                                gccgccAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTg                                                                                                              >  7:106463/1‑148 (MQ=255)                                  cgccAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCAc                                                                                                           >  2:238394/1‑149 (MQ=255)                                      aaGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGAc                                                                                                       >  5:176272/1‑149 (MQ=255)                                           aCGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCt                                                                                                  >  1:80430/1‑149 (MQ=255)                                            cGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTg                                                                                                 <  8:106463/149‑1 (MQ=255)                                             gAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGc                                                                                                <  3:8769/149‑1 (MQ=255)                                                 cTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACg                                                                                            >  4:170189/1‑149 (MQ=255)                                                   gTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACgg                                                                                           >  6:176931/1‑148 (MQ=255)                                                    ttCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCa                                                                                         <  1:238394/149‑1 (MQ=255)                                                              aTCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCAc                                                                               >  5:316278/1‑149 (MQ=255)                                                              aTCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCAc                                                                               >  2:324886/1‑149 (MQ=255)                                                              aTCTGGTCGAGGTTGGGCGAGAAAGTAATGCCGCCGATTGTGCAAGAGCAGGCGCGGGGCGAAGACAATAAGACCGGCCAAATGGCCGGCGT‑TTAAAGAATGCTAATCAAACCGCCTGCACCGACCCGCAGCAACGGCAGGCCAAGCa                                                                                >  7:220865/1‑148 (MQ=255)                                                                   gTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACg                                                                          >  6:185791/1‑149 (MQ=255)                                                                         gTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGa                                                                    <  2:80430/149‑1 (MQ=255)                                                                              gCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgttgt                                                               <  5:185791/149‑1 (MQ=255)                                                                              gCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgttgt                                                               >  1:159207/1‑149 (MQ=255)                                                                                 agaAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACTCCAGGCgttggtgc                                                            <  5:370399/149‑5 (MQ=255)                                                                                  gaAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgtt                                                                 <  5:103358/143‑1 (MQ=255)                                                                                  gaAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAgtt                                                                 >  6:103358/1‑143 (MQ=255)                                                                                  gaAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGc                                                           <  3:201256/149‑1 (MQ=255)                                                                                   aaaGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCCCCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCt                                                          >  5:120884/1‑149 (MQ=255)                                                                                    aaGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTg                                                         <  6:277014/149‑1 (MQ=255)                                                                                     aGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGt                                                        >  1:105864/1‑149 (MQ=255)                                                                                     aGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGt                                                        <  6:316278/149‑1 (MQ=255)                                                                                     aGTCATCCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATAACACCGCCTCCACCGCCGCCCTGCCACGGCAGCCCCAGCACCAACGCCTGGACTTGTTCGCAGt                                                        >  6:268365/1‑149 (MQ=18)                                                                                      ttcttgtcgccgTTTGTGCCAGTTCAGGCGCGGGGCGACGACAATTATACCGGCCAAATTGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCAGCTACGCGCTGCCACGGCAGGCCCAGCACAAACGCCTGGAGTTGTTCGCTGTc                                                       <  8:220865/142‑1 (MQ=18)                                                                                      gTCATGCCGCCGTTTGTGCGAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCTAATGGCCGGCGT‑TTACAGCATGCTAATCACAGCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTc                                                       <  5:268365/149‑1 (MQ=255)                                                                                       tCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTc                                                       >  8:194682/1‑148 (MQ=255)                                                                                         aTGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACgg                                                                                           >  4:104857/1‑110 (MQ=255)                                                                                         aTGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACgg                                                                                           <  3:104857/110‑1 (MQ=255)                                                                                            ccgccgATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCCCCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTc                                                 >  6:353603/1‑149 (MQ=255)                                                                                               ccgATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGAc                                              >  1:220753/1‑149 (MQ=255)                                                                                                cgATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGAcc                                             >  4:63348/1‑149 (MQ=255)                                                                                                      tgCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAg                                       >  3:164740/1‑149 (MQ=255)                                                                                                           gagCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTg                                  <  7:309100/149‑1 (MQ=255)                                                                                                            ccataacgcgGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTg                                  <  8:221600/142‑1 (MQ=255)                                                                                                                       ggCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCggg                      <  5:176931/149‑1 (MQ=255)                                                                                                                       ggCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCCGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCggg                      <  1:324886/149‑1 (MQ=255)                                                                                                                                  gTAGACCGGCCAAATGGCCGGCATGAAACAGCATGATAAACACACAGCCAGCACCGACGAGCAGCCACGGCAGGCCCAGCACCACCGCGTGTAGTTGTGAGCTGGGGAGTGCGACGGCGCAGCGGGGGCGAGGGCCGGGCGGGTGGagg            >  7:224854/1‑147 (MQ=255)                                                                                                                                   aacgccGGCCATTTGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACtt          >  2:224244/5‑149 (MQ=255)                                                                                                                                     gACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTg         <  8:224854/148‑1 (MQ=255)                                                                                                                                     gACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGc        <  2:220753/149‑1 (MQ=255)                                                                                                                                      gccGGCCATTTGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGcc       >  6:42105/2‑149 (MQ=255)                                                                                                                                      aCCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGCCCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGATCTTGcg       >  8:239516/1‑148 (MQ=255)                                                                                                                                       ccGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCt      <  2:105864/149‑1 (MQ=255)                                                                                                                                         ggCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTg    <  8:185462/149‑1 (MQ=255)                                                                                                                                           ccAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGt  >  4:146656/1‑149 (MQ=255)                                                                                                                                            cAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCt      >  6:365050/1‑144 (MQ=255)                                                                                                                                            cAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCt      <  5:365050/144‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                   CTCGACCTGGGATTGCAATTGCAACGCCTGTGCCGCCAAGGCACGAAGCTGTTCCGGGGTCATCTGGTCGAGGTTGGGCGAGAAAGTCATGCCGCCGATTGTGCCAGAGCAGGCGCGGGGCGACGACAAGTAGACCGGCCAAATGGCCGGCGT‑TTACAGCATGCTAATCACACCGCCTGCACCGACGCGCTGCCACGGCAGGCCCAGCACCAACGCCTGGAGTTGTTCGCTGTCGAGTTCGACCTCGCAGCCGTGGCGAATGCCGGGCCAGTGGAACTTGCCTTGGT  >  NC_002947/4475446‑4475732

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑