Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,650,666	A→T	5.7%	I176N (ATC→AAC)	PP_4959 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,650,666	0	A	T	5.7%	211.4 / 2.3	70	I176N (ATC→AAC)	PP_4959	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base T (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

GATCTGCCGCCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTG  >  NC_002947/5650518‑5650812                                                                                                                                                      |                                                                                                                                                    gATCTGCCGCCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACGACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGga                                                                                                                                                     >  6:146329/1‑149 (MQ=255)          ccccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCAt                                                                                                                                           >  4:80102/1‑149 (MQ=255)           cccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGGGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATg                                                                                                                                          <  2:107055/149‑1 (MQ=255)            ccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATgg                                                                                                                                         <  2:279853/149‑1 (MQ=255)            ccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATgg                                                                                                                                         <  6:376601/149‑1 (MQ=255)             cTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACGACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATggg                                                                                                                                        <  5:146329/149‑1 (MQ=255)             cTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATggg                                                                                                                                        >  1:301706/1‑149 (MQ=255)                gTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGCCCGGATTGGCAAAACGGATGATG‑CCATGGGCAt                                                                                                                                     <  8:284420/149‑1 (MQ=255)                 tAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATc                                                                                                                                    <  4:107769/149‑1 (MQ=255)                    gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAac                                                                                                                                 <  5:140401/149‑1 (MQ=255)                     ccTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCacc                                                                                                                                 >  2:306973/1‑146 (MQ=255)                     ccTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CAATGGGCATCAac                                                                                                                                 >  3:344058/1‑148 (MQ=255)                       ttGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAacca                                                                                                                               >  1:247285/1‑149 (MQ=255)                         gTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAaccac                                                                                                                              >  1:162194/1‑147 (MQ=255)                          tAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAAccacca                                                                                                                            <  6:141783/148‑1 (MQ=255)                            gAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCAt                                                                                                                                     <  2:242882/137‑1 (MQ=255)                            gAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGaa                                                                                                                         >  5:367055/1‑149 (MQ=255)                            gAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCAGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCAt                                                                                                                                     >  1:242882/1‑137 (MQ=255)                             aaGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAAt                                                                                                                        <  3:80102/149‑1 (MQ=255)                              aGTCCGAATCTCGACACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGTAAATGGCCGGATTGGCAAAACGTATGATG‑CCATGGGCATCAACCACCAGAAt                                                                                                                        <  4:344058/148‑1 (MQ=255)                                      tCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTcggc                                                                                                               >  4:73702/1‑149 (MQ=255)                                         cgcCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGTTGATG‑CCCATGGCATCAtcc                                                                                                                                <  5:59354/129‑1 (MQ=255)                                         cgcCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGTTGATG‑CCCATGGCATCAtcc                                                                                                                                >  6:59354/1‑129 (MQ=255)                                           ccACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCAt                                                                                                          >  3:17065/1‑149 (MQ=255)                                                 gCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGTTGATGTCCAT‑GGCATCAtccg                                                                                                                               >  2:31616/1‑121 (MQ=255)                                                 gCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                    <  5:335992/149‑1 (MQ=255)                                                 gCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCa                                                                                                    <  2:301706/149‑1 (MQ=255)                                                 gCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGCCAAAACGGTTGATGTCCAT‑GGCATCAtccg                                                                                                                               <  1:31616/122‑2 (MQ=255)                                                   tGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAg                                                                                                   >  2:92043/1‑148 (MQ=255)                                                      tACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATg                                                                                               <  1:126665/149‑1 (MQ=255)                                                         t‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGAc                                                                                            >  3:163510/1‑149 (MQ=255)                                                               aGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCAtt                                                                                       >  7:173576/1‑149 (MQ=255)                                                               aGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCAtt                                                                                       <  8:53569/149‑1 (MQ=255)                                                               aGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCAt                                                                                        <  4:194880/148‑1 (MQ=255)                                                                gTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCAtt                                                                                       >  4:242440/1‑148 (MQ=255)                                                                 ttGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGCCGGCAATTTCCAGGATGGACGCATTg                                                                                      >  6:222889/1‑148 (MQ=255)                                                                     acaACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGAAAAACGGATGATG‑CCATGGGAATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGc                                                                                 >  3:202201/1‑149 (MQ=255)                                                                     acaACGTCGAGCAATTGGCTGACCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGCAGGCATTTTCCAGGATGGACGCATTGAACGc                                                                                 >  7:82196/1‑149 (MQ=255)                                                                         cGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTg                                                                             <  1:92043/149‑1 (MQ=255)                                                                             gAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGcc                                                                          <  4:17065/148‑1 (MQ=255)                                                                             gAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCg                                                                         >  7:13782/1‑149 (MQ=255)                                                                              aGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGt                                                                        >  4:300236/1‑149 (MQ=255)                                                                                      ggTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGTAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCt                                                                <  6:318034/149‑1 (MQ=255)                                                                                        tGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCa                                                              <  2:247285/149‑1 (MQ=255)                                                                                          cccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                             >  6:36881/1‑148 (MQ=255)                                                                                            cTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGAc                                                          >  3:44590/1‑149 (MQ=255)                                                                                              gCAGGTGTTGGACCTGTGCCGGCAGCAATCGGGATATGGCCGGATTGGGAAAACGGATGTTG‑GCATGGGCATGAACCAGCAGTATGCCTTCGGCGGCATTTTCCAGGATGGACGGTTTGAAGTCCCTGGCCGTGTCCAGATCACGACt                                                         <  4:202201/148‑1 (MQ=255)                                                                                                aGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAg                                                      >  3:188375/1‑149 (MQ=255)                                                                                                    gtTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgctgc                                                  >  3:140901/1‑149 (MQ=255)                                                                                                     ttGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGACTTCGGCGGCATTTTCCAGGATGCACGCATTGAACGCCCTGGCCTTCTCCAGCTCACGACTCAgatgct                                                 >  1:90755/1‑149 (MQ=255)                                                                                                       gCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGcc                                                                          <  2:131865/122‑1 (MQ=255)                                                                                                       gCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                >  4:332308/1‑148 (MQ=255)                                                                                                          ccGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcag                                             >  1:241121/1‑148 (MQ=255)                                                                                                                ccgccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATgcgc                                      >  6:329729/1‑149 (MQ=255)                                                                                                                  gccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCg                                    <  5:36301/149‑1 (MQ=255)                                                                                                                   ccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATg                                                                                               >  2:147008/1‑89 (MQ=255)                                                                                                                   ccAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATg                                                                                               <  1:147008/89‑1 (MQ=255)                                                                                                                     aGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGtt                                 <  3:300236/149‑1 (MQ=255)                                                                                                                      gCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTg                                <  6:150621/149‑1 (MQ=255)                                                                                                                        aaaCGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTgcg                              <  4:44590/149‑1 (MQ=255)                                                                                                                          aCGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgct                            <  3:242440/149‑1 (MQ=255)                                                                                                                              gAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGctgctg                        <  8:173576/149‑1 (MQ=255)                                                                                                                               aaaTGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGAATCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTgg                       >  1:42021/1‑149 (MQ=255)                                                                                                                                 aTGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                     <  5:36881/149‑1 (MQ=255)                                                                                                                                    gCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGa                   <  8:13782/148‑1 (MQ=255)                                                                                                                                      cGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGagca                <  4:188375/149‑1 (MQ=255)                                                                                                                                              cAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCTGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTCGACCAGCAgcgcaacc        <  5:69035/149‑5 (MQ=255)                                                                                                                                              cAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTg         >  2:1865/1‑148 (MQ=255)                                                                                                                                               aaaaCGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTggg       >  2:109558/1‑148 (MQ=255)                                                                                                                                                  aCGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCt    >  7:183371/1‑149 (MQ=255)                                                                                                                                                    ggATGATG‑CCATGGGCATCAACCAGCAGAATGCCTTCGGCGGCATTTTGCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCATCTCACGACTCAGGTGCTGAAGCATGCTCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTg  <  8:82196/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GATCTGCCGCCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACT‑GGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAAATGGCCGGATTGGCAAAACGGATGATG‑CCATGGGCATCAACCACCAGAATGCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTG  >  NC_002947/5650518‑5650812

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