Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,650,675	T→A	6.7%	H173L (CAT→CTT) ‡	PP_4959 ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,650,675	0	T	A	6.7%	219.6 / 3.6	73	H173L (CAT→CTT) ‡	PP_4959	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (40/28);  new base A (2/3);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAT  >  NC_002947/5650527‑5650823                                                                                                                                                      |                                                                                                                                                      ccccTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCGAACGGGAT‑ATGGCCGGATTGGCAAAACGGATGATGCCAt                                                                                                                                                       >  6:332813/1‑149 (MQ=255)       gTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCAt                                                                                                                                                 <  8:216106/149‑1 (MQ=255)        tAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATc                                                                                                                                                <  6:224318/149‑1 (MQ=255)           gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAac                                                                                                                                             <  3:355178/149‑1 (MQ=255)           gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAac                                                                                                                                             >  5:20025/1‑149 (MQ=255)           gCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAac                                                                                                                                             <  5:362205/149‑1 (MQ=255)            ccTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCAtcc                                                                                                                                            <  5:314542/149‑1 (MQ=255)            ccTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGAGTTGCACCGGTGCCGCCAGCAAACGGCAA‑ATCTCCGGATTGGCAAAACGGATGATGCCATGGGCATCAac                                                                                                                                             >  4:7131/1‑148 (MQ=255)             cTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAACATGGCCGGATTGGCAAAACGGATGATGCCCATGGCATCAtcc                                                                                                                                            <  1:270942/149‑1 (MQ=255)                gTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCTCAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGTATGATGCCATGGGCATCAaccacc                                                                                                                                         <  2:97241/148‑1 (MQ=255)                gTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGACATGGTCATCAaacac                                                                                                                                          >  6:266683/1‑147 (MQ=255)                    aaGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACTGATGATGCCATGGGCATCAACCACCAGAAt                                                                                                                                    >  8:322131/1‑149 (MQ=255)                         cGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCtt                                                                                                                              >  3:89086/1‑149 (MQ=255)                         cGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCtt                                                                                                                              >  6:182507/1‑149 (MQ=255)                         cGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCtt                                                                                                                              >  6:341654/1‑149 (MQ=255)                          gACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTc                                                                                                                             <  1:259064/149‑1 (MQ=255)                             tCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTcggc                                                                                                                          >  4:18654/1‑149 (MQ=255)                               gcgcCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTcggcg                                                                                                                         <  6:20025/148‑1 (MQ=255)                                  ccACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCAATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCAt                                                                                                                     >  4:69406/1‑149 (MQ=255)                                         cTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAg                                                                                                              <  3:179930/149‑1 (MQ=255)                                             tACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTcggcg                                                                                                                         >  6:183440/1‑134 (MQ=255)                                             tACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTcggcg                                                                                                                         <  5:183440/134‑1 (MQ=255)                                               cTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGAGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGa                                                                                                        >  1:172407/1‑149 (MQ=255)                                                 ggTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAacca                                                                                                                                           >  4:19031/1‑113 (MQ=255)                                                 ggTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAacca                                                                                                                                           <  3:19031/113‑1 (MQ=255)                                                  gTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACg                                                                                                      >  7:72257/1‑148 (MQ=255)                                                      gTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCAtt                                                                                                  >  5:28269/1‑148 (MQ=255)                                                       ttGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCtt                                                                                                                              >  6:37372/1‑119 (MQ=255)                                                       ttGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCtt                                                                                                                              <  5:37372/119‑1 (MQ=255)                                                             aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATgg                                                                                                         >  6:380422/1‑134 (MQ=255)                                                             aaCGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATgg                                                                                                         <  5:380422/134‑1 (MQ=255)                                                              aCGTCGAGCACATGGGTCCCCTGCTGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGTATGATGCCATGTGCATCAACCACCATAATG‑CCTTCGGCGGCATTATCCAGGATCGACGCATTGACGGCCCt                                                                                         <  2:307767/149‑1 (MQ=255)                                                               cGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCTAGGATGGACGCATTGAAGGCCCTg                                                                                        >  4:150693/1‑149 (MQ=255)                                                                  cGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGcc                                                                                     <  5:130954/149‑1 (MQ=255)                                                                        aTTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCAAGGATGGACGCATTGAAGGCCCTGGCCGtctc                                                                                >  6:49799/1‑148 (MQ=255)                                                                             gTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTc                                                                          <  2:172407/149‑1 (MQ=255)                                                                              tGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCa                                                                         >  5:316145/1‑149 (MQ=255)                                                                                cccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCAc                                                                        >  1:63471/1‑148 (MQ=255)                                                                                 ccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGcc                                                                                     >  6:97014/1‑134 (MQ=255)                                                                                 ccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGcc                                                                                     <  5:97014/134‑1 (MQ=255)                                                                                 ccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGa                                                                      >  4:22235/1‑149 (MQ=255)                                                                                  cTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGAc                                                                     <  1:364370/149‑1 (MQ=255)                                                                                      aGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAg                                                                 <  5:266683/149‑1 (MQ=255)                                                                                           ttGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgctgct                                                            >  4:164353/1‑149 (MQ=255)                                                                                           ttGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAgctgc                                                             <  3:7131/148‑1 (MQ=255)                                                                                             gCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                           >  4:207355/1‑148 (MQ=255)                                                                                             gCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                           >  6:36054/1‑148 (MQ=255)                                                                                             gCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                           <  7:322131/148‑1 (MQ=255)                                                                                              cACCGGTGCCGCCAGCAAACTGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGGATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCtgctgc                                                          <  5:182507/148‑1 (MQ=255)                                                                                              cACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGt                                                                                   <  2:328501/123‑1 (MQ=255)                                                                                              cACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCtgctgc                                                          <  4:89086/148‑1 (MQ=255)                                                                                                 cGGTGCCGCCAGCAAACGTGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTgcagca                                                      <  3:69406/149‑1 (MQ=255)                                                                                                        gccAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCg                                               >  3:355704/1‑149 (MQ=255)                                                                                                         ccAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCACCCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGCAGGCCCTGGCCGTCTCCAGCTCACGACTCCGCTGCTGCAGCATCCGcctg                                              >  4:292699/1‑147 (MQ=255)                                                                                                            gCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTg                                           >  5:71806/1‑149 (MQ=255)                                                                                                                 cGGGAA‑ATGGCCGTATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgctg                                      <  5:22572/149‑1 (MQ=255)                                                                                                                  gggAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATGTCCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgctg                                      >  1:300708/1‑149 (MQ=255)                                                                                                                  gggAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCgctgc                                     >  1:143801/1‑149 (MQ=255)                                                                                                                    gAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGctgctg                                   <  2:63471/149‑1 (MQ=255)                                                                                                                        aTGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                                >  6:340470/1‑149 (MQ=255)                                                                                                                         tGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                                >  1:50073/1‑148 (MQ=255)                                                                                                                          ggCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                                >  2:295446/1‑147 (MQ=255)                                                                                                                             cGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGagca                           <  5:332813/149‑1 (MQ=255)                                                                                                                                     cAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTgg                   >  3:79802/1‑149 (MQ=255)                                                                                                                                     cAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTg                    >  5:174974/1‑148 (MQ=255)                                                                                                                                      aaaaCGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGa                  >  4:226476/1‑149 (MQ=255)                                                                                                                                          tggttgatgCCCATGGCATCATCCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCtt              >  6:341558/5‑149 (MQ=255)                                                                                                                                           ggttgatgCCCATGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                                <  6:154960/127‑1 (MQ=255)                                                                                                                                           ggttgatgCCCATGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTc                                >  5:154960/4‑130 (MQ=255)                                                                                                                                            gatgatGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTgg            >  3:280347/1‑149 (MQ=255)                                                                                                                                                 tgCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCtt        >  2:305095/1‑148 (MQ=255)                                                                                                                                                   ccATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGa      <  8:323460/148‑1 (MQ=255)                                                                                                                                                   ccATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAg     >  5:212166/1‑149 (MQ=255)                                                                                                                                                     aTGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGa   >  8:93430/1‑149 (MQ=255)                                                                                                                                                      tGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAt  >  6:135667/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CCCCTGGTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAT  >  NC_002947/5650527‑5650823

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