Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,030,801	A→T	6.2%	T287S (ACT→TCT)	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,801	0	A	T	6.2%	188.5 / 3.1	65	T287S (ACT→TCT)	PP_0887	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (26/34);  new base T (2/2);  total (29/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 = 7.41e-01

ACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCG  >  NC_002947/1030654‑1030945                                                                                                                                                     |                                                                                                                                                 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                                                                                     >  7:98790/1‑149 (MQ=255)                                                             aTCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGa                                                                                     >  2:130066/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                                                                       <  8:98790/149‑1 (MQ=255)                                                                           ccAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggc                                                                       >  4:182193/1‑149 (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                                                               >  5:230645/1‑125 (MQ=255)                                                                                                           aCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAg                                                               <  6:230645/125‑1 (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            >  5:67544/1‑149 (MQ=255)                                                                                                                                      aGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAg            <  2:324631/149‑1 (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)                                                                                                                                                     |                                                                                                                                                 ACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCG  >  NC_002947/1030654‑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: ‑