Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,030,799	A→T	7.4%	E286V (GAG→GTG) ‡	PP_0887 →	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,030,799	0	A	T	7.4%	186.7 / 3.5	64	E286V (GAG→GTG) ‡	PP_0887	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (26/33);  new base T (2/3);  total (28/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 = 1.00e+00

GCTACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCG  >  NC_002947/1030651‑1030945                                                                                                                                                      |                                                                                                                                                   gCTACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAaga                                                                                                                                                    <  3:154557/149‑1 (MQ=255)    aCCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGAc                                                                                                                                                 >  1:324631/1‑149 (MQ=255)     cccAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAgaga                                                                                                                                                  <  5:328252/147‑1 (MQ=255)      ccAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGAcgt                                                                                                                                               >  4:304344/1‑147 (MQ=255)       cAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACttt                                                                                                                                              <  8:11762/149‑1 (MQ=255)       cAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACttt                                                                                                                                              <  2:61176/149‑1 (MQ=255)           ccgccgGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGTCTCTTTGCac                                                                                                                                          <  3:227606/149‑2 (MQ=255)               cgGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTTGCAAAGTCTCTTTGCacagc                                                                                                                                       <  3:160963/149‑5 (MQ=255)                 gTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTg                                                                                                                                    >  1:325422/1‑149 (MQ=255)                     cgccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgc                                                                                                                                <  3:67762/149‑1 (MQ=255)                     cgccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCAAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgc                                                                                                                                >  6:190776/1‑149 (MQ=255)                      gccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGTACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcg                                                                                                                               <  3:45234/149‑1 (MQ=255)                        ccATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcgcg                                                                                                                             >  5:199369/1‑149 (MQ=255)                                  ggTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcgc                                                                                                                              >  4:52915/1‑138 (MQ=255)                                  ggTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcgc                                                                                                                              <  3:52915/138‑1 (MQ=255)                                  ggTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGc                                                                                                                   <  7:257326/149‑1 (MQ=255)                                  ggTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACGTTGCAACAGCTGGTGCGCGCAGCCGAAg                                                                                                                    <  3:304344/148‑1 (MQ=255)                                    tACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTg                                                                                                                                             <  1:308653/121‑1 (MQ=255)                                    tACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTg                                                                                                                                             >  2:308653/1‑121 (MQ=255)                                    tACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGcc                                                                                                                  >  3:133841/1‑148 (MQ=255)                                     aCGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCaa                                                                                                                >  4:34499/1‑149 (MQ=255)                                              gCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑AGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGcc                                                                                                        <  7:90940/148‑1 (MQ=255)                                               cGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTg                                                                                                      <  3:34499/149‑1 (MQ=255)                                               cGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTg                                                                                                      <  2:40561/149‑1 (MQ=255)                                                gACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTgg                                                                                                     <  6:96021/149‑1 (MQ=255)                                                        gCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTg                                                                                             <  1:306341/149‑1 (MQ=255)                                                           tGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAg                                                                                          <  8:32724/149‑1 (MQ=255)                                                             cTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCa                                                                                        <  2:325422/149‑1 (MQ=255)                                                                aTCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGa                                                                                     >  2:130066/1‑149 (MQ=255)                                                                aTCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGa                                                                                     >  7:98790/1‑149 (MQ=255)                                                                 tCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                                    <  5:62211/149‑1 (MQ=255)                                                                   aaCGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAtgt                                                                                  >  6:265632/1‑149 (MQ=255)                                                                   aaCGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                                    <  6:42207/147‑1 (MQ=255)                                                                   aaCGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                                    >  5:42207/1‑147 (MQ=255)                                                                             gcCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTTCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                                    <  4:247302/137‑1 (MQ=255)                                                                             gcCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAACCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                                    >  3:247302/1‑137 (MQ=255)                                                                              ccAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggc                                                                       >  4:182193/1‑149 (MQ=255)                                                                              ccAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggc                                                                       <  8:98790/149‑1 (MQ=255)                                                                               cAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGTGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggct                                                                      <  4:76049/149‑1 (MQ=255)                                                                                     cacaCCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCa                                                                                        >  8:108836/1‑125 (MQ=255)                                                                                     cacaCCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCa                                                                                        <  7:108836/125‑1 (MQ=255)                                                                                          cGACCAGCTGTTGGATGAGGACCTTCAAATGATCCAGGACCAGGTAAAGCTG‑TGCAAAGTGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGgcgc                                                           <  5:190776/149‑1 (MQ=255)                                                                                                gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGCCTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAAca                                                     >  8:263080/1‑148 (MQ=255)                                                                                                gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCg                                                                                               >  4:61276/1‑107 (MQ=255)                                                                                                gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCg                                                                                               <  3:61276/107‑1 (MQ=255)                                                                                                gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAAcc                                                     >  3:277829/1‑149 (MQ=255)                                                                                                             gaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGAt                                        <  2:206886/149‑1 (MQ=255)                                                                                                              aCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAg                                                               <  6:230645/125‑1 (MQ=255)                                                                                                              aCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAg                                                               >  5:230645/1‑125 (MQ=255)                                                                                                               ccTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgc                                      <  4:133841/149‑1 (MQ=255)                                                                                                                cTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgcg                                     >  5:320977/1‑149 (MQ=255)                                                                                                                  gCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGcc                                   >  5:272634/1‑149 (MQ=255)                                                                                                                   cAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgcgc                                    <  4:277829/147‑1 (MQ=255)                                                                                                                   cAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGccc                                  <  6:199369/149‑1 (MQ=255)                                                                                                                   cAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGccc                                  >  6:94125/1‑149 (MQ=255)                                                                                                                          tccAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGc                            >  2:148018/1‑148 (MQ=255)                                                                                                                            cAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAg                         <  5:265632/149‑1 (MQ=255)                                                                                                                                aCCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTAc                     <  4:273435/149‑1 (MQ=255)                                                                                                                                aCCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTAc                     <  3:182193/149‑1 (MQ=255)                                                                                                                                    ggTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGc                        >  6:69528/1‑142 (MQ=255)                                                                                                                                    ggTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGc                        <  5:69528/142‑1 (MQ=255)                                                                                                                                         aGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAg            <  2:324631/149‑1 (MQ=255)                                                                                                                                         aGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAg            >  5:67544/1‑149 (MQ=255)                                                                                                                                          gctgttgCAAAGTCTCTTTGC‑ACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgc           >  2:64400/6‑149 (MQ=255)                                                                                                                                              ttgCAAAGTCTCTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGcc        >  8:208646/2‑149 (MQ=255)                                                                                                                                                 cAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTgc     <  6:272634/149‑1 (MQ=255)                                                                                                                                                    agagACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGcgcg  >  7:285243/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCTACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCG  >  NC_002947/1030651‑1030945

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