Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,807,246 T→C 7.1% K216E (AAG→GAG)  PP_4234 ← bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,807,2460TC 7.1% 159.9 / 2.3 56K216E (AAG→GAG) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base C (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GCCAGTCCTCGGGCGGGTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTG  >  NC_002947/4807102‑4807393
                                                                                                                                                  |                                                                                                                                                   
gccagTCCTCGGGCGGGTTGTCCAGCTGCTCGGCCATGCGCTTCGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAgg                                                                                                                                                 >  4:236274/1‑149 (MQ=255)
 ccagTCCTCGGGCGGGTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGc                                                                                                                                                >  8:142273/1‑149 (MQ=255)
        tCGGGCGGGTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCAGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGaa                                                                                                                                         >  5:308396/1‑149 (MQ=255)
         cgggcgggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAAc                                                                                                                                        >  2:12530/1‑149 (MQ=255)
          gggcgggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGaa                                                                                                                                         <  5:17442/147‑1 (MQ=255)
          gggcgggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGaa                                                                                                                                         <  3:130839/147‑1 (MQ=255)
            gcgggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTTCGCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTCCAGCCTGAAGAACTc                                                                                                                                      <  1:89089/149‑1 (MQ=255)
             cgggTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAt                                                                                                                                    <  2:98259/149‑1 (MQ=255)
                gTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtcgt                                                                                                                                 >  6:228143/1‑149 (MQ=255)
                       aGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGaa                                                                                                                          >  1:158386/1‑149 (MQ=255)
                          tgctCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACgg                                                                                                                       >  6:220326/1‑149 (MQ=255)
                            ctCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACgg                                                                                                                       <  7:57598/147‑1 (MQ=255)
                                 ccATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGgcg                                                                                                                >  5:349674/1‑149 (MQ=255)
                                   aTGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCa                                                                                                                                                   <  3:116799/112‑1 (MQ=255)
                                   aTGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCa                                                                                                                                                   >  4:116799/1‑112 (MQ=255)
                                            gCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCacccaccc                                                                                                      <  1:12530/148‑1 (MQ=255)
                                            gCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCg                                                                                                     >  2:166993/1‑149 (MQ=255)
                                            gCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCg                                                                                                     >  3:93853/1‑149 (MQ=255)
                                              gTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCa                                                                                                   <  2:358175/149‑1 (MQ=255)
                                              gTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCa                                                                                                   >  6:324705/1‑149 (MQ=255)
                                                  aacaacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTCCAGCCTGAAGAACTccat                                                                                                                                   <  4:267979/113‑4 (MQ=255)
                                                  aacaacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTCCAGCCTGAAGAACTccat                                                                                                                                   >  3:267979/1‑110 (MQ=255)
                                                   acaacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGccacccac                                                                                                        <  5:215968/139‑1 (MQ=255)
                                                   acaacCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGccacccac                                                                                                        >  6:215968/1‑139 (MQ=255)
                                                     aacCCCTGCTCGGTCAGGTCGCCGTGCTGGGTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCg                                                                                            <  6:234834/149‑1 (MQ=255)
                                                        cccTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGaa                                                                                         <  6:308396/149‑1 (MQ=255)
                                                             cTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTg                                                                                     <  5:228143/148‑1 (MQ=255)
                                                              tCGGTCATGTCGCCGTCCTGTTAGCGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGATTGGTGGGGAACAGCATGGAGGTCTTCAGGCTGGAGAAGTCATCGTCGTCGAACGGCGTGGCGCCACCCTCCCGCACGAAGCGGAAGCTGc                                                                                    >  6:62323/1‑148 (MQ=255)
                                                              tCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCAACCACCCGCACGAAGCGGAAGCTGc                                                                                    >  3:49277/1‑148 (MQ=255)
                                                                    aGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAgg                                                                             <  4:93853/149‑1 (MQ=255)
                                                                    aGGTAGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCTCACGAAGCGGAAGAAGCCATAAgg                                                                             >  8:126538/1‑149 (MQ=255)
                                                                       tCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTc                                                                          <  1:166993/149‑1 (MQ=255)
                                                                          ccGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtc                                                                                                                                   <  6:248504/89‑1 (MQ=255)
                                                                          ccGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCAtc                                                                                                                                   >  5:248504/1‑89 (MQ=255)
                                                                             tGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTg                                                                    <  3:236274/149‑1 (MQ=255)
                                                                                 ggTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTa                                                                 >  6:290097/1‑148 (MQ=255)
                                                                                  gTTGGGCTGCTTCGACCACACCGGCTT‑GCGNGTTGAGAATGGTGGGGAACAGCATGGAGTTGTTCAGGCTGGAGAACTCTTCGTCGTCGACCGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTa                                                                 >  1:32592/1‑148 (MQ=255)
                                                                                    tGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAgc                                                             >  5:5372/1‑149 (MQ=255)
                                                                                      ggCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGATCGGCGTGGCGCCGCCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAgcgg                                                           >  6:352494/1‑149 (MQ=255)
                                                                                           ctTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGgcgc                                                        <  6:349674/147‑1 (MQ=255)
                                                                                               gACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTggg                                                   <  5:220326/148‑1 (MQ=255)
                                                                                                                  gTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGc                                                                                    <  5:160582/98‑1 (MQ=255)
                                                                                                                  gTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGc                                                                                    >  6:160582/1‑98 (MQ=255)
                                                                                                                      agaATTGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGTAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCctgct                              <  7:126538/148‑1 (MQ=255)
                                                                                                                      agaATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCctgct                              >  3:134785/1‑148 (MQ=255)
                                                                                                                       gaATGGTGGGGAACAGCATGTAGTTCTTCAGGATGGAGAAGTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGa                            <  2:158386/149‑1 (MQ=255)
                                                                                                                               gggAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCgg                                                                                           >  8:2288/1‑78 (MQ=255)
                                                                                                                               gggAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCgg                                                                                           <  7:2288/78‑1 (MQ=255)
                                                                                                                                   aCAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGc                <  5:290097/149‑1 (MQ=255)
                                                                                                                                     aGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCtt              <  5:324705/149‑1 (MQ=255)
                                                                                                                                      acgatgAGTTCTCCAGCCTGAAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTg             >  2:212453/6‑149 (MQ=255)
                                                                                                                                          ggAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCt         >  7:192752/1‑149 (MQ=255)
                                                                                                                                          ggAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCt         >  6:86656/1‑149 (MQ=255)
                                                                                                                                            aGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGaa                       <  1:339321/133‑1 (MQ=255)
                                                                                                                                            aGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGaa                       >  2:339321/1‑133 (MQ=255)
                                                                                                                                            aGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTg       <  5:86656/149‑1 (MQ=255)
                                                                                                                                              ttcttcAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGc     <  8:192752/149‑1 (MQ=255)
                                                                                                                                                 ttcAGGCTGTAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGc                   <  7:297970/132‑1 (MQ=255)
                                                                                                                                                 ttcAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGc                   >  8:297970/1‑132 (MQ=255)
                                                                                                                                                 ttcAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTg  >  2:145505/1‑149 (MQ=255)
                                                                                                                                                 ttcAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTg  >  8:290018/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCCAGTCCTCGGGCGGGTTGTCCAGCTGCTCGGCCATGCGCTTGGCGTAGAACAACCCCTGCTCGGTCAGGTCGCCGTGCTGGTTGGGCTGCTTCGACCACACCGGCTT‑GCG‑GTTGAGAATGGTGGGGAACAGCATGGAGTTCTTCAGGCTGGAGAACTCATCGTCGTCGAACGGCGTGGCGCCACCCACCCGCACGAAGCGGAAGCTGCCATAAGGGTCGCGCTGGTACAGCGGCGCCTGGGCACCGTAACGGATCGCCTGCTGATTGAACTGCAGCTTGCGCTTGGCCTG  >  NC_002947/4807102‑4807393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑