Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,403,828	G→C	5.4%	T70S (ACC→AGC)	PP_3883 ←	holin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,403,828	0	G	C	5.4%	298.5 / 6.3	92	T70S (ACC→AGC)	PP_3883	holin
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base C (3/2);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

CCATTTTCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATC  >  NC_002947/4403681‑4403969                                                                                                                                                    |                                                                                                                                               ccATTTTCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGt                                                                                                                                               >  2:306241/1‑149 (MQ=255)  cATTTTCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTg                                                                                                                                              <  7:184898/149‑1 (MQ=255)     tttCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGCCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGca                                                                                                                                            >  8:240892/1‑148 (MQ=255)      ttCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGca                                                                                                                                            >  7:184834/1‑147 (MQ=255)       tCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGcaca                                                                                                                                          <  3:325224/148‑1 (MQ=255)        cATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAg                                                                                                                                         <  1:64547/148‑1 (MQ=255)         aTGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCg                                                                                                                                       <  6:181756/149‑1 (MQ=255)            gggCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGAc                                                                                                                                     <  1:292785/148‑1 (MQ=255)             ggCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTg                                                                                                                                   >  4:165489/1‑149 (MQ=255)              gCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGc                                                                                                                                  <  6:294689/149‑1 (MQ=255)               ccccGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGCTGTGCACCAACTGcc                                                                                                                                 <  4:188593/149‑1 (MQ=255)                cccGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGCTGTGCACCAACTGCCa                                                                                                                                <  4:226186/149‑1 (MQ=255)                 ccGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAg                                                                                                                               >  5:74950/1‑149 (MQ=255)                          ccGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTcc                                                                                                                        >  3:228005/1‑147 (MQ=255)                             gAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGcc                                                                                                                   >  2:335319/1‑149 (MQ=255)                               ggCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCa                                                                                                                 <  4:183670/149‑1 (MQ=255)                                    tCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCa                                                                                                                 >  1:351893/1‑144 (MQ=255)                                    tCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCa                                                                                                                 <  2:351893/144‑1 (MQ=255)                                     cAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGaa                                                                                                              >  2:71688/1‑146 (MQ=255)                                      aggaggGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACaga                                                                                                          <  4:372028/149‑1 (MQ=255)                                        gaggGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACagag                                                                                                         <  4:408091/148‑1 (MQ=255)                                        cggggCCTATTCCAGTGGTGCCGGTGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACtga                                                                                                          <  2:79879/145‑1 (MQ=255)                                             ccTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGt                                                                                                   >  6:351364/1‑149 (MQ=255)                                              cTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGTGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTc                                                                                                  >  5:227306/1‑149 (MQ=255)                                              cTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGACCAGAGTGATGTc                                                                                                  >  2:424329/1‑149 (MQ=255)                                               tATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCAGACTTCTCCAGCCGCTGCTTGAAGGAGTTGTGGCAGTAGGTGCACAGCGACAGCCACTTGGTCCGAGCCCAGAACAGAGTCATGTcg                                                                                                 >  4:330004/1‑149 (MQ=255)                                                   ccAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTcg‑ccg                                                                                             <  1:335319/148‑1 (MQ=255)                                                    cAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAAc                                                                                                             <  6:73666/132‑1 (MQ=255)                                                    cAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAAc                                                                                                             >  5:73666/1‑132 (MQ=255)                                                     aGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCgg                                                                                          >  5:429439/1‑149 (MQ=255)                                                      gtggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGt                                                                                         >  8:21912/1‑149 (MQ=255)                                                      gtggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGt                                                                                         >  2:306556/1‑149 (MQ=255)                                                       tggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTg                                                                                        <  5:351364/149‑1 (MQ=255)                                                         gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                      <  6:368851/149‑1 (MQ=255)                                                         gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                      >  4:323369/1‑149 (MQ=255)                                                         gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                      >  1:362020/1‑149 (MQ=255)                                                            ccGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGc                                                                                    <  8:324520/148‑1 (MQ=255)                                                            ccGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCt                                                                                   <  3:165489/149‑1 (MQ=255)                                                              gggggTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTg                                                                                  >  8:91443/1‑148 (MQ=255)                                                              gggggTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGa                                                                                 <  6:272818/149‑1 (MQ=255)                                                                gggTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAAt                                                                               <  3:323369/149‑1 (MQ=255)                                                                  gTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAAtgt                                                                             >  7:381669/1‑149 (MQ=255)                                                                    ttaGCGGCAGGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGa                                                                           >  7:376544/3‑149 (MQ=255)                                                                     cAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGAt                                                                          >  3:132747/1‑149 (MQ=255)                                                                      aGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATc                                                                         >  3:361617/1‑149 (MQ=255)                                                                      aGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATc                                                                         >  6:212824/1‑149 (MQ=255)                                                                          gCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATcgac                                                                      >  1:91793/1‑148 (MQ=255)                                                                              gCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAg                                                                 >  3:214931/1‑149 (MQ=255)                                                                               ccGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACa                                                                  <  4:132747/147‑1 (MQ=255)                                                                                 gCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGa                                                                           <  1:93195/136‑1 (MQ=255)                                                                                 gCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGa                                                                           >  2:93195/1‑136 (MQ=255)                                                                                     aCGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCGTCCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCg                                                            <  1:71688/148‑1 (MQ=255)                                                                                      cGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGc                                                         <  2:362020/149‑1 (MQ=255)                                                                                       gTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCt                                                        >  5:329029/1‑149 (MQ=255)                                                                                       gTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCt                                                        >  2:116141/1‑149 (MQ=255)                                                                                        tCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGCATGGGAGCGACGCCAGTGGCGGGCtt                                                       >  6:333481/1‑149 (MQ=255)                                                                                         cGCAACCAGCCTCGCGCCCGGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGCGATCGACGACAGTGGCGGGCTTg                                                      <  8:184834/149‑1 (MQ=255)                                                                                          gCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGc                                                     >  5:33829/1‑149 (MQ=255)                                                                                           cAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGTGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGATAGTGGCGGGCTTGCt                                                    <  6:227306/149‑1 (MQ=255)                                                                                           cAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCt                                                    >  1:431998/1‑149 (MQ=255)                                                                                                 gCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCtt                                                       <  2:279678/140‑1 (MQ=255)                                                                                                 gCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCtt                                                       >  1:279678/1‑140 (MQ=255)                                                                                                  ccTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCt                                             <  5:205708/149‑1 (MQ=255)                                                                                                      gcgcCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACt                                         <  7:91443/149‑1 (MQ=255)                                                                                                         cccTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGa                                      <  4:361617/149‑1 (MQ=255)                                                                                                          ccTGACTTCTCCAGGCGCTGCTTGTCGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGGCCGGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAg                                     <  6:74950/149‑1 (MQ=255)                                                                                                          ccTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAg                                     <  7:21912/149‑1 (MQ=255)                                                                                                          ccTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAg                                     <  2:431998/149‑1 (MQ=255)                                                                                                           cTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAgg                                    <  7:240892/149‑1 (MQ=255)                                                                                                              aCTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGc                                   <  1:116141/147‑1 (MQ=255)                                                                                                                ttCTCCAGGCGCTGCTTCTAGGAGTTGTGGCAGTTGGAGCACATCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCT‑CCGCTGTGCGGCTGAATGTGAACGACGACAGTGGCGGGCATCCTCAGGCCTCACAGGAGGCAGc                                >  8:398608/1‑148 (MQ=255)                                                                                                                  ctcCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGaa                             >  6:214732/1‑149 (MQ=255)                                                                                                                   tcCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGaa                             >  4:228860/1‑148 (MQ=255)                                                                                                                     gaGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACTGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCTGGCTTGCTCAGGCCTGACTGGAGGCAGCGAaca                           <  7:398608/147‑1 (MQ=255)                                                                                                                     cAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAaca                           >  3:108827/1‑148 (MQ=255)                                                                                                                          gctgctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGc                     <  2:91793/149‑1 (MQ=255)                                                                                                                           ctgctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCg                    >  5:202032/1‑149 (MQ=255)                                                                                                                              ctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCgg                   <  3:228860/147‑1 (MQ=255)                                                                                                                              ctTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGt                 <  5:212824/149‑1 (MQ=255)                                                                                                                               ttGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGtt                >  6:224689/1‑149 (MQ=255)                                                                                                                                      ccaactgGCAGTCGCTGTGCACCAACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGCGAGGGTTCTTgcg          >  7:51417/6‑148 (MQ=255)                                                                                                                                      aGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTgcgc         <  7:64599/149‑1 (MQ=255)                                                                                                                                      aGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGc                                               <  1:16598/111‑1 (MQ=255)                                                                                                                                      aGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGc                                               >  2:16598/1‑111 (MQ=255)                                                                                                                                       gTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTgcgc         >  8:363955/1‑148 (MQ=255)                                                                                                                                         actgGCAGTCGCTGTGCACCAACTGCCAGTTGGTCTGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGTCGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGCACACAGAGCGGGTTCTTGCGCAAc      >  4:211219/3‑149 (MQ=255)                                                                                                                                         actgGCAGTCGCTGTGCACCAACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAAc      >  6:126426/3‑149 (MQ=255)                                                                                                                                          gtgGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAAc      >  6:80168/1‑148 (MQ=255)                                                                                                                                           tgGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAAcc     <  6:202032/148‑1 (MQ=255)                                                                                                                                             gCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATc  >  5:385307/1‑149 (MQ=255)                                                                                                                                              cAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCGTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATc  <  4:371250/148‑1 (MQ=255)                                                                                                                                              cAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATc  >  3:371250/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                               CCATTTTCATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATC  >  NC_002947/4403681‑4403969

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