Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,025,714	G→T	6.9%	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	6.9%	209.7 / 4.3	67	Q238H (CAG→CAT) ‡	wbpZ	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (27/35);  new base T (3/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

ATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGG  >  NC_002947/2025574‑2025859                                                                                                                                             |                                                                                                                                                   aggagagagcTCTTTTTCCTGTTTGTATGCGCGTTGCGTTACTACAAGGGGCTTGATTACCTGCTAGAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGa                                                                                                                                           <  4:233660/147‑1 (MQ=255) aTGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGa                                                                                                                                           <  4:289352/149‑1 (MQ=255)      cgagcgCTTTTTGCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGATGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTC‑AGCAGCCTGCTGCTTGAGCtg                                                                                                                                       <  2:266821/148‑2 (MQ=255)      cgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCa                                                                                                                                      <  3:303525/149‑1 (MQ=255)      cgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTC‑AGCAGCCTGCTGCTTGAGCtgt                                                                                                                                      <  6:171502/149‑3 (MQ=255)        agcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctg                                                                                                                                                >  7:70084/1‑138 (MQ=255)        agcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctg                                                                                                                                                <  8:70084/138‑1 (MQ=255)          cgcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGCGCTCAGGCt                                                                                                                                  >  7:300660/1‑149 (MQ=255)           gcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGATAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCtt                                                                                                                                 <  8:362050/149‑1 (MQ=255)             tttttCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGa                                                                                                                               >  6:72142/1‑149 (MQ=255)                   cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAg                                                                                                                                          <  7:83540/132‑1 (MQ=255)                    tGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGtt                                                                                                                        >  4:363220/1‑149 (MQ=255)                     gTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGt                                                                                                                         >  7:154961/1‑147 (MQ=255)                     gTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGt                                                                                                                         <  8:154961/147‑1 (MQ=255)                      tttGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCa                                                                                                                      <  5:72142/149‑1 (MQ=255)                      tttGTAGGCGCGTGGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTc                                                                                                                       <  2:299468/148‑1 (MQ=255)                         gTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCCGGCTTGACAATGTTCAGtt                                                                                                                   >  2:332576/1‑149 (MQ=255)                               gcgTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTg                                                                                                               <  1:235930/147‑1 (MQ=255)                               gcgTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTg                                                                                                               >  2:235930/1‑147 (MQ=255)                               gcgTTGCGGTACTACAAGGGCCTGGATTCCCTGCTACAGGCTGCGTCTATCAACCGGCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTggg                                                                                                             <  4:230270/149‑1 (MQ=255)                                    gCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGc                                                                                                         <  7:18019/148‑1 (MQ=255)                                    gCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGc                                                                                                         <  7:23574/148‑1 (MQ=255)                                          ctacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAAGGTTCAGTTCCTGGGCAGCCTGGCa                                                                                                   <  2:321261/148‑1 (MQ=255)                                             cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGAc                                                                                                >  7:339388/1‑148 (MQ=255)                                             cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGAAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGAc                                                                                                >  4:307646/1‑148 (MQ=255)                                               aGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGTCTGGCAGACGAt                                                                                             <  2:160148/149‑1 (MQ=255)                                               aGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGAt                                                                                             <  5:258961/149‑1 (MQ=255)                                               aGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGAt                                                                                             <  5:114770/149‑1 (MQ=255)                                                   gCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGa                                                                                           >  5:118718/1‑147 (MQ=255)                                                           aCCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAgcggc                                                                                  <  6:118718/148‑1 (MQ=255)                                                                      gCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCt                                                                      >  1:144237/1‑149 (MQ=255)                                                                            tCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTa                                                                <  3:289739/149‑1 (MQ=255)                                                                            tCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTa                                                                >  3:267665/1‑149 (MQ=255)                                                                              tATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACag                                                                                        <  7:4939/123‑1 (MQ=255)                                                                              tATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACag                                                                                        >  8:4939/1‑123 (MQ=255)                                                                                tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAgcggcg                                                                                 >  3:98248/1‑128 (MQ=255)                                                                                tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAgcggcg                                                                                 <  4:98248/128‑1 (MQ=255)                                                                                tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                               <  8:339388/146‑1 (MQ=255)                                                                                 cAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                >  3:259259/1‑96 (MQ=255)                                                                                 cAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                                <  4:259259/96‑1 (MQ=255)                                                                                     cGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGa                                                                                                                               >  5:237852/1‑77 (MQ=255)                                                                                     cGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGa                                                                                                                               <  6:237852/77‑1 (MQ=255)                                                                                            cGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTcca                                                >  5:8663/1‑149 (MQ=255)                                                                                              gTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                              <  3:363220/149‑1 (MQ=255)                                                                                                    atcatcGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                                 <  2:351635/124‑1 (MQ=255)                                                                                                    atcatcGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                                 >  1:351635/1‑124 (MQ=255)                                                                                                         cggcggTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGtt                                   >  3:264107/1‑149 (MQ=255)                                                                                                           gcggTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTcc                                 <  6:8663/149‑1 (MQ=255)                                                                                                                   ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGtt                         >  6:138905/1‑149 (MQ=255)                                                                                                                   ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGtt                         >  5:175245/1‑149 (MQ=255)                                                                                                                     ccTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTc                        <  4:267665/148‑1 (MQ=255)                                                                                                                         gAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCCGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGc                     <  1:332576/147‑1 (MQ=255)                                                                                                                         gAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCAt                   >  1:238942/1‑149 (MQ=255)                                                                                                                          aaGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTgg                                                                                                     >  4:5131/1‑66 (MQ=255)                                                                                                                          aaGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTgg                                                                                                     <  3:5131/66‑1 (MQ=255)                                                                                                                            cctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTtc                >  2:11028/5‑149 (MQ=255)                                                                                                                            cctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTtc                >  8:112333/5‑149 (MQ=255)                                                                                                                             ctgaGCTC‑AGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTtct               >  3:190148/4‑149 (MQ=255)                                                                                                                               cAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTctct             <  2:144237/149‑1 (MQ=255)                                                                                                                               cAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTctct             <  5:138905/149‑1 (MQ=255)                                                                                                                                aGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctg            <  7:289547/149‑1 (MQ=255)                                                                                                                                   tCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCAttt                 >  6:75399/1‑141 (MQ=255)                                                                                                                                   tCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCAttt                 <  5:75399/141‑1 (MQ=255)                                                                                                                                    cAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGCCGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgctg         >  2:181247/1‑148 (MQ=255)                                                                                                                                    cAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgctg         >  6:202434/1‑148 (MQ=255)                                                                                                                                     aaGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgc           <  2:27020/145‑1 (MQ=255)                                                                                                                                     aaGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgc           >  1:27020/1‑145 (MQ=255)                                                                                                                                          gcagGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCgg  >  1:92383/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                   ATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGC‑TGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGGAAGCGG  >  NC_002947/2025574‑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: ‑