Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,403,835	G→C	5.3%	L68V (CTG→GTG)	PP_3883 ←	holin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,403,835	0	G	C	5.3%	307.6 / 6.3	94	L68V (CTG→GTG)	PP_3883	holin
Reads supporting (aligned to +/- strand):  ref base G (46/42);  new base C (3/2);  total (50/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 = 3.84e-01

CATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGCCTG  >  NC_002947/4403688‑4403982                                                                                                                                                    |                                                                                                                                                     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                                                                                                                          >  5:73666/1‑132 (MQ=255)                                             cAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAAc                                                                                                                          <  6:73666/132‑1 (MQ=255)                                              aGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCgg                                                                                                       >  5:429439/1‑149 (MQ=255)                                               gtggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGt                                                                                                      >  2:306556/1‑149 (MQ=255)                                               gtggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGt                                                                                                      >  8:21912/1‑149 (MQ=255)                                                tggtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTg                                                                                                     <  5:351364/149‑1 (MQ=255)                                                  gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                                   >  1:362020/1‑149 (MQ=255)                                                  gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                                   >  4:323369/1‑149 (MQ=255)                                                  gtgCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCg                                                                                                   <  6:368851/149‑1 (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                                                                     >  2:116141/1‑149 (MQ=255)                                                                                gTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCt                                                                     >  5:329029/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)                                                                                                                                              tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCa          <  2:186949/146‑1 (MQ=255)                                                                                                                                              tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCa          >  1:186949/1‑146 (MQ=255)                                                                                                                                              tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCAc         >  5:185377/1‑147 (MQ=255)                                                                                                                                              tGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGa       <  6:329029/149‑1 (MQ=255)                                                                                                                                                cacaGCGACTGCCAGTTGGTCCGAGTCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGc     >  3:10043/1‑149 (MQ=255)                                                                                                                                                cacaGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGc     >  4:106893/1‑149 (MQ=255)                                                                                                                                                cacaGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGc     <  1:306556/149‑1 (MQ=255)                                                                                                                                                   aGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGAGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGCctg  >  5:220605/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGCCTG  >  NC_002947/4403688‑4403982

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