Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,030,791:1	+T	5.6%	coding (849/1353 nt)	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,791	1	.	T	5.6%	296.2 / 7.9	72	coding (849/1353 nt)	PP_0887	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base . (36/32);  new base T (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

ACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCG  >  NC_002947/1030654‑1030937                                                                                                                                           |                                                                                                                                                    aCCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGAc                                                                                                                                         >  6:169634/1‑149 (MQ=255)   ccAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACtt                                                                                                                                       <  8:10584/149‑1 (MQ=255)     aGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTg                                                                                                                                     >  4:235467/1‑149 (MQ=255)     aGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTg                                                                                                                                     >  3:304081/1‑149 (MQ=255)           ccgGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAg                                                                                                                               <  7:5203/149‑1 (MQ=255)            cgGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGc                                                                                                                              <  4:201462/149‑1 (MQ=255)              gTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGTCTCTTTGC‑ACAGCttt                                                                                                                           <  3:188237/149‑3 (MQ=255)               tGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTT‑GCAAAGTCTCTTTGC‑ACAGCttt                                                                                                                           <  7:124118/149‑3 (MQ=255)                  cgccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTgc                                                                                                                        <  4:131388/149‑1 (MQ=255)                    cccATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTgcgc                                                                                                                      >  8:85446/1‑149 (MQ=255)                    cccATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTgcgc                                                                                                                      <  3:310230/149‑1 (MQ=255)                      cATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGTCTCTTTGCa                                                                                                                                   <  2:197181/134‑1 (MQ=255)                      cATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGTCTCTTTGCa                                                                                                                                   >  1:197181/1‑134 (MQ=255)                       aTGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCa                                                                                                                   <  4:67741/149‑1 (MQ=255)                        tGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTT‑GCAAAGTCTCTTTGCaca                                                                                                                                 >  2:207552/1‑133 (MQ=255)                        tGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTT‑GCAAAGTCTCTTTGCaca                                                                                                                                 <  1:207552/135‑3 (MQ=255)                         gAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGc                                                                                                                 >  5:73981/1‑149 (MQ=255)                            cTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCg                                                                                                               >  1:241916/1‑148 (MQ=255)                               ggTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGc                                                                                                           >  2:248251/1‑149 (MQ=255)                                 tACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGcc                                                                                                          >  2:247642/1‑148 (MQ=255)                                    gccgcTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAAcc                                                                                                      >  1:322359/1‑149 (MQ=255)                                        cTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGgcg                                                                                                  <  4:221975/149‑1 (MQ=255)                                           gCGACCATGAGCGTGCTGATAAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCt                                                                                               >  1:306630/1‑149 (MQ=255)                                              aCCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGc                                                                                            <  6:73981/149‑1 (MQ=255)                                                cATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCa                                                                                          <  2:74982/149‑1 (MQ=255)                                                     gCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTg                                                                                     <  1:248251/149‑1 (MQ=255)                                                      cGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTgg                                                                                    <  3:235467/149‑1 (MQ=255)                                                        tGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAg                                                                                  <  7:85446/149‑1 (MQ=255)                                                           tGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACa                                                                                                                                >  4:32323/1‑100 (MQ=255)                                                           tGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACa                                                                                                                                <  3:32323/100‑1 (MQ=255)                                                              tCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                            >  7:16135/1‑149 (MQ=255)                                                                aaCGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAtgt                                                                          >  3:110562/1‑149 (MQ=255)                                                                aaCGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAtgt                                                                          >  8:74540/1‑149 (MQ=255)                                                                 aCGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAtgtg                                                                         <  4:304081/149‑1 (MQ=255)                                                                         cgcCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctg                                                                 <  4:110562/149‑1 (MQ=255)                                                                          gcCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctgg                                                                >  5:118774/1‑149 (MQ=255)                                                                           ccAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggc                                                               <  2:306630/149‑1 (MQ=255)                                                                               gACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                            <  3:61274/132‑1 (MQ=255)                                                                               gACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                            >  4:61274/1‑132 (MQ=255)                                                                                aCCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGAc                                                          >  6:28179/1‑149 (MQ=255)                                                                                  cacaCCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGAc                                                          >  6:297167/1‑147 (MQ=255)                                                                                  cacaCCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGa                                                        <  2:241916/149‑1 (MQ=255)                                                                                    cacCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAgg                                                      <  2:664/149‑1 (MQ=255)                                                                                     acCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAgg                                                      >  1:267500/1‑149 (MQ=255)                                                                                     acCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCACGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGgc                                                     >  7:287672/1‑149 (MQ=255)                                                                                            cGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAAc                                              >  2:18837/1‑149 (MQ=255)                                                                                            cGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTCGCCTTCGTGGAGCAGGATCTGACCGCCTGGCTGGACGAGGCGCTCAAc                                              >  8:289073/1‑149 (MQ=255)                                                                                             gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAAcc                                             >  8:192059/1‑149 (MQ=255)                                                                                             gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAAcc                                             <  1:247642/149‑1 (MQ=255)                                                                                                   gCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTgg                                        >  1:222274/1‑148 (MQ=255)                                                                                                     aggaggACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCAc                                     >  3:211566/1‑149 (MQ=255)                                                                                                      ggaggaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCAcc                                    >  5:236741/1‑149 (MQ=255)                                                                                                        aggaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTg                                  <  7:192059/149‑1 (MQ=255)                                                                                                          gaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGAt                                <  5:169634/149‑1 (MQ=255)                                                                                                           aCCTGCAAATACTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGAt                                <  8:16135/148‑1 (MQ=255)                                                                                                             cTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgcg                             <  4:211566/149‑1 (MQ=255)                                                                                                               gCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAAc                                                                                                       <  4:265299/73‑1 (MQ=255)                                                                                                               gCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAAc                                                                                                       >  3:265299/1‑73 (MQ=255)                                                                                                                cAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGccc                          <  1:7902/149‑1 (MQ=255)                                                                                                                 aaaTCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCg                         >  2:82612/1‑149 (MQ=255)                                                                                                                    tcctccAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCgg                        >  6:269320/1‑147 (MQ=255)                                                                                                                    tcctccAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGaa                      <  7:74540/149‑1 (MQ=255)                                                                                                                       tccAGTACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGc                    >  2:267107/1‑148 (MQ=255)                                                                                                                       tccAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGc                    >  6:328428/1‑148 (MQ=255)                                                                                                                         cAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCa                  <  5:328428/148‑1 (MQ=255)                                                                                                                            gACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTa              >  2:268809/1‑149 (MQ=255)                                                                                                                            gACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTa              <  2:222274/149‑1 (MQ=255)                                                                                                                            gACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTCATGCGCCCGGAAGCCAGcaa              >  2:276786/1‑147 (MQ=255)                                                                                                                                   tAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTc       <  5:269320/149‑1 (MQ=255)                                                                                                                                      aGCTGTT‑GCAAAGTCTCTTTGC‑ACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAg    >  1:251342/1‑149 (MQ=255)                                                                                                                                       gCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCACAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgc   >  8:300111/1‑149 (MQ=255)                                                                                                                                        cTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgcg  >  6:332532/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                                    ACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCG  >  NC_002947/1030654‑1030937

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