Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,608,997 C→G 7.6% E201Q (GAG→CAG)  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,9970CG 7.6% 241.2 / 7.9 79E201Q (GAG→CAG) waaA3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand):  ref base C (36/37);  new base G (3/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

GCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAG  >  NC_002947/5608852‑5609144
                                                                                                                                                 |                                                                                                                                                   
gCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCCCACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAgg                                                                                                                                                  >  5:216258/1‑149 (MQ=255)
gCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAgg                                                                                                                                                  <  1:48369/149‑1 (MQ=255)
   tCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGc                                                                                                                                                >  2:96160/1‑148 (MQ=255)
   tCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGc                                                                                                                                                <  1:96160/148‑1 (MQ=255)
    cTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGt                                                                                                                                              >  7:333004/1‑149 (MQ=255)
       cACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGt                                                                                                                                              <  1:261586/146‑1 (MQ=255)
        aCCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCGGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                          <  1:113805/149‑2 (MQ=255)
              gtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGt                                                                                                                                    >  2:193259/1‑149 (MQ=255)
                     cTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGaa                                                                                                                             <  1:191109/149‑1 (MQ=255)
                        gcggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgcg                                                                                                                          >  4:262935/1‑149 (MQ=255)
                          ggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCt                                                                                                                        >  8:319392/1‑149 (MQ=255)
                           gcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTg                                                                                                                       >  2:227810/1‑149 (MQ=255)
                                  acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCg                                                                                                                >  7:80776/1‑149 (MQ=255)
                                  acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCg                                                                                                                >  8:303694/1‑149 (MQ=255)
                                   cacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCgg                                                                                                               >  4:165205/1‑149 (MQ=255)
                                      aGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGttt                                                                                                            >  3:302153/1‑149 (MQ=255)
                                      aGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGttt                                                                                                            >  1:132441/1‑149 (MQ=255)
                                       gggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTc                                                                                                           <  8:333004/149‑1 (MQ=255)
                                        ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  6:155840/149‑1 (MQ=255)
                                        ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  7:303694/149‑1 (MQ=255)
                                        ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                          >  5:323343/1‑116 (MQ=255)
                                        ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                          <  6:323343/117‑2 (MQ=255)
                                        ggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATAGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  8:328946/149‑1 (MQ=255)
                                           gattgccgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                        <  3:346494/144‑1 (MQ=255)
                                            cTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCt                                                                                                      <  5:116588/149‑1 (MQ=255)
                                             tGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                      <  3:108497/116‑1 (MQ=255)
                                             tGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                      >  4:108497/1‑116 (MQ=255)
                                              ggCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTg                                                                                                     >  6:48408/1‑148 (MQ=255)
                                                ccgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGGCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCAc                                                                                                  <  2:19162/149‑1 (MQ=255)
                                                  gccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACgg                                                                                                <  1:193259/149‑1 (MQ=255)
                                                     ccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGTGCTGGGCCTCGGTTTCGGTCTGCACGGCGa                                                                                             <  1:227810/149‑1 (MQ=255)
                                                      cccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                            <  5:48408/149‑1 (MQ=255)
                                                      cccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACGTGCACGCCCTGAGCGCGTGc                                                                                                                                            <  2:239137/101‑1 (MQ=255)
                                                      cccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACGTGCACGCCCTGAGCGCGTGc                                                                                                                                            >  1:239137/1‑101 (MQ=255)
                                                       ccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                            >  1:43902/1‑148 (MQ=255)
                                                        cAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCa                                                                                          <  3:262935/149‑1 (MQ=255)
                                                          aCTGTTCACGCAAAGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                        <  6:216258/149‑1 (MQ=255)
                                                             gTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          <  6:29922/128‑1 (MQ=255)
                                                             gTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                          >  5:29922/1‑128 (MQ=255)
                                                              ttCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCa                                                                                    >  2:18187/1‑149 (MQ=255)
                                                                cACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                                 <  1:192681/118‑1 (MQ=255)
                                                                cACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                                 >  2:192681/1‑118 (MQ=255)
                                                                 aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGc                                                                                               <  2:197865/135‑1 (MQ=255)
                                                                 aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGc                                                                                               >  1:197865/1‑135 (MQ=255)
                                                                      aTGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGccc                                                                            <  6:107760/149‑1 (MQ=255)
                                                                      aTGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGcc                                                                            <  7:319392/149‑1 (MQ=255)
                                                                          ccTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGc                                                                         >  2:315611/1‑148 (MQ=255)
                                                                                 cgTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCg                                                                 <  1:142815/149‑1 (MQ=255)
                                                                                    ggCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGt                                                              >  7:210934/1‑149 (MQ=255)
                                                                                     gcagcaGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGt                                                              <  3:351385/148‑1 (MQ=255)
                                                                                      cagcagCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCa                                                            <  4:302153/149‑1 (MQ=255)
                                                                                        gcagcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGc                                                          >  6:170193/1‑149 (MQ=255)
                                                                                            cTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTgg                                                      >  2:153996/1‑149 (MQ=255)
                                                                                             tGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCgg                                                                    >  6:193404/1‑134 (MQ=255)
                                                                                             tGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCgg                                                                    <  5:193404/134‑1 (MQ=255)
                                                                                                 cgtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAg                                                 >  4:91858/1‑149 (MQ=255)
                                                                                                  gtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgc                                                >  1:352378/1‑149 (MQ=255)
                                                                                                  gtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgc                                                >  8:284605/1‑149 (MQ=255)
                                                                                                         gTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCa                                           >  8:304117/1‑147 (MQ=255)
                                                                                                         gTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATa                                         <  3:91858/149‑1 (MQ=255)
                                                                                                         gTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATa                                         <  5:170193/149‑1 (MQ=255)
                                                                                                          tCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAg                                        <  1:153996/149‑1 (MQ=255)
                                                                                                          tCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAg                                        >  4:342246/1‑149 (MQ=255)
                                                                                                              gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCg                                                                     >  1:154848/1‑116 (MQ=255)
                                                                                                              gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCg                                                                     <  2:154848/116‑1 (MQ=255)
                                                                                                               tCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACg                                   <  8:80776/149‑1 (MQ=255)
                                                                                                                cGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACgg                                  <  6:347111/149‑1 (MQ=255)
                                                                                                                     ttGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                               >  6:86170/1‑147 (MQ=255)
                                                                                                                     ttGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGcc                               >  7:155490/1‑147 (MQ=255)
                                                                                                                           gAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCt                        >  6:228204/1‑148 (MQ=255)
                                                                                                                            aGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCa                      <  3:342246/149‑1 (MQ=255)
                                                                                                                            aGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCa                      <  3:165205/149‑1 (MQ=255)
                                                                                                                               cGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                              <  2:137628/122‑1 (MQ=255)
                                                                                                                               cGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                              >  1:137628/1‑122 (MQ=255)
                                                                                                                                    gggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcg              >  1:373483/4‑149 (MQ=255)
                                                                                                                                         cACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCAt         <  2:132441/149‑1 (MQ=255)
                                                                                                                                         cACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCAt         >  3:81266/1‑149 (MQ=255)
                                                                                                                                           cgcgCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg       <  5:228204/149‑1 (MQ=255)
                                                                                                                                                tCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAg  >  5:73640/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAG  >  NC_002947/5608852‑5609144

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: