Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,025,706 Δ1 bp 6.7% coding (706/1110 nt) wbpZ → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,025,7060A. 6.7% 236.4 / 10.7 75coding (706/1110 nt)wbpZglycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base . (3/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

ACTGGCGCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGG  >  NC_002947/2025560‑2025853
                                                                                                                                                  |                                                                                                                                                    
aCTGGCGCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAag                                                                                                                                                    >  8:248785/1‑149 (MQ=255)
 cTGGCGCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAagc                                                                                                                                                   <  3:91458/149‑1 (MQ=255)
      gcgcACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGcagcag                                                                                                                                              <  6:294970/149‑1 (MQ=255)
           cgcaTGAGCGAGCGCTTTTTCCTGTTTGTCGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc                                                                                                                                         <  7:96091/149‑1 (MQ=255)
              aTGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgct                                                                                                                                       >  6:156912/1‑147 (MQ=255)
                gagcgagcGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGc                                                                                                                                   >  5:12424/1‑149 (MQ=255)
                 agcgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCt                                                                                                                                  <  5:174463/149‑1 (MQ=255)
                 agcgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAGCAGCCTGCTGCTTGAGCt                                                                                                                                  <  1:283143/149‑1 (MQ=255)
                   cgagcgCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCa                                                                                                                                <  5:156912/149‑1 (MQ=255)
                        gcTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAAGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCtt                                                                                                                           >  1:280124/1‑149 (MQ=255)
                                cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGt                                                                                                                   >  6:318866/1‑149 (MQ=255)
                                cTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGACTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGt                                                                                                                   >  3:201980/1‑149 (MQ=255)
                                 tGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc                                                                                                                                         <  5:176026/127‑1 (MQ=255)
                                 tGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGgctgc                                                                                                                                         >  6:176026/1‑127 (MQ=255)
                                    ttGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCa                                                                                                                <  5:318866/148‑1 (MQ=255)
                                     tGTAGGCGCGGTGCGGTACTACAAGGGGCGGGATTACCTGCTACAGGCTGCGTCTATGAAACGTCTGCCGGTAGTGATGATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAg                                                                                                               <  4:273290/148‑1 (MQ=255)
                                        aGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTg                                                                                                                          >  6:128641/1‑134 (MQ=255)
                                        aGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTg                                                                                                                          <  5:128641/134‑1 (MQ=255)
                                         ggCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                          <  3:248619/149‑1 (MQ=255)
                                         ggCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                          <  1:6097/149‑1 (MQ=255)
                                          gcgcGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCaa                                                                                                                                                     <  8:142443/106‑1 (MQ=255)
                                          gcgcGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCaa                                                                                                                                                     >  7:142443/1‑106 (MQ=255)
                                           cgcgGTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTgg                                                                                                        <  6:117843/149‑1 (MQ=255)
                                                   ggTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTg                                                                                                >  6:408576/1‑149 (MQ=255)
                                                   ggTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTg                                                                                                >  8:148398/1‑149 (MQ=255)
                                                       ctacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCa                                                                                             <  2:280124/148‑1 (MQ=255)
                                                        tacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAag                                                                                                                                                    >  7:35098/1‑93 (MQ=255)
                                                        tacAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAag                                                                                                                                                    <  8:35098/93‑1 (MQ=255)
                                                        tacAAGGGGCTGGATTAACTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGAGCAGCCTGGCAGa                                                                                           >  5:116625/1‑149 (MQ=255)
                                                          cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                              <  8:23328/144‑1 (MQ=255)
                                                          cAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGc                                                                                              >  7:23328/1‑144 (MQ=255)
                                                           aaGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCTGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGa                                                                                        >  5:37215/1‑149 (MQ=255)
                                                           aaGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGa                                                                                        <  7:248785/149‑1 (MQ=255)
                                                                   ggATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACagag                                                                                >  5:60768/1‑149 (MQ=255)
                                                                   ggATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACagag                                                                                >  6:35463/1‑149 (MQ=255)
                                                                          cTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGc                                                                                                      >  5:211309/1‑120 (MQ=255)
                                                                          cTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGc                                                                                                      <  6:211309/120‑1 (MQ=255)
                                                                            gCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgctgc                                                                       >  3:207151/1‑149 (MQ=255)
                                                                                   gCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCGGGCAGACGATGACAGAGCGGCGCTGCTGGAGCt                                                                >  4:349910/1‑149 (MQ=255)
                                                                                        gTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGtt                                                           <  4:201980/149‑1 (MQ=255)
                                                                                            aTCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTgtgt                                                            <  6:377702/144‑1 (MQ=255)
                                                                                            aTCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTgtgt                                                            >  5:377702/1‑144 (MQ=255)
                                                                                             tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                         <  5:35463/146‑1 (MQ=255)
                                                                                             tCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACAATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCa                                                      <  2:46322/149‑1 (MQ=255)
                                                                                               aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgctgct                                                                      <  5:48040/131‑1 (MQ=255)
                                                                                               aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCgctgct                                                                      >  6:48040/1‑131 (MQ=255)
                                                                                               aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCtgt                                                              <  2:255220/139‑1 (MQ=255)
                                                                                               aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCtgt                                                              >  1:255220/1‑139 (MQ=255)
                                                                                               aaCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAtt                                                    >  3:197176/1‑149 (MQ=255)
                                                                                                aCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                     <  6:391143/147‑1 (MQ=255)
                                                                                                aCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                     >  5:391143/1‑147 (MQ=255)
                                                                                                aCCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTc                                                   <  2:2160/149‑1 (MQ=255)
                                                                                                   gTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                          >  4:168257/1‑91 (MQ=255)
                                                                                                   gTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCt                                                                                                          <  3:168257/91‑1 (MQ=255)
                                                                                                    tCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTAt                                               >  7:323215/1‑149 (MQ=255)
                                                                                                        ctggTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTAtt                                              <  6:37215/144‑1 (MQ=255)
                                                                                                        ccGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                         >  2:397862/1‑135 (MQ=255)
                                                                                                        ccGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTAt                                                         <  1:397862/135‑1 (MQ=255)
                                                                                                         cGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGTTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTcca                                          >  3:356497/1‑149 (MQ=255)
                                                                                                         cGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTcca                                          >  8:351442/1‑149 (MQ=255)
                                                                                                          ggTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccat                                         <  8:323215/149‑1 (MQ=255)
                                                                                                           gTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGTTGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                        <  1:302773/149‑1 (MQ=255)
                                                                                                                        gcggTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGAGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTcc                           <  6:116625/149‑1 (MQ=255)
                                                                                                                          ggTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGa                         <  4:207151/149‑1 (MQ=255)
                                                                                                                              gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                        >  3:49675/1‑130 (MQ=255)
                                                                                                                              gCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatc                                        <  4:49675/130‑1 (MQ=255)
                                                                                                                                ttCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGGTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGtt                   >  7:283428/1‑149 (MQ=255)
                                                                                                                                   cTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                     >  2:41611/1‑112 (MQ=255)
                                                                                                                                   cTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCAt                                                     <  1:41611/112‑1 (MQ=255)
                                                                                                                                   cTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCgg                <  5:408576/149‑1 (MQ=255)
                                                                                                                                      gAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCa              <  6:60768/148‑1 (MQ=255)
                                                                                                                                        gcctgaGCTCAGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatccca                                     <  5:41065/118‑1 (MQ=38)
                                                                                                                                        gcctgaGCTCAGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTccatccca                                     >  6:41065/6‑123 (MQ=38)
                                                                                                                                          ctgaGCTCAGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCTTTtct         >  6:347271/4‑149 (MQ=255)
                                                                                                                                          ctgaGCTCAGCAGCCTGCTGCTTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTtc          >  4:204106/4‑148 (MQ=255)
                                                                                                                                            cAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCCTTCGGCATTTctct       <  7:351442/149‑1 (MQ=255)
                                                                                                                                               cTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTctgct    >  4:191897/1‑149 (MQ=255)
                                                                                                                                                  aaGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTgg  >  5:80418/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACTGGCGCGCACGCATGAGCGAGCGCTTTTTCCTGTTTGTAGGCGCGTTGCGGTACTACAAGGGGCTGGATTACCTGCTACAGGCTGCGTCTATCAACCGTCTGCCGGTAGTGATCATCGGCGGTGGCTTCCTCGAAGCGCAGCTCAAGCAGCAGGCTGCTGAGCTCAGGCTTGACAATGTTCAGTTCCTGGGCAGCCTGGCAGACGATGACAGAGCGGCGCTGCTGGAGCTGTGTTATGCATTCGTATTTCCATCCCATCTGCGTTCCGAGTCGTTCGGCATTTCTCTGCTGG  >  NC_002947/2025560‑2025853

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: