Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,403,837	G→A	5.2%	S67L (TCG→TTG)	PP_3883 ←	holin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,403,837	0	G	A	5.2%	294.2 / 3.1	95	S67L (TCG→TTG)	PP_3883	holin
Reads supporting (aligned to +/- strand):  ref base G (48/42);  new base A (3/2);  total (51/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 = 2.80e-01

ATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGCCTGCT  >  NC_002947/4403689‑4403984                                                                                                                                                     |                                                                                                                                                     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                                                                      >  1:279678/1‑140 (MQ=255)                                                                                         gCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCtt                                                                      <  2:279678/140‑1 (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                                                    <  2:431998/149‑1 (MQ=255)                                                                                                  ccTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAg                                                    <  7:21912/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                                                              >  2:16598/1‑111 (MQ=255)                                                                                                                              aGTTGTGGCAGTCGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGc                                                              <  1:16598/111‑1 (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)                                                                                                                                                    cGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGCctgct  >  2:33690/1‑149 (MQ=255)                                                                                                                                                    cGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGCctgct  <  3:330004/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ATGGGGCCCCGCAGCGCCCGGAGGCTTTCAGGAGGGCCTATTCCAGTGGTGCCGGGGGTCCAGCGGCCTGCCGCTAACGTCGCAACCAGCCTCGCGCCCTGACTTCTCCAGGCGCTGCTTGTAGGAGTTGTGGCAGTTGGTGCACAGCGACTGCCAGTTGGTCCGAGCCCAGAACAGAGTCATGTCG‑CCGCGGTGCGGCTGAATGTGATCGACGACAGTGGCGGGCTTGCTCAGGCCTGACTGGAGGCAGCGAACACAGAGCGGGTTCTTGCGCAACCATCCTTCACGAGCCTGCT  >  NC_002947/4403689‑4403984

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