Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,608,994 G→C 8.2% R202G (CGC→GGC)  waaA ← 3‑deoxy‑D‑manno‑octulosonic acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,608,9940GC 8.2% 227.9 / 9.1 74R202G (CGC→GGC) waaA3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (3/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

AGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTG  >  NC_002947/5608850‑5609139
                                                                                                                                                |                                                                                                                                                 
aGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCt                                                                                                                                                 >  3:7563/1‑147 (MQ=255)
  gCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAgg                                                                                                                                             <  5:34774/149‑1 (MQ=255)
        tCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTg                                                                                                                                        >  1:288411/1‑148 (MQ=255)
         cACCGTCGTGGGTGCTGGCGTCAATCCACTCAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGt                                                                                                                                         <  5:253835/146‑1 (MQ=255)
         cACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCTACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGc                                                                                                                                       <  5:281050/148‑1 (MQ=255)
         cACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCg                                                                                                                                          <  5:167941/145‑1 (MQ=255)
           ccGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCAcc                                                                                                                                    >  1:124944/1‑149 (MQ=255)
            cgtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCAcc                                                                                                                                    >  8:162505/1‑148 (MQ=255)
            cgtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaagt                                                                                                                                   <  2:288411/149‑4 (MQ=255)
                  gggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCa                                                                                                                             >  7:147660/1‑149 (MQ=255)
                  gggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCa                                                                                                                             <  7:162505/149‑1 (MQ=255)
                  gggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                     >  6:93674/1‑140 (MQ=255)
                  gggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                     <  5:93674/141‑2 (MQ=255)
                       cTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTTCACGCGCTCAGGGCGTGCACCCAGGTCACGGaa                                                                                                                        <  8:197664/149‑1 (MQ=255)
                                  ccACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACgg                                                                                                                          >  8:129952/1‑136 (MQ=255)
                                  ccACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACgg                                                                                                                          <  7:129952/136‑1 (MQ=255)
                                  ccACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACATGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCt                                                                                                             >  4:138702/1‑149 (MQ=255)
                                   cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                            <  1:203911/149‑1 (MQ=255)
                                     cacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCgg                                                                                                          >  7:149656/1‑149 (MQ=255)
                                          ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCTAACTTGATCGAGCCGGTGAGTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                     <  2:24861/149‑1 (MQ=255)
                                          ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                     >  6:74402/1‑149 (MQ=255)
                                          ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGGTTCg                                                                                                     >  2:201275/1‑149 (MQ=255)
                                          ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGGCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                     <  5:311116/149‑1 (MQ=255)
                                            cgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                   >  2:15364/1‑149 (MQ=255)
                                             gcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTc                                                                                                  >  3:167345/1‑149 (MQ=255)
                                             gcTGGCCGCCCCCAAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTc                                                                                                  >  2:227872/1‑149 (MQ=255)
                                                 gccgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGc                                                                                               <  8:149656/148‑1 (MQ=255)
                                                 gccgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGc                                                                                               <  1:15364/148‑1 (MQ=255)
                                                   cgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACg                                                                                            >  5:48034/1‑149 (MQ=255)
                                                    gccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACgg                                                                                           >  5:92139/1‑149 (MQ=255)
                                                                ttCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCTTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGc                                                                                               <  7:292558/133‑1 (MQ=255)
                                                                ttCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGc                                                                                               >  8:292558/1‑133 (MQ=255)
                                                                 gcACGCACTGCCCGCGCGGGGGGCTGCCGCTGTTCGGCGATGGTCATGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCAt                                                                              <  8:203564/148‑1 (MQ=255)
                                                                   aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCAtct                                                                            >  4:150270/1‑149 (MQ=255)
                                                                      cAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCgg                                                                         <  3:254438/149‑1 (MQ=255)
                                                                                    gTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                           <  4:7563/149‑1 (MQ=255)
                                                                                       gcagcaGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCg                                                            <  7:201437/145‑1 (MQ=255)
                                                                                       gcagcaGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCg                                                            >  8:201437/1‑145 (MQ=255)
                                                                                        cagcagCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCa                                                       <  5:211002/149‑1 (MQ=255)
                                                                                          gcagcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGc                                                     <  3:150270/149‑1 (MQ=255)
                                                                                            agcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGc                                                     <  3:135143/147‑1 (MQ=255)
                                                                                            agcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGc                                                     >  4:135143/1‑147 (MQ=255)
                                                                                             gcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTg                                                  <  4:141403/149‑1 (MQ=255)
                                                                                             gcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTg                                                  >  5:216845/1‑149 (MQ=255)
                                                                                              cTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTgg                                                 <  2:124944/149‑1 (MQ=255)
                                                                                               tGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCggg                                                          <  6:77063/139‑1 (MQ=255)
                                                                                               tGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCggg                                                          >  5:77063/1‑139 (MQ=255)
                                                                                                gTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGc                                                >  4:229496/1‑148 (MQ=255)
                                                                                                    gtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGc                                                >  5:313830/1‑144 (MQ=255)
                                                                                                    gtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGc                                                <  6:313830/144‑1 (MQ=255)
                                                                                                    gtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgc                                           >  2:222416/1‑149 (MQ=255)
                                                                                                       gATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                         <  4:167345/148‑1 (MQ=255)
                                                                                                        agggTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGcc                                        >  1:198429/3‑148 (MQ=255)
                                                                                                           gTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATa                                    <  5:70777/149‑1 (MQ=255)
                                                                                                             cAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCACTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGc                                  >  7:252015/1‑149 (MQ=255)
                                                                                                                gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGcc                                        <  6:152447/140‑1 (MQ=255)
                                                                                                                gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGcc                                        >  5:152447/1‑140 (MQ=255)
                                                                                                                gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCAc                               >  6:2375/1‑149 (MQ=255)
                                                                                                                gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCAc                               <  8:147660/149‑1 (MQ=255)
                                                                                                                   gAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACggg                            <  6:92139/149‑1 (MQ=255)
                                                                                                                    aaCTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGc                           <  2:2148/149‑1 (MQ=255)
                                                                                                                      cTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                          >  3:243122/1‑148 (MQ=255)
                                                                                                                       ttGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCGTCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                          >  1:137051/1‑147 (MQ=255)
                                                                                                                       ttGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCACTATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGa                        <  1:227872/149‑1 (MQ=255)
                                                                                                                           tCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcgc                    >  3:7517/1‑149 (MQ=255)
                                                                                                                           tCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcgc                    >  8:76773/1‑149 (MQ=255)
                                                                                                                                ccGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGa               >  8:11336/1‑149 (MQ=255)
                                                                                                                                  ggTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGAc              >  5:202678/1‑148 (MQ=255)
                                                                                                                                      aCTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcg         <  1:222416/149‑1 (MQ=255)
                                                                                                                                       ggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAg                       <  7:43800/132‑1 (MQ=255)
                                                                                                                                       ggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAg                       >  8:43800/3‑134 (MQ=255)
                                                                                                                                        ttGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgcg       >  1:218327/1‑149 (MQ=255)
                                                                                                                                            aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                         <  2:235122/111‑1 (MQ=255)
                                                                                                                                            aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                         >  1:235122/1‑111 (MQ=255)
                                                                                                                                             cgcgCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg  <  4:7517/149‑1 (MQ=255)
                                                                                                                                             cgcgCTCAGGGCGTGCACCCAGGGCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg  <  3:229496/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
AGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTG  >  NC_002947/5608850‑5609139

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: