Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,650,676	G→T	5.6%	H173N (CAT→AAT) ‡	PP_4959 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,650,676	0	G	T	5.6%	294.3 / 5.9	90	H173N (CAT→AAT) ‡	PP_4959	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (50/35);  new base T (3/2);  total (53/37)
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.02e-01

CCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGA  >  NC_002947/5650529‑5650822                                                                                                                                                    |                                                                                                                                                   ccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCTGCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                   <  4:155509/149‑1 (MQ=255)  cTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATggg                                                                                                                                                  >  3:139031/1‑149 (MQ=255)     gTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCAt                                                                                                                                               >  4:288022/1‑149 (MQ=255)        ggCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCat                                                                                                                                            <  2:77014/149‑2 (MQ=255)         gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCatc                                                                                                                                           <  6:355020/149‑3 (MQ=255)         gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                   >  1:328429/1‑141 (MQ=255)         gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                   <  2:328429/141‑1 (MQ=255)         gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAac                                                                                                                                           >  3:198902/1‑149 (MQ=255)           cTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAGACGGATGATGCCATGGGCATCAacca                                                                                                                                         >  3:1372/1‑149 (MQ=255)              gTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAaccac                                                                                                                                        <  4:307749/147‑1 (MQ=255)                  aaGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAAt                                                                                                                                  >  7:200615/1‑149 (MQ=255)                   aGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATg                                                                                                                                 <  4:198902/149‑1 (MQ=255)                    gTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGc                                                                                                                                >  6:238192/1‑149 (MQ=255)                     tCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATg                                                                                                                                 <  1:292864/147‑1 (MQ=255)                     tCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATg                                                                                                                                 >  2:292864/1‑147 (MQ=255)                        gACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTc                                                                                                                            <  7:145136/149‑1 (MQ=255)                                 cacaGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACAGGTGCCGCCAGCAAACGCGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCAtt                                                                                                                   >  7:116850/1‑149 (MQ=255)                                  acaGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCAttt                                                                                                                  >  4:135601/1‑149 (MQ=255)                                    aGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTc                                                                                                                >  3:316562/1‑149 (MQ=255)                                     ggCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTcc                                                                                                               <  3:818/149‑1 (MQ=255)                                      gCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTc                                                                                                                >  5:13541/1‑147 (MQ=255)                                      gCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                              <  7:47842/149‑1 (MQ=255)                                           tACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAacc                                                                                                                                          <  5:7336/116‑1 (MQ=255)                                           tACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAacc                                                                                                                                          >  6:7336/1‑116 (MQ=255)                                           tACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATg                                                                                                         >  5:273923/1‑149 (MQ=255)                                              tGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                   <  5:214459/104‑1 (MQ=255)                                              tGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATgg                                                                                                                                                   >  6:214459/1‑104 (MQ=255)                                               ggTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAATGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACg                                                                                                     >  7:161511/1‑149 (MQ=255)                                               ggTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACg                                                                                                     >  2:42978/1‑149 (MQ=255)                                                gTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCTTTTTCCAGGATGGACg                                                                                                     >  5:107956/1‑148 (MQ=255)                                                 tCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGTAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCa                                                                                                   >  6:350497/1‑149 (MQ=255)                                                       gCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAg                                                                                             >  6:186621/1‑149 (MQ=255)                                                        cacaACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGa                                                                                               <  5:201164/146‑1 (MQ=255)                                                        cacaACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGa                                                                                               >  6:201164/1‑146 (MQ=255)                                                        cacaACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAg                                                                                             <  3:235048/148‑1 (MQ=255)                                                        cacaACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAg                                                                                             <  6:134181/148‑1 (MQ=255)                                                         acaACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGc                                                                                           >  5:296536/1‑149 (MQ=255)                                                          caACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGcc                                                                                          <  5:34312/149‑1 (MQ=255)                                                           aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTcg                                                                                                                           >  3:92488/1‑115 (MQ=255)                                                           aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTcg                                                                                                                           <  4:92488/115‑1 (MQ=255)                                                           aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGccc                                                                                         >  6:181472/1‑149 (MQ=255)                                                               tCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAATGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGc                                                                                     <  8:161511/149‑1 (MQ=255)                                                               tCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGc                                                                                     <  6:273923/149‑1 (MQ=255)                                                                     aaTTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGtctc                                                                               <  5:238192/149‑1 (MQ=255)                                                                       ttGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGtctc                                                                               >  7:63892/1‑147 (MQ=255)                                                                       ttGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGtctc                                                                               >  7:11899/1‑147 (MQ=255)                                                                       ttGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCa                                                                             <  8:200615/149‑1 (MQ=255)                                                                          ggTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGc                                                                           <  3:288022/148‑1 (MQ=255)                                                                          ggTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCt                                                                          <  6:13541/149‑1 (MQ=255)                                                                           gTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCt                                                                          >  3:288486/1‑148 (MQ=255)                                                                             gCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                                       >  3:161123/1‑149 (MQ=255)                                                                              cccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                                       >  2:251437/1‑148 (MQ=255)                                                                                         ttGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgctgct                                                           >  5:123859/1‑149 (MQ=255)                                                                                           gCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAGACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                          <  4:1372/148‑1 (MQ=255)                                                                                            cACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgca                                                        <  4:288486/149‑1 (MQ=255)                                                                                             aCCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcag                                                       >  7:63537/1‑149 (MQ=255)                                                                                             aCCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcag                                                       >  4:178346/1‑149 (MQ=255)                                                                                             aCCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcag                                                       >  2:102650/1‑149 (MQ=255)                                                                                              ccGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcag                                                       >  2:50238/1‑148 (MQ=255)                                                                                                 gTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACTGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATg                                                   <  8:116850/149‑1 (MQ=255)                                                                                                     cgccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgctg                                                             >  3:230667/1‑135 (MQ=255)                                                                                                     cgccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgctg                                                             <  4:230667/135‑1 (MQ=255)                                                                                                      gccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGGATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCg                                              >  4:177482/1‑149 (MQ=255)                                                                                                      gccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                                       <  5:360746/124‑1 (MQ=255)                                                                                                      gccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                                       >  6:360746/1‑124 (MQ=255)                                                                                                       ccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCCCGCCTCAGCTGCTGCAGCATGCGCCgg                                             >  8:351533/1‑149 (MQ=255)                                                                                                       ccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCgg                                             >  4:42315/1‑149 (MQ=255)                                                                                                           cAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTgc                                         >  6:329515/1‑149 (MQ=255)                                                                                                                  gAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgctgcag                                  >  4:255846/1‑147 (MQ=255)                                                                                                                     aTGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                               <  3:178346/149‑1 (MQ=255)                                                                                                                     aTGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                               <  1:251437/149‑1 (MQ=255)                                                                                                                         ccGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGagc                           >  8:263149/1‑149 (MQ=255)                                                                                                                                 gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGCAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTg                   >  5:266562/1‑149 (MQ=255)                                                                                                                                 gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCt                    <  5:186621/148‑1 (MQ=255)                                                                                                                                 gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTg                   >  7:90380/1‑149 (MQ=255)                                                                                                                                 gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTg                   >  7:280158/1‑149 (MQ=255)                                                                                                                                 gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTg                   <  8:63892/149‑1 (MQ=255)                                                                                                                                 gCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTCGACCAGCAgcgcaa                    <  5:234615/148‑3 (MQ=255)                                                                                                                                   aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGa                 >  5:218493/1‑149 (MQ=255)                                                                                                                                   aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGa                 >  5:338034/1‑149 (MQ=255)                                                                                                                                   aaaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGa                 >  4:200174/1‑149 (MQ=255)                                                                                                                                    aaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGAc                <  8:63537/149‑1 (MQ=255)                                                                                                                                    aaaCGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGAc                <  1:50238/149‑1 (MQ=255)                                                                                                                                       cGGATGATGCCATGGGCATCAACCACCAGGATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCtt             <  3:177482/149‑1 (MQ=255)                                                                                                                                        ggttgatgCCCATGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTg            >  4:114656/4‑149 (MQ=255)                                                                                                                                        ggATGATGGCATGGGCATCAAACACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTg            <  6:107956/149‑1 (MQ=255)                                                                                                                                        ggATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTg            <  3:135601/149‑1 (MQ=255)                                                                                                                                         gttgatgCCCATGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTgg           >  8:128842/3‑149 (MQ=255)                                                                                                                                         gttgatgCCCATGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTgg           >  2:220498/3‑149 (MQ=255)                                                                                                                                         gatgatGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTgg           >  1:56680/1‑149 (MQ=255)                                                                                                                                              tgCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCtt       >  7:324432/1‑148 (MQ=255)                                                                                                                                              tgCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTg      <  3:255846/149‑1 (MQ=255)                                                                                                                                                ccATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGa     <  3:42315/148‑1 (MQ=255)                                                                                                                                                 cATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAg    >  1:370198/1‑148 (MQ=255)                                                                                                                                                 cATGGGCATCAACCACCAGAATGCCTTCGACGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAg    >  4:166138/1‑148 (MQ=255)                                                                                                                                                  aTGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGa  >  2:220902/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGA  >  NC_002947/5650529‑5650822

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