Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,030,801	A→T	10.8%	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	10.8%	201.7 / 10.0	74	T287S (ACT→TCT)	PP_0887	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base T (4/4);  total (39/35)
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.71e-01

ACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGG  >  NC_002947/1030654‑1030949                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                               >  1:197181/1‑134 (MQ=255)                      cATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGTCTCTTTGCa                                                                                                                                               <  2:197181/134‑1 (MQ=255)                       aTGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCa                                                                                                                               <  4:67741/149‑1 (MQ=255)                        tGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTT‑GCAAAGTCTCTTTGCaca                                                                                                                                             <  1:207552/135‑3 (MQ=255)                        tGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTT‑GCAAAGTCTCTTTGCaca                                                                                                                                             >  2:207552/1‑133 (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                                                                                        >  4:61274/1‑132 (MQ=255)                                                                               gACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGAt                                                                                        <  3:61274/132‑1 (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                                                                                                                   >  3:265299/1‑73 (MQ=255)                                                                                                               gCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAAc                                                                                                                   <  4:265299/73‑1 (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:222274/149‑1 (MQ=255)                                                                                                                            gACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTa                          >  2:268809/1‑149 (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)                                                                                                                                         tgtt‑gCAAAGTCTCTTTGC‑ACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgcgc             >  4:273445/4‑149 (MQ=255)                                                                                                                                               cAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTgc         <  4:59715/149‑1 (MQ=255)                                                                                                                                                   gagaCTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGa     >  6:78044/1‑149 (MQ=255)                                                                                                                                                     gaCTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACg   <  2:267500/149‑1 (MQ=255)                                                                                                                                                      aCTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACgg  >  7:210288/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑T‑GCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGG  >  NC_002947/1030654‑1030949

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