Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,650,676	G→T	6.8%	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	6.8%	235.7 / 6.1	74	H173N (CAT→AAT) ‡	PP_4959	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (41/28);  new base T (2/3);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.72e-01

GTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAT  >  NC_002947/5650533‑5650823                                                                                                                                                 |                                                                                                                                                     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                                                                                                                         <  5:183440/134‑1 (MQ=255)                                       tACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTcggcg                                                                                                                         >  6:183440/1‑134 (MQ=255)                                         cTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGAGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGa                                                                                                        >  1:172407/1‑149 (MQ=255)                                           ggTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAacca                                                                                                                                           <  3:19031/113‑1 (MQ=255)                                           ggTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAacca                                                                                                                                           >  4:19031/1‑113 (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                                                                                     <  5:97014/134‑1 (MQ=255)                                                                           ccTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGcc                                                                                     >  6:97014/1‑134 (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                                                           >  6:36054/1‑148 (MQ=255)                                                                                       gCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGctgctg                                                           >  4:207355/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)                                                                                                                                                 |                                                                                                                                                     GTAGGCCTTGTAGAAGTCCGACTCGCGCCACAGGCTGGTACTGGTCAGTTGCACAACGTCGAGCAATTGGGTGCCCTGCAGGTGTTGCACCGGTGCCGCCAGCAAACGGGAA‑ATGGCCGGATTGGCAAAACGGATGATGCCATGGGCATCAACCACCAGAATG‑CCTTCGGCGGCATTTTCCAGGATGGACGCATTGAAGGCCCTGGCCGTCTCCAGCTCACGACTCAGCTGCTGCAGCATGCGCCGGTTGCGCTGCTGGTCGAGCAGCGCCTGGACCTTGGGCTTGAGGAT  >  NC_002947/5650533‑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: ‑