Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,025,714	G→T	8.3%	Q238H (CAG→CAT) ‡	wbpZ →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,025,714	0	G	T	8.3%	225.9 / 8.6	72	Q238H (CAG→CAT) ‡	wbpZ	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (32/34);  new base T (4/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

GCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGG  >  NC_002947/2025566‑2025859                                                                                                                                                     |                                                                                                                                                   gcgcACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGcagcag                                                                                                                                                    <  6:294970/149‑1 (MQ=255)      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                                                                                                                                >  6:128641/1‑134 (MQ=255)                                   aGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTg                                                                                                                                <  5:128641/134‑1 (MQ=255)                                    ggCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                <  1:6097/149‑1 (MQ=255)                                    ggCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                <  3:248619/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                                                                                                    <  8:23328/144‑1 (MQ=255)                                                     cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                                    >  7:23328/1‑144 (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                                                               <  1:397862/135‑1 (MQ=255)                                                                                                   ccGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                               >  2:397862/1‑135 (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)                                                                                                                                                     |                                                                                                                                                   GCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGG  >  NC_002947/2025566‑2025859

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