Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,025,719:1	+T	8.6%	coding (719/1110 nt)	wbpZ →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,025,719	1	.	T	8.6%	270.3 / 19.3	70	coding (719/1110 nt)	wbpZ	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (32/32);  new base T (4/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATGTA  >  NC_002947/2025571‑2025866                                                                                                                                                      |                                                                                                                                                    cgcaTGAGCGAGCGCTTTTTCCTGTTTGTCGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc                                                                                                                                                      <  7:96091/149‑1 (MQ=255)    aTGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc‑t                                                                                                                                                    >  6:156912/1‑147 (MQ=255)      gagcgagcGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGc                                                                                                                                                >  5:12424/1‑149 (MQ=255)       agcgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCt                                                                                                                                               <  5:174463/149‑1 (MQ=255)       agcgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTC‑AGCAGCCTGCTGCTTGAGCt                                                                                                                                               <  1:283143/149‑1 (MQ=255)         cgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCa                                                                                                                                             <  5:156912/149‑1 (MQ=255)              gcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAAGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCtt                                                                                                                                        >  1:280124/1‑149 (MQ=255)                      cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGt                                                                                                                                >  6:318866/1‑149 (MQ=255)                      cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGACTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGt                                                                                                                                >  3:201980/1‑149 (MQ=255)                       tGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc                                                                                                                                                      <  5:176026/127‑1 (MQ=255)                       tGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc                                                                                                                                                      >  6:176026/1‑127 (MQ=255)                          ttGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCa                                                                                                                             <  5:318866/148‑1 (MQ=255)                           tGTAGGCGCGGTGCGGTACTACAAGGGGCGGGATTACCTGCTACAGGCTGCGTCTATGAAACGTCTGCCGGTAGTGATGATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAg                                                                                                                            <  4:273290/148‑1 (MQ=255)                              aGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTg                                                                                                                                       <  5:128641/134‑1 (MQ=255)                              aGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTg                                                                                                                                       >  6:128641/1‑134 (MQ=255)                               ggCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                       <  3:248619/149‑1 (MQ=255)                               ggCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                       <  1:6097/149‑1 (MQ=255)                                 cgcgGTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTgg                                                                                                                     <  6:117843/149‑1 (MQ=255)                                         ggTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTg                                                                                                             >  8:148398/1‑149 (MQ=255)                                         ggTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTg                                                                                                             >  6:408576/1‑149 (MQ=255)                                             ctacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCa                                                                                                          <  2:280124/148‑1 (MQ=255)                                              tacAAGGGGCTGGATTAACTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGAGCAGCCTGGCAGa                                                                                                        >  5:116625/1‑149 (MQ=255)                                                cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                                           >  7:23328/1‑144 (MQ=255)                                                cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                                           <  8:23328/144‑1 (MQ=255)                                                 aaGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCTGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGa                                                                                                     >  5:37215/1‑149 (MQ=255)                                                 aaGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGa                                                                                                     <  7:248785/149‑1 (MQ=255)                                                         ggATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACagag                                                                                             >  5:60768/1‑149 (MQ=255)                                                         ggATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACagag                                                                                             >  6:35463/1‑149 (MQ=255)                                                                cTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGc                                                                                                                   >  5:211309/1‑120 (MQ=255)                                                                cTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGc                                                                                                                   <  6:211309/120‑1 (MQ=255)                                                                  gCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgctgc                                                                                    >  3:207151/1‑149 (MQ=255)                                                                         gCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCGGGCAGACGATGACAGAGCGGCGCTGCTGGAGCt                                                                             >  4:349910/1‑149 (MQ=255)                                                                              gTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                                        <  4:201980/149‑1 (MQ=255)                                                                                  aTCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTgtgt                                                                         <  6:377702/144‑1 (MQ=255)                                                                                  aTCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTgtgt                                                                         >  5:377702/1‑144 (MQ=255)                                                                                   tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                                      <  5:35463/146‑1 (MQ=255)                                                                                   tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACAATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCa                                                                   <  2:46322/149‑1 (MQ=255)                                                                                     aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgctgct                                                                                   >  6:48040/1‑131 (MQ=255)                                                                                     aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgctgct                                                                                   <  5:48040/131‑1 (MQ=255)                                                                                     aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCtgt                                                                           >  1:255220/1‑139 (MQ=255)                                                                                     aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCtgt                                                                           <  2:255220/139‑1 (MQ=255)                                                                                     aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAtt                                                                 >  3:197176/1‑149 (MQ=255)                                                                                      aCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                                  <  6:391143/147‑1 (MQ=255)                                                                                      aCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                                  >  5:391143/1‑147 (MQ=255)                                                                                      aCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTc                                                                <  2:2160/149‑1 (MQ=255)                                                                                         gTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                       >  4:168257/1‑91 (MQ=255)                                                                                         gTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                       <  3:168257/91‑1 (MQ=255)                                                                                          tCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTAt                                                            >  7:323215/1‑149 (MQ=255)                                                                                              ctggTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTAtt                                                           <  6:37215/144‑1 (MQ=255)                                                                                              ccGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                                      >  2:397862/1‑135 (MQ=255)                                                                                              ccGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                                      <  1:397862/135‑1 (MQ=255)                                                                                               cGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGTTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTcca                                                       >  3:356497/1‑149 (MQ=255)                                                                                               cGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTcca                                                       >  8:351442/1‑149 (MQ=255)                                                                                                ggTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccat                                                      <  8:323215/149‑1 (MQ=255)                                                                                                 gTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGTTGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                                     <  1:302773/149‑1 (MQ=255)                                                                                                              gcggTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGAGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTcc                                        <  6:116625/149‑1 (MQ=255)                                                                                                                ggTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGa                                      <  4:207151/149‑1 (MQ=255)                                                                                                                    gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                                     <  4:49675/130‑1 (MQ=255)                                                                                                                    gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                                     >  3:49675/1‑130 (MQ=255)                                                                                                                      ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGGTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGtt                                >  7:283428/1‑149 (MQ=255)                                                                                                                         cTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                                  >  2:41611/1‑112 (MQ=255)                                                                                                                         cTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                                  <  1:41611/112‑1 (MQ=255)                                                                                                                         cTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCgg                             <  5:408576/149‑1 (MQ=255)                                                                                                                            gAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCa                           <  6:60768/148‑1 (MQ=255)                                                                                                                              gcctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatccca                                                  >  6:41065/6‑123 (MQ=38)                                                                                                                              gcctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatccca                                                  <  5:41065/118‑1 (MQ=38)                                                                                                                                ctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCTTTtct                      >  6:347271/4‑149 (MQ=255)                                                                                                                                ctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTtc                       >  4:204106/4‑148 (MQ=255)                                                                                                                                  cAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCCTTCGGCATTTctct                    <  7:351442/149‑1 (MQ=255)                                                                                                                                     gctcagcagcCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgc                  >  2:319108/5‑149 (MQ=255)                                                                                                                                     cTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgct                 >  4:191897/1‑149 (MQ=255)                                                                                                                                        aaGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTgg               >  5:80418/1‑148 (MQ=255)                                                                                                                                           cagcagGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGc           >  6:202059/1‑149 (MQ=255)                                                                                                                                              cagGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCgg         >  8:309902/1‑148 (MQ=255)                                                                                                                                                 gctgc‑tGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGGATTGGTATTTCCTTCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAAt     <  4:197176/149‑1 (MQ=255)                                                                                                                                                   tgc‑tgAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGTGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATGt   >  2:204626/1‑149 (MQ=255)                                                                                                                                                   tgc‑tgAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATg    >  2:293615/1‑148 (MQ=255)                                                                                                                                                     c‑tgAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATGTa  <  1:293615/148‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGGCAATGTA  >  NC_002947/2025571‑2025866

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