Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,795,556 Δ1 bp 6.5% coding (732/1284 nt) puuB ← gamma‑glutamylputrescine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,795,5560T. 6.5% 198.3 / 7.2 62coding (732/1284 nt)puuBgamma‑glutamylputrescine oxidase
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base . (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GTCAGGCGGTAGTAGTCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATAGATGATGCCACCCAGCGACTCGACC  >  NC_002947/2795413‑2795703
                                                                                                                                               |                                                                                                                                                   
gTCAGGCGGTAGTAGTCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTAcc                                                                                                                                                  <  7:216048/147‑1 (MQ=255)
    ggCGGTAGTAGTCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTg                                                                                                                                            <  3:277855/149‑1 (MQ=255)
     gCGGTAGTAGTCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGc                                                                                                                                           >  4:267052/1‑149 (MQ=255)
     gCGGTAGTAGTCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGc                                                                                                                                           >  6:107793/1‑149 (MQ=255)
      cGGTAGTAGTCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGcc                                                                                                                                          <  5:90015/149‑1 (MQ=255)
       ggTAGTAGGCGAGCAGGTTGTTGCAGTCCTCGACGCTGTAGTCCTGCGGCAGCTGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGGGTGGAGTTGGGGGCCAGGTCGGGTACGAGGTTGgcc                                                                                                                                         >  4:8540/1‑149 (MQ=255)
           gtagtCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAgg                                                                                                                                      >  1:198648/1‑148 (MQ=255)
               tCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGc                                                                                                                                 <  1:121127/149‑1 (MQ=255)
               tCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGc                                                                                                                                 <  3:267052/149‑1 (MQ=255)
               tCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGc                                                                                                                                 >  6:167777/1‑149 (MQ=255)
                cGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCg                                                                                                                                >  8:235949/1‑149 (MQ=255)
                 gAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGt                                                                                                                               >  3:257026/1‑149 (MQ=255)
                      ggTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCg                                                                                                                          >  5:18997/1‑149 (MQ=255)
                          gTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTTGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGccc                                                                                                                                         <  8:246481/130‑1 (MQ=255)
                          gTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGccc                                                                                                                                         >  7:246481/1‑130 (MQ=255)
                          gTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCaa                                                                                                                      <  7:235949/149‑1 (MQ=255)
                              cAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACgat                                                                                                                  <  6:18997/149‑1 (MQ=255)
                                   cTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGatgatg                                                                                                              >  7:180399/1‑148 (MQ=255)
                                   cTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGG‑ACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGaa                                                                                                            >  5:335623/1‑149 (MQ=255)
                                   cTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGG‑ACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGaa                                                                                                            >  3:91631/1‑149 (MQ=255)
                                      gACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACGAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGaa                                                                                                            <  8:274761/147‑1 (MQ=255)
                                      gACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACt                                                                                                          >  2:10906/1‑149 (MQ=255)
                                       aCGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACtt                                                                                                         >  6:97879/1‑149 (MQ=255)
                                           agtagtCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGgcg                                                                                                     >  4:162843/1‑149 (MQ=255)
                                             tagtCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGcgcg                                                                                                   <  5:107793/149‑1 (MQ=255)
                                               gtCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGG‑ACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGa                                                                                                 <  4:91631/148‑1 (MQ=255)
                                               gtCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGG‑ACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGa                                                                                                 <  6:335623/148‑1 (MQ=255)
                                                  cTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCt                                                                                              >  6:54501/1‑149 (MQ=255)
                                                           cagGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCg                                                                                     <  2:292447/149‑1 (MQ=255)
                                                           cagGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCg                                                                                     <  8:180399/149‑1 (MQ=255)
                                                              ggTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGGGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCgggg                                                                                  <  3:8540/149‑1 (MQ=255)
                                                                   gtgcCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGc                                                                             <  2:198648/149‑1 (MQ=255)
                                                                     gcCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTg                                                                                                        <  5:104548/120‑1 (MQ=255)
                                                                     gcCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTg                                                                                                        >  6:104548/1‑120 (MQ=255)
                                                                     gcCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCAc                                                                           >  4:191651/1‑149 (MQ=255)
                                                                               cGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCCGGTTGCCCAGGTAGGCGGTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCTGTGTATGCACGACCGTGtt                                                                  <  2:277162/148‑1 (MQ=255)
                                                                                gCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTgg                                                                >  7:132499/1‑149 (MQ=255)
                                                                                  ccAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCa                                                                                                                                                 <  3:187776/66‑1 (MQ=255)
                                                                                  ccAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCa                                                                                                                                                 >  4:187776/1‑66 (MQ=255)
                                                                                  ccAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCa                                                              >  5:168485/1‑149 (MQ=255)
                                                                                       ggCGCGGAGGGGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCGCGGGTACCAGGTTGCCCAGGTCGGCGGTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACcgcg                                                         <  6:244176/149‑1 (MQ=255)
                                                                                         cTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGtt                                                       <  8:132499/149‑1 (MQ=255)
                                                                                             gtggtgATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGAt                                                   >  6:25109/1‑149 (MQ=255)
                                                                                                  gATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCa                                              >  6:338827/1‑149 (MQ=255)
                                                                                                    tCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACg                                            >  3:156000/1‑149 (MQ=255)
                                                                                                          gCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCggg                                      <  3:191651/149‑1 (MQ=255)
                                                                                                          gCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCggg                                      <  6:168485/149‑1 (MQ=255)
                                                                                                           cGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGTCGCTGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCgg                                       <  5:257272/147‑1 (MQ=255)
                                                                                                           cGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCgg                                       <  5:97879/147‑1 (MQ=255)
                                                                                                           cGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCgg                                       <  5:25109/147‑1 (MQ=255)
                                                                                                            gTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGt                                    >  5:250914/1‑149 (MQ=255)
                                                                                                            gTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGt                                    <  5:54501/149‑1 (MQ=255)
                                                                                                                      gcatggACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATa                           <  4:257026/148‑1 (MQ=255)
                                                                                                                             cTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATAGATGATGc                   >  1:320668/1‑149 (MQ=255)
                                                                                                                                 gcggcCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATAGATGATGCCAcc               <  5:167777/149‑1 (MQ=255)
                                                                                                                                    gcCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGtt                                                                  <  6:204941/95‑1 (MQ=255)
                                                                                                                                    gcCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGtt                                                                  >  5:204941/1‑95 (MQ=255)
                                                                                                                                    gcCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATAGATGATGCCACCCAg            <  1:10906/149‑1 (MQ=255)
                                                                                                                                        gCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATAGATGATGCCACCCAGCGAc        <  3:82777/149‑1 (MQ=255)
                                                                                                                                          tCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATAGATGATGCCACCCAGCGACTc      >  3:244645/1‑149 (MQ=255)
                                                                                                                                            gcgTACCAGGTTGCCCAGGTAGGCGTAGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGAAGCGCACGGCCGGGGTTTGCTCATAGATGATGCCACCCACCGACTCGa    >  8:97975/3‑149 (MQ=255)
                                                                                                                                              gTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTGGCTCATAGCTGATGCCACCCGGTGACTCGccc  >  6:256925/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GTCAGGCGGTAGTAGTCGAGCAGGTAGTTGCAGTCCTCGACGCAGTAGTCCTGCGGCAGCAGGGTGCGTGCCAGTTCTTCGCCCAGTGGCTCGGTGGTGATCACCTGCGTGCCGCATGGCATGGACTTGGCGGCCAGCTCGGGTACCAGGTTGCCCAGGTAGGCGTTGCCGGCAACGATGATGAACTTGGCGCGGACCTTGCCCTGCGGGGTATGCACGACCGGGTTGGCACCGCGTTCGATGCGCACGGCCGGGGTTTGCTCATAGATGATGCCACCCAGCGACTCGACC  >  NC_002947/2795413‑2795703

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