Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,201,109	C→G	5.6%	A136P (GCG→CCG)	PP_4582 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,201,109	0	C	G	5.6%	272.5 / 3.9	83	A136P (GCG→CCG)	PP_4582	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (35/43);  new base G (3/2);  total (38/45)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GGTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACAT  >  NC_002947/5200964‑5201253                                                                                                                                                   |                                                                                                                                                  ggTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCt                                                                                                                                                 <  2:243522/148‑1 (MQ=255) ggTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCt                                                                                                                                                 <  6:160816/148‑1 (MQ=255) ggTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTg                                                                                                                                                >  8:30225/1‑149 (MQ=255)  gTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTgg                                                                                                                                               <  6:269150/149‑1 (MQ=255)  gTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTgg                                                                                                                                               <  8:369631/149‑1 (MQ=255)     gggACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGt                                                                                                                                             <  7:290568/148‑1 (MQ=255)     gggACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGt                                                                                                                                             <  2:79751/148‑1 (MQ=255)      ggACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGc                                                                                                                                           >  2:225051/1‑149 (MQ=255)        aCGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCg                                                                                                                                         <  2:388235/149‑1 (MQ=255)        aCGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCGGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGcc                                                                                                                                          >  1:60397/1‑148 (MQ=255)              gCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTgg                                                                                                                                   <  3:52449/149‑1 (MQ=255)               cACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGt                                                                                                                                  <  6:15017/149‑1 (MQ=255)                   cATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATc                                                                                                                              >  4:271841/1‑149 (MQ=255)                    aTGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCAt                                                                                                                               <  2:321765/147‑1 (MQ=255)                    aTGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATcc                                                                                                                             >  7:79706/1‑149 (MQ=255)                    aTGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATcc                                                                                                                             <  5:279908/149‑1 (MQ=255)                        ccAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCgg                                                                                                                         <  8:205238/149‑1 (MQ=255)                          agcagTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGcc                                                                                                                       >  6:7660/1‑149 (MQ=255)                             agTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGc                                                                                                                    <  2:49147/149‑1 (MQ=255)                                  acacaTAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACg                                                                                                               >  8:205977/1‑149 (MQ=255)                                  acacaTAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACg                                                                                                               >  8:200463/1‑149 (MQ=255)                                  acacaTAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACg                                                                                                               >  4:5324/1‑149 (MQ=255)                                   cacaTAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGc                                                                                                              >  1:167402/1‑149 (MQ=255)                                    acaTAGGCGTGTCTCCGTTTACAGACGGGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGc                                                                                                              >  1:223414/1‑149 (MQ=255)                                    acaTAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGcc                                                                                                             >  1:235117/1‑149 (MQ=255)                                            gtgtCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTg‑gcggcg                                                                                                    >  4:217873/1‑149 (MQ=255)                                            gggtCTCCGTTTACAGAC‑GGGGCTCGAGCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCGCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTg‑gcggcg                                                                                                    <  2:199169/147‑1 (MQ=255)                                              gtCTCCGTTTACAGAC‑GGGGCTCGGCCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTg‑gcggcggc                                                                                                  <  8:240828/149‑1 (MQ=255)                                               tctcCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTg‑gcggcggc                                                                                                  <  8:79706/148‑1 (MQ=255)                                               tctcCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCt                                                                                                 >  3:308536/1‑149 (MQ=255)                                                ctcCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCtt                                                                                                >  3:137896/1‑149 (MQ=255)                                                 tcCGTTTACAGAC‑GGGGCTCGACCCGGTCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTc                                                                                               >  4:295789/1‑149 (MQ=255)                                                  ccGTTTACAGAC‑GGGGCTAGACCCGTCCTGCGACACGAGAAGCGGCTGGCGGCAAGTCTTCACTGTAACGTTCGCTGCGCTCGGCAGCCAGCGGCTCCTGGGTGCCGAGCTGGTCATCCTCGGACCTCTGACGCCTG‑GCGGCGGCTTc                                                                                               <  2:235117/148‑1 (MQ=255)                                                   cGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCgg                                                                                             >  2:175148/1‑149 (MQ=255)                                                               ggggCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACgg                                                                                  >  7:168878/1‑149 (MQ=255)                                                                 ggCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgcg                                                                                <  1:175148/149‑1 (MQ=255)                                                                  gCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCACCCAGGCGCCGCTGGCTGcc                                                                                                                                          <  3:167733/91‑1 (MQ=37)                                                                  gCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCACCCAGGCGCCGCTGGCTGcc                                                                                                                                          >  4:167733/1‑91 (MQ=37)                                                                   cTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCt                                                                                                                   <  2:333033/113‑1 (MQ=255)                                                                   cTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCt                                                                                                                   >  1:333033/1‑113 (MQ=255)                                                                      gACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgcgtgcgt                                                                           <  7:30225/149‑1 (MQ=255)                                                                      gACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgcgtgc                                                                             <  7:288740/147‑1 (MQ=255)                                                                       aCCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTGTGCGGCGGCTTCGGCCTGCTCACGgcgtgcgt                                                                           <  1:241559/149‑1 (MQ=255)                                                                       aCCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgcgtgcctg                                                                          >  2:194963/1‑149 (MQ=255)                                                                         ccGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgcgtgcgtgcg                                                                        >  6:162577/1‑149 (MQ=255)                                                                         ccGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgcgtgcgtgcg                                                                        >  7:12891/1‑149 (MQ=255)                                                                         ccGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGACTTCGGCCTGCTCACGgcgtgcgtgcg                                                                        >  3:335461/1‑149 (MQ=255)                                                                          cGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgcgtgcgtgcgt                                                                       <  7:315172/149‑1 (MQ=255)                                                                            gCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCgtgcgtgcgtgc                                                                     <  3:5324/149‑1 (MQ=255)                                                                            gCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCgtgcgtgcgtgc                                                                     <  3:217873/149‑1 (MQ=255)                                                                              cTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgc                                                                                 >  7:106414/1‑135 (MQ=255)                                                                              cTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGgc                                                                                 <  8:106414/135‑1 (MQ=255)                                                                              cTGCGCCACGCGCAGAGGTTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTgcgc                                                                   <  4:173103/149‑1 (MQ=255)                                                                               tGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCt                                                                  >  2:27718/1‑149 (MQ=255)                                                                                  gcCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCtttt                                                               <  1:225051/149‑1 (MQ=255)                                                                                    cACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTcc                                                             <  2:167402/149‑1 (MQ=255)                                                                                       gcgcAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCt                                                                                                 <  2:64692/110‑1 (MQ=255)                                                                                       gcgcAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCt                                                                                                 >  1:64692/1‑110 (MQ=255)                                                                                        cgcAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCg                                                         <  2:223414/149‑1 (MQ=255)                                                                                              ggCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAgc                                                   <  7:200463/149‑1 (MQ=255)                                                                                              ggCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAgc                                                   <  7:205977/149‑1 (MQ=255)                                                                                               gCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAgcg                                                  <  2:60397/149‑1 (MQ=255)                                                                                               gCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAgcg                                                  >  8:344962/1‑149 (MQ=255)                                                                                                           tcgtcACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGaa                                      <  5:7660/149‑1 (MQ=255)                                                                                                             gtcACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAag                                     <  3:271841/148‑1 (MQ=255)                                                                                                             gtcACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAag                                     <  3:295789/148‑1 (MQ=255)                                                                                                                  ggTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAg                               >  7:271862/1‑149 (MQ=255)                                                                                                                      aCGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTa                             >  7:175277/1‑147 (MQ=255)                                                                                                                        ggTCGCTGCGCTCGGCAGCCAGCGGCGGCTGGGTGCCGAGCTGGTCATCCTCTGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCTCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAgg                         <  6:375655/149‑1 (MQ=255)                                                                                                                            gCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGACTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGa                     <  4:335461/149‑1 (MQ=255)                                                                                                                              tGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAg                   >  8:87436/1‑149 (MQ=255)                                                                                                                               cagcTCGGCACCCAGGCGCCGCTGGCTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAgg                  >  7:341560/3‑149 (MQ=255)                                                                                                                                cgcTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAgg                  <  4:137896/148‑1 (MQ=255)                                                                                                                                  cTCGGCACCCAGGCGCCGCTGGCTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGa               >  7:270397/1‑149 (MQ=255)                                                                                                                                  cTCGACAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGa               <  2:169306/149‑1 (MQ=255)                                                                                                                                     ggCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATc             <  2:390236/148‑1 (MQ=255)                                                                                                                                     ggCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATc             >  1:390236/1‑148 (MQ=255)                                                                                                                                      gCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAgg                                  >  4:138043/1‑126 (MQ=255)                                                                                                                                      gCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAgg                                  <  3:138043/126‑1 (MQ=255)                                                                                                                                      gCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCAt           <  1:27718/149‑1 (MQ=255)                                                                                                                                           cagcGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTAACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAACTACAGGCCGAAGGCGATCATGCCg       >  1:133567/1‑148 (MQ=255)                                                                                                                                            agcGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTa     >  5:237183/1‑149 (MQ=255)                                                                                                                                               ggcgCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACAt  >  5:294366/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GGTAGGGACGTTGGCACTCATGCCCAGCAGTTCACACATAGGCGTGTCTCCGTTTACAGAC‑GGGGCTCGACCCGGCCTGCGCCACGCGCAGCGGCTGGCGGCAAGTCGTCACGGTAACGGTCGCTGCGCTCGGCAGCCAGCGGCGCCTGGGTGCCGAGCTGGTCATCCTCGGCCCGCTGACGCCTG‑GCGGCGGCTTCGGCCTGCTCACGGCGTGCGTGCGTGCGCTTTTCCAGCGGCCAGCGCAGCAGCACGAAGAGGAAGTACAGGCCGAAGGCGATCATGCCGTACAT  >  NC_002947/5200964‑5201253

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