Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	820,893	A→C	5.0%	Q88H (CAA→CAC) ‡	PP_0707 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	820,893	0	A	C	5.0%	317.8 / 4.1	99	Q88H (CAA→CAC) ‡	PP_0707	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/51);  new base C (2/3);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GGCTGCCGAGGCGCCCGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTGGTAGCC  >  NC_002947/820746‑821040                                                                                                                                                    |                                                                                                                                                    ggcTGCCGAGGCGCCCGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAg                                                                                                                                                    >  4:370373/1‑149 (MQ=255) ggcTGCCGAGGCGCCCGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAg                                                                                                                                                    <  3:198403/149‑1 (MQ=255)    tGCCGAGGCGCCCGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCt                                                                                                                                                  >  3:73283/1‑148 (MQ=255)       cGAGGCGCCCGCACTGGCCGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCTCGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaac                                                                                                                                              <  8:63118/149‑3 (MQ=255)              cccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAg                                                                                                                                       >  3:1181/1‑149 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaa                                                                                                                                               >  4:93591/1‑139 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaa                                                                                                                                               >  4:88509/1‑139 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaa                                                                                                                                               <  3:93591/140‑2 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaa                                                                                                                                               <  3:88509/140‑2 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGt                                                                                                                                      <  5:186288/149‑1 (MQ=255)                  cACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGtgcg                                                                                                                                   <  7:103462/149‑1 (MQ=255)                   aCTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGtgcgc                                                                                                                                  >  6:296633/1‑149 (MQ=255)                   aCTGGAAGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAAAGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGtgcgc                                                                                                                                  >  6:22500/1‑149 (MQ=255)                       gACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGtgcgctgc                                                                                                                               <  8:208616/148‑1 (MQ=255)                       gACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGtgcgctgc                                                                                                                               <  8:122056/148‑1 (MQ=255)                          gAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCAc                                                                                                                           >  2:168549/1‑149 (MQ=255)                          gAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCAc                                                                                                                           >  6:241723/1‑149 (MQ=255)                          gAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCAc                                                                                                                           >  3:267147/1‑149 (MQ=255)                                gCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGa                                                                                                                     <  6:271028/149‑1 (MQ=255)                                 ccAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGAt                                                                                                                    <  5:241723/149‑1 (MQ=255)                                  cAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGAtg                                                                                                                   >  4:45170/1‑149 (MQ=255)                                  cAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGAtg                                                                                                                   >  7:346810/1‑149 (MQ=255)                                   aGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGAtgt                                                                                                                  <  1:254748/149‑1 (MQ=255)                                     ccTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTgg                                                                                                                >  1:471/1‑149 (MQ=255)                                     ccTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTgg                                                                                                                <  7:297758/149‑1 (MQ=255)                                     ccTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTgg                                                                                                                >  7:262175/1‑149 (MQ=255)                                     ccTCGAACAACTGAGCGACCGCGTGGACCTCATCCGCCAAGGCGGCGCCAGCGAGGCTAGCGACGACGTGCTGTCGCAACTGCGCCGGGGCGCGATGCAGGTTCAGCGTCTAGCTGATGCAGTGAGTGCGCTGCGCACTGCCGCTGTgg                                                                                                                >  8:297758/1‑149 (MQ=255)                                      cTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTgg                                                                                                                >  7:233310/1‑148 (MQ=255)                                       tCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTggg                                                                                                               >  5:123790/1‑148 (MQ=255)                                         gAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGaaa                                                                                                            <  2:25138/149‑1 (MQ=255)                                               cTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTTAGTGCGCTGCGCACTGCCGATGTGGGAAAACTcga                                                                                                      <  2:334766/149‑1 (MQ=255)                                                   gCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGAGGACGTGCTGTCGCAACTGCGGCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCgacgac                                                                                                  <  4:1181/149‑1 (MQ=255)                                                    cGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGAcc                                                                                                 <  3:258351/149‑1 (MQ=255)                                                    cGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGAcc                                                                                                 <  3:263266/149‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGaaaa                                                                                                           <  6:133124/137‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGaaaa                                                                                                           >  5:133124/1‑137 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                               <  6:93714/149‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                               >  1:173024/1‑149 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                               <  4:73283/149‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                               <  6:88543/149‑1 (MQ=255)                                                             ggACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGaa                                                                                          <  3:370373/147‑1 (MQ=255)                                                             ggACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAg                                                                                         <  4:267147/148‑1 (MQ=255)                                                               aCCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGt                                                                                       >  3:35369/1‑148 (MQ=255)                                                                 cTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCAt                                                                                    >  5:171944/1‑149 (MQ=255)                                                                 cTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCAt                                                                                    >  5:210950/1‑149 (MQ=255)                                                                 cTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCAt                                                                                    >  1:181869/1‑149 (MQ=255)                                                                  tCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATc                                                                                   >  3:16351/1‑149 (MQ=255)                                                                   cATCCGGCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATcg                                                                                  >  4:287599/1‑149 (MQ=255)                                                                    aTCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATcgg                                                                                 >  3:319373/1‑149 (MQ=255)                                                                    aTCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATcgg                                                                                 >  5:300565/1‑149 (MQ=255)                                                                    aTCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATcg                                                                                  <  8:346810/148‑1 (MQ=255)                                                                    aTCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATc                                                                                   <  5:111630/147‑1 (MQ=255)                                                                               cGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGc                                                                      <  3:45170/149‑1 (MQ=255)                                                                                gTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCTCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGcc                                                                     <  5:22500/149‑1 (MQ=255)                                                                                gTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGcc                                                                     <  2:39476/149‑1 (MQ=255)                                                                                        cGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAgg                                                             >  6:8561/1‑149 (MQ=255)                                                                                        cGAGGCTAACGACGACGGGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAgg                                                             <  4:224013/149‑1 (MQ=255)                                                                                          aGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGcc                                                           <  1:168549/149‑1 (MQ=255)                                                                                           ggCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCg                                                          <  8:262175/149‑1 (MQ=255)                                                                                            gctggtgacgacgTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGc                                                         <  5:370470/143‑1 (MQ=255)                                                                                              tAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGa                                                                                           >  4:329512/1‑113 (MQ=255)                                                                                              tAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGa                                                                                           <  3:329512/113‑1 (MQ=255)                                                                                                acgacgacgTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTa                                                       <  4:319373/147‑1 (MQ=255)                                                                                                acgacgacgTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTa                                                       <  8:234658/147‑1 (MQ=255)                                                                                                      acgTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACg                                               <  5:8561/149‑1 (MQ=255)                                                                                                          gCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACgcag                                            >  2:360297/1‑148 (MQ=255)                                                                                                           cTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACgcagca                                          <  3:216367/149‑1 (MQ=255)                                                                                                           cTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACgcag                                            >  8:295739/1‑147 (MQ=255)                                                                                                           cTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACgcag                                            >  2:19428/1‑147 (MQ=255)                                                                                                              tCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAgcg                                       >  1:297291/1‑149 (MQ=255)                                                                                                                gCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGACGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAgcgc                                      <  8:23313/148‑1 (MQ=255)                                                                                                                gCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAgcgc                                      >  6:164002/1‑148 (MQ=255)                                                                                                                gCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCa                                     <  2:303172/149‑1 (MQ=255)                                                                                                                gCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGACATCGGCCCGGCCCAGCCGGACGAGGCCGCACCCCGGACGCAGCAccgc                                      >  4:59736/1‑148 (MQ=255)                                                                                                                  aaCTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCaaa                                   <  7:295739/149‑1 (MQ=255)                                                                                                                   aCTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGa                                                               >  1:67011/1‑120 (MQ=255)                                                                                                                   aCTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGa                                                               <  2:67011/120‑1 (MQ=255)                                                                                                                      gcgcCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGa                                                 <  5:311004/131‑1 (MQ=255)                                                                                                                      gcgcCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGa                                                 >  6:311004/1‑131 (MQ=255)                                                                                                                          cTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCCGCCGGACGCGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGa                           >  8:175487/1‑149 (MQ=255)                                                                                                                          cTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCTTTGCGCTgct                           <  4:35369/149‑3 (MQ=255)                                                                                                                          cTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTg                             >  1:157503/1‑147 (MQ=255)                                                                                                                          cTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGa                           <  6:300565/149‑1 (MQ=255)                                                                                                                            gggCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAgg                         >  8:83216/1‑149 (MQ=255)                                                                                                                             ggCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGc                        <  6:171944/149‑1 (MQ=255)                                                                                                                               cgcgATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGcc                       >  5:71097/1‑148 (MQ=255)                                                                                                                                gcgATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCga                     <  3:287599/149‑1 (MQ=255)                                                                                                                                gcgATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCga                     <  6:210950/149‑1 (MQ=255)                                                                                                                                 cgATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGGCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCgag                    <  1:346711/149‑1 (MQ=255)                                                                                                                                       aGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAg              >  6:27206/1‑149 (MQ=255)                                                                                                                                         gTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGgc            >  3:97841/1‑149 (MQ=255)                                                                                                                                          ttCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGgcg           >  1:317569/1‑149 (MQ=255)                                                                                                                                            cAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGgcgc          <  5:27206/148‑1 (MQ=255)                                                                                                                                             aGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCt         <  5:296633/148‑1 (MQ=255)                                                                                                                                             aGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTg        <  8:233310/149‑1 (MQ=255)                                                                                                                                               cGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                               >  1:320894/1‑60 (MQ=255)                                                                                                                                               cGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                               <  2:320894/60‑1 (MQ=255)                                                                                                                                               cGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTGGt      <  1:360297/149‑1 (MQ=255)                                                                                                                                                gTCAAGCTTATGCATTGAGGGCGCTGCGCACTGCCGATGGGGGAAAACTCGACGACCACCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTGGTa     <  3:59736/149‑1 (MQ=255)                                                                                                                                                   aaGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTGGTAGcc  >  1:360667/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGCTGCCGAGGCGCCCGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTGGTAGCC  >  NC_002947/820746‑821040

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