Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,451,601	G→C	11.7%	A20G (GCC→GGC)	PP_2147 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,451,601	0	G	C	11.7%	238.9 / 16.3	84	A20G (GCC→GGC)	PP_2147	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (40/34);  new base C (6/4);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTCATCTTCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGTACCAG  >  NC_002947/2451456‑2451749                                                                                                                                                   |                                                                                                                                                     ttCATCTTCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCg                                                                                                                                                    >  2:171073/1‑148 (MQ=255)   cATCTTCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGaa                                                                                                                                                  <  8:310030/148‑1 (MQ=255)   cATCTTCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAg                                                                                                                                                 >  3:307706/1‑149 (MQ=255)    aTCTTCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGCGGCGACAGCGCTACACACTTCGGGGGCGATCGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCAACCTGATACAGGCCATCGGCGAAgg                                                                                                                                                >  8:284740/1‑149 (MQ=255)        tCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGc                                                                                                                                            >  3:168657/1‑149 (MQ=255)        tCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGCCGAAGGcctg                                                                                                                                            >  6:228168/1‑146 (MQ=255)        tCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGCCGAAGGcctg                                                                                                                                            >  6:223056/1‑146 (MQ=255)        tCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCAATTCTACCTGATACAGGCCTTCGCCGAAGGCAGc                                                                                                                                            >  5:51160/1‑149 (MQ=255)          cAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTg                                                                                                                                          >  3:342020/1‑149 (MQ=255)           aGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGCCGAAGGcctgtat                                                                                                                                         <  5:228168/149‑7 (MQ=255)           aGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGCCGAAGGcctgtat                                                                                                                                         <  5:223056/149‑7 (MQ=255)              gCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCTGCTGGGGc                                                                                                                                      <  4:307706/149‑1 (MQ=255)                tGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGaa                                                                                                                                                  >  5:259355/1‑135 (MQ=255)                tGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTATACCTGATACAGGCCTTCGGCGaa                                                                                                                                                  <  6:259355/135‑1 (MQ=255)                 ggggggCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAgg                                                                                                                                                >  4:11239/1‑136 (MQ=255)                 ggggggCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAgg                                                                                                                                                <  3:11239/136‑1 (MQ=255)                     ggCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAgc                                                                                                                               <  4:112633/149‑1 (MQ=255)                        gcTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGACTTCGGCGAAGCCAGCTGGGGCGAACAGCGcc                                                                                                                            >  2:191943/1‑149 (MQ=255)                                  tttCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAggg                                                                                                                  >  2:173541/1‑149 (MQ=255)                                   ttCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGt                                                                                                                 >  7:242629/1‑149 (MQ=255)                                      gCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAg                                                                                                              >  5:201326/1‑149 (MQ=255)                                       ccAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAgg                                                                                                             >  8:290434/1‑149 (MQ=255)                                           gCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAAc                                                                                                         <  4:148262/149‑1 (MQ=255)                                           gCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACACGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAAc                                                                                                         <  5:212819/149‑1 (MQ=255)                                             gTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACg                                                                                                        >  4:257852/1‑148 (MQ=255)                                               ggggCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGAt                                                                                                     <  4:215564/149‑1 (MQ=255)                                                gggCGCCAGCGCTCCTGGCTTCGGGGGCGATCGGGGTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATg                                                                                                    <  5:284250/149‑1 (MQ=255)                                                gggCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATg                                                                                                    <  7:201535/149‑1 (MQ=255)                                                  gcgcCAGCGCTCGAGTCTTCGGGGGCGTACGGGCTGGTGAAGGCTGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGc                                                                                                  <  6:51160/149‑1 (MQ=255)                                                      cAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGCTACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGaa                                                                                                                                   >  1:76454/1‑113 (MQ=255)                                                      cAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGaa                                                                                                                                   <  2:76454/112‑1 (MQ=255)                                                       aGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCa                                                                                             <  1:173541/149‑1 (MQ=255)                                                        gcgcTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCAt                                                                                            <  7:290434/149‑1 (MQ=255)                                                          gcTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGa                                                                                          <  1:171073/149‑1 (MQ=255)                                                          gcTCCAGGCTTCGGGTGCGACCGGGCTGGTGAAGGCGTGC‑ACGCGGTTCTGCCGCACCAGAATCATTTCTACCTGATTCAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGAGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGa                                                                                          <  8:71145/149‑1 (MQ=255)                                                            tCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGt                                                                                        >  3:21001/1‑149 (MQ=255)                                                            tCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGt                                                                                        >  4:89287/1‑149 (MQ=255)                                                            tCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGt                                                                                        >  4:278954/1‑149 (MQ=255)                                                              cAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTc                                                                                       <  3:257852/148‑1 (MQ=255)                                                              cAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCt                                                                                      >  1:222389/1‑149 (MQ=255)                                                               aGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCtt                                                                                     <  6:201326/149‑1 (MQ=255)                                                                   ttCGGGGTCGATCGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCAACCTGATACAGGGCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGTTCAGGCAACGGATGCCACGCATGAGTCTTCCtt                                                                                 <  4:39786/149‑1 (MQ=255)                                                                   ttCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGa                                                                                                      >  3:9168/1‑128 (MQ=255)                                                                   ttCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGa                                                                                                      <  4:9168/128‑1 (MQ=255)                                                                     cGGGGGGGAACGGGCGGGTGAAGGCGGGG‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTTATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCt                                                                                  <  7:148961/146‑1 (MQ=255)                                                                     cGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCt                                                                                  <  8:242629/146‑1 (MQ=255)                                                                     cGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTa                                                                                <  4:21001/148‑1 (MQ=255)                                                                      gggggCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTag                                                                               >  4:61502/1‑148 (MQ=255)                                                                      gggggCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTag                                                                               <  3:61502/148‑1 (MQ=255)                                                                             aaCGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCtg                                                                       >  3:337150/1‑149 (MQ=255)                                                                                      gTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATATCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCCCGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTcgg                                                              >  5:123358/1‑147 (MQ=255)                                                                                              gggC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTa                                                                                >  2:241420/1‑123 (MQ=255)                                                                                              gggC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTa                                                                                <  1:241420/123‑1 (MQ=255)                                                                                               ggC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGaa                                                     >  1:54945/1‑149 (MQ=255)                                                                                                   aCGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAg                                                  >  4:10978/1‑149 (MQ=255)                                                                                                     gCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAgcg                                                <  3:150837/149‑1 (MQ=255)                                                                                                       tGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAgcgc                                               >  1:94007/1‑148 (MQ=255)                                                                                                        gTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCt                                             >  5:149427/1‑149 (MQ=255)                                                                                                         ttCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAgcgc                                               <  2:94007/146‑1 (MQ=255)                                                                                                            tGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGCACAGCGCCTCCAg                                         >  2:265199/1‑149 (MQ=255)                                                                                                            tGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAg                                         >  4:295357/1‑149 (MQ=255)                                                                                                              ccGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAgg                                        <  4:168657/148‑1 (MQ=255)                                                                                                              ccGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCACCTCCAgg                                        <  2:54945/148‑1 (MQ=255)                                                                                                                 cACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCTAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTAt                                    >  5:257836/1‑149 (MQ=255)                                                                                                                 cACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTAt                                    >  6:240943/1‑149 (MQ=255)                                                                                                                     aGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGa                                 >  4:303671/1‑148 (MQ=255)                                                                                                                         tCATTTCTACCTGATACAGGCGTTCGGCGAAGGCAGCTGGGGCGAACAGCGCGAGCAGCAGGGTCAGGCAACGGATGGCAGGGATGAGTCTTCCTGAGGGAGGCTGTGGGGTCAGGGGCTCGGACGGCGCCTCGAGGGTATTGAAACGt                            >  4:288788/1‑149 (MQ=255)                                                                                                                           aTTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCCGGCTGTGGGGTCAGGctct                                                         >  4:260967/1‑118 (MQ=255)                                                                                                                           aTTTCTACCTGAAACAGGCCTTCGGCGAAGGCAGGTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTtc                          <  7:284740/149‑1 (MQ=255)                                                                                                                            tttCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTtcg                         >  8:295966/1‑149 (MQ=255)                                                                                                                                 aCCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGTGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCggg                    >  4:31971/1‑149 (MQ=255)                                                                                                                                     cccagctgCCTTCGCGGAAGGGAGCTGGGGCTAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTATGCAGGCTGTGGGGGCAGGCGCTCGAACAGCGCCTCCAGTGTATTGAAACGTTCTTCGGGGc                  <  4:134634/140‑1 (MQ=255)                                                                                                                                     cccagctgCCTTCGCCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGt                >  1:305079/8‑149 (MQ=255)                                                                                                                                        agctgCCTTCGCCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTc              >  8:320785/5‑148 (MQ=255)                                                                                                                                        aCAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGCGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTc              >  4:349979/1‑148 (MQ=255)                                                                                                                                         gctgCCTTCGCCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAg                                                              <  7:85502/96‑1 (MQ=255)                                                                                                                                         gctgCCTTCGCCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAg                                                              >  8:85502/4‑99 (MQ=255)                                                                                                                                         cAGGCCTTCGGCGAAGGCAGCTTGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCa            <  1:191943/149‑1 (MQ=255)                                                                                                                                          aGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCAt           <  3:303671/149‑1 (MQ=255)                                                                                                                                            gCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCAAtag         >  6:117503/1‑147 (MQ=255)                                                                                                                                             ccTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCgg        <  3:10978/149‑1 (MQ=255)                                                                                                                                               ttCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGt       >  7:61456/1‑148 (MQ=255)                                                                                                                                               ttCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGTa      >  4:165124/1‑149 (MQ=255)                                                                                                                                                 cGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGt       <  8:61456/146‑1 (MQ=255)                                                                                                                                                  ggCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGTAc     >  6:31175/1‑147 (MQ=255)                                                                                                                                                   gCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGTACCAg  <  6:149427/149‑1 (MQ=255)                                                                                                                                                   gCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGTACCAg  >  8:41009/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TTCATCTTCCAGTGCTGGGGGGCGCTCGGCATCTTTCGCCAGGCGTGGGGCGCCAGCGCTCCAGTCTTCGGGGGCGAACGGGCTGGTGAAGGCGGGC‑ACGCTGTTCTGCCGCACCAGAATCATTTCTACCTGATACAGGCCTTCGGCGAAGGCAGCTGGGGCGAACAGCGCCAGCAGCAGGGTCAGGCAACGGATGGCACGCATGAGTCTTCCTTAGGCAGGCTGTGGGGTCAGGCGCTCGAACAGCGCCTCCAGGGTATTGAAACGTTCGTCGGGGCGTTCCATCGGTACCAG  >  NC_002947/2451456‑2451749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑