Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	782,840	T→G	6.0%	E1029A (GAG→GCG)	PP_0672 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	782,840	0	T	G	6.0%	213.5 / 3.8	67	E1029A (GAG→GCG)	PP_0672	sensory box protein
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base G (2/2);  total (31/36)
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.54e-01

CAGTTCCTCGATCACCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGG  >  NC_002947/782695‑782980                                                                                                                                                   |                                                                                                                                             cAGTTCCTCGATCACCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGCTTCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTc                                                                                                                                             >  1:110624/1‑148 (MQ=255)   gTTCCTCGATCACCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCcatg                                                                                                                                         <  5:40613/149‑1 (MQ=255)      ccTCGATCACCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCcatggca                                                                                                                                      >  5:202763/1‑149 (MQ=255)              aCCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAg                                                                                                                              >  8:177134/1‑149 (MQ=255)              aCCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAg                                                                                                                              >  5:56098/1‑149 (MQ=255)                cGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGt                                                                                                                            >  3:205020/1‑149 (MQ=255)                      gAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCcatgg                                                                                                                                        >  8:353576/1‑131 (MQ=255)                      gAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCcatgg                                                                                                                                        <  7:353576/131‑1 (MQ=255)                       aaGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTc                                                                                                                     >  5:53875/1‑149 (MQ=255)                        aGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTGGCCGTCGCCGCTGAACGGCGGCTGGTAAGAGAGGCGGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCGTTcc                                                                                                                    >  5:97087/1‑149 (MQ=255)                         gTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTc                                                                                                                     >  3:253326/1‑147 (MQ=255)                             cccGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTc                                                                                                               <  6:134385/149‑1 (MQ=255)                              ccGGCGGCACCAGGCCTCGGGTGGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTc                                                                                                               >  7:57454/1‑148 (MQ=255)                                  cggcACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGgc                                                                                                          >  7:117156/1‑149 (MQ=255)                                   ggcACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGgcg                                                                                                         >  7:187115/1‑149 (MQ=255)                                     cACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCt                                                                                                       >  5:372466/1‑149 (MQ=255)                                        cAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCGCcat                                                                                                                                          >  7:153589/1‑111 (MQ=255)                                        cAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCGCcat                                                                                                                                          <  8:153589/111‑1 (MQ=255)                                          ggCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAgg                                                                                                  <  6:29798/149‑1 (MQ=255)                                          ggCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAgg                                                                                                  <  6:24373/149‑1 (MQ=255)                                                     gggtGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGtt                                                                                       <  1:22482/149‑1 (MQ=255)                                                      ggtGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTc                                                                                      <  3:244101/149‑1 (MQ=255)                                                       gtGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCa                                                                                     >  4:10154/1‑149 (MQ=255)                                                           cGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTc                                                                                 <  6:372466/149‑1 (MQ=255)                                                            gCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCg                                                                                >  8:311788/1‑149 (MQ=255)                                                             ccAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCt                                                                                                                                              >  6:332273/1‑86 (MQ=255)                                                             ccAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCt                                                                                                                                              <  5:332273/86‑1 (MQ=255)                                                              cAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAg                                                                                                   >  2:226512/1‑128 (MQ=255)                                                              cAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAg                                                                                                   <  1:226512/128‑1 (MQ=255)                                                              cAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGc                                                                              <  8:57454/149‑1 (MQ=255)                                                                    cagcagCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGt                                                                         <  8:287406/149‑1 (MQ=255)                                                                     agcagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCg                                                                                >  5:158657/1‑141 (MQ=255)                                                                     agcagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCg                                                                                <  6:158657/141‑1 (MQ=255)                                                                      gcagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAgg                                                                                                  >  4:273355/1‑122 (MQ=255)                                                                      gcagcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAgg                                                                                                  <  3:273355/122‑1 (MQ=255)                                                                        agcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGaa                                                                     <  6:56098/149‑1 (MQ=255)                                                                        agcGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGaa                                                                     >  6:347451/1‑149 (MQ=255)                                                                         gcgcTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATAGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAAc                                                                    >  7:192455/1‑149 (MQ=255)                                                                           gcTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTg                                                                  >  7:241962/1‑149 (MQ=255)                                                                           gcTTCGGCGCAGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTg                                                                  >  4:104431/1‑149 (MQ=255)                                                                               cGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGaa                                                               <  6:53875/148‑1 (MQ=255)                                                                               cGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAg                                                              >  8:199789/1‑149 (MQ=255)                                                                                    ccGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAgttgtt                                                         <  5:347451/149‑1 (MQ=255)                                                                                              gcTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCgcgc                                               <  8:187115/149‑1 (MQ=255)                                                                                                        gccgcTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCg                                      <  7:311788/148‑1 (MQ=255)                                                                                                        gccgcTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGt                                     <  8:117156/149‑1 (MQ=255)                                                                                                        gccgcTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGt                                     <  7:65421/149‑1 (MQ=255)                                                                                                        gccgcTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGt                                     <  4:253326/149‑1 (MQ=255)                                                                                                         ccgcTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTg                                    <  3:104431/149‑1 (MQ=255)                                                                                                           gcTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTgg                                   >  2:345141/1‑148 (MQ=255)                                                                                                               aaCTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACAt                              <  1:345141/149‑1 (MQ=255)                                                                                                                aCTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACAt                              <  2:110624/148‑1 (MQ=255)                                                                                                                 cTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATgg                            <  6:202763/149‑1 (MQ=255)                                                                                                                    cGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCg                          <  1:363485/148‑1 (MQ=255)                                                                                                                      gCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGgtg                       <  7:199789/149‑1 (MQ=255)                                                                                                                      gCTGGTAATAGAGGAGGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGgtg                       <  3:348728/149‑1 (MQ=255)                                                                                                                          gTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCgg                   <  3:336231/149‑1 (MQ=255)                                                                                                                             aTAGAGGATGAATTCGTTCTGCTCCATAGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCg                 <  8:192455/148‑1 (MQ=255)                                                                                                                                    aTGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTtca          >  7:179799/1‑148 (MQ=255)                                                                                                                                         ttCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTtcatca       >  3:313227/1‑146 (MQ=255)                                                                                                                                         ttCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTtcatca       <  4:313227/146‑1 (MQ=255)                                                                                                                                          tCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTg   >  8:200215/1‑149 (MQ=255)                                                                                                                                            gTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAgttgtt                                                         >  8:118563/1‑93 (MQ=255)                                                                                                                                            gTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAgttgtt                                                         <  7:118563/93‑1 (MQ=255)                                                                                                                                            gTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTgg  <  7:177134/148‑1 (MQ=255)                                                                                                                                             ttCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGt                                  >  7:244182/1‑115 (MQ=255)                                                                                                                                             ttCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGt                                  <  8:244182/115‑1 (MQ=255)                                                                                                                                              ccTGCGCCATGCCATGGAGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACt                                                                   <  2:36562/80‑1 (MQ=39)                                                                                                                                              ccTGCGCCATGCCATGGAGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACt                                                                   >  1:36562/2‑81 (MQ=39)                                                                                                                                                   |                                                                                                                                             CAGTTCCTCGATCACCGGGATGAAGTCGCCCGGCGGCACCAGGCCTCGGGTCGGGTGACGCCAGC‑GCAGCAGCGCTTCGGCGCCGGTCAGGCGCTTGCCGTCGCCGCTGAACTGCGGCTGGTAATAGAGGATGAATTCGTTCTGCTCCATGGCATGGCGCAGGTCGCTTTCCAGCTCCAGGCGCTCCAGGGCGCTGGCGTTCATTTCGGCCTGGTAGAACTGGAAGTTGTTCTTGCCGCGCTCCTTGGCGTGGTACATGGCGGTGTCGGCGTTCTTCATCAGCTGG  >  NC_002947/782695‑782980

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