Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	820,889	G→T	5.1%	R87L (CGT→CTT)	PP_0707 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	820,889	0	G	T	5.1%	321.3 / 4.7	98	R87L (CGT→CTT)	PP_0707	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/50);  new base T (2/3);  total (45/53)
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.96e-01

GGCTGCCGAGGCGCCCGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTGGTA  >  NC_002947/820746‑821037                                                                                                                                                |                                                                                                                                                     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                                                                                                                                            <  3:93591/140‑2 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaa                                                                                                                                            <  3:88509/140‑2 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaa                                                                                                                                            >  4:88509/1‑139 (MQ=255)               ccGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCTTGACGCTGaa                                                                                                                                            >  4:93591/1‑139 (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                                                                                                                        >  6:241723/1‑149 (MQ=255)                          gAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCAc                                                                                                                        >  2:168549/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                                                                                                             <  7:297758/149‑1 (MQ=255)                                     ccTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTgg                                                                                                             >  1:471/1‑149 (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:263266/149‑1 (MQ=255)                                                    cGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGAcc                                                                                              <  3:258351/149‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGaaaa                                                                                                        <  6:133124/137‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGaaaa                                                                                                        >  5:133124/1‑137 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                            <  6:88543/149‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                            <  4:73283/149‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                            <  6:93714/149‑1 (MQ=255)                                                      aCCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAg                                                                                            >  1:173024/1‑149 (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                                                                              >  5:300565/1‑149 (MQ=255)                                                                    aTCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATcgg                                                                              >  3:319373/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                                                                                        <  3:329512/113‑1 (MQ=255)                                                                                              tAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGa                                                                                        >  4:329512/1‑113 (MQ=255)                                                                                                acgacgacgTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTa                                                    <  8:234658/147‑1 (MQ=255)                                                                                                acgacgacgTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTa                                                    <  4:319373/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                  <  6:210950/149‑1 (MQ=255)                                                                                                                                gcgATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCga                  <  3:287599/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)                                                                                                                                                |                                                                                                                                                     GGCTGCCGAGGCGCCCGCACTGGACGAAAACGCCAGCCTCGAACAACTGAGCGACCGCTTGGACCTCATCCGCCAAGGCGTCACCAGCGAGGCTAACGACGACGTGCTGTCGCAACTGCGCCTGGGCGCGATGCAGGTTCAGCGTCAAGCTGATGCATTGAGTGCGCTGCGCACTGCCGATGTGGGAAAACTCGACGACCAGCTGAAGGTCATCGGCCCGGCCCAGCCGGACGAGGCCGCTACCCTGACGCAGCAGCGCAAAGCGCTGGAGGCCGAGAAAAAGGCGCTGGTA  >  NC_002947/820746‑821037

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