Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,806,116 A→T 7.1% R592R (CGT→CGA 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,806,1160AT 7.1% 219.2 / 5.3 75R592R (CGT→CGAPP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/35);  new base T (2/3);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

ATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCT  >  NC_002947/4805974‑4806264
                                                                                                                                              |                                                                                                                                                     
aTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCg                                                                                                                                                 <  4:281678/149‑1 (MQ=255)
 tGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAAc                                                                                                                                                      <  6:226083/143‑1 (MQ=255)
 tGGAGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAAc                                                                                                                                                      >  5:226083/1‑143 (MQ=255)
   gCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGac                                                                                                                                               <  7:112014/148‑2 (MQ=255)
    cGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAgg                                                                                                                                              >  2:310492/1‑148 (MQ=255)
        gCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAAt                                                                                                                                         >  1:217681/1‑148 (MQ=255)
        gCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAAt                                                                                                                                         <  2:217681/148‑1 (MQ=255)
        gCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATc                                                                                                                                        <  6:308539/149‑1 (MQ=255)
        gCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATc                                                                                                                                        <  7:239676/149‑1 (MQ=255)
         ccAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCg                                                                                                                                       <  8:207298/149‑1 (MQ=255)
             cTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGa                                                                                                                                   <  2:49869/149‑1 (MQ=255)
                            cgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGcccacgag                                                                                                                                            >  8:42931/1‑120 (MQ=255)
                            cgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgag                                                                                                                                           <  7:42931/126‑2 (MQ=255)
                                ttGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                                <  7:158709/149‑1 (MQ=255)
                                ttGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                                <  1:310492/149‑1 (MQ=255)
                                    gcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGac                                                                                                                                               <  7:201652/115‑2 (MQ=255)
                                    gcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGac                                                                                                                                               >  8:201652/1‑114 (MQ=255)
                                    gcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGa                                                                                                             >  8:122439/1‑148 (MQ=255)
                                     cAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACa                                                                                                           <  4:12982/149‑1 (MQ=255)
                                       gcgcGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATg                                                                                                         <  6:325920/149‑1 (MQ=255)
                                               tCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa                                                                                                                                                >  2:308078/1‑103 (MQ=255)
                                               tCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa                                                                                                                                                <  1:308078/103‑1 (MQ=255)
                                                ccTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGcc                                                                                                 >  6:94119/1‑148 (MQ=255)
                                                 cTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCa                                                                                               >  6:37595/1‑149 (MQ=255)
                                                   gCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGc                                                                                             >  4:164766/1‑149 (MQ=255)
                                                         gctggctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTg                                                                                       >  2:294633/1‑149 (MQ=255)
                                                          ctggctgATCAGGTACGGGTTTCCTCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGc                                                                                      <  5:6845/149‑1 (MQ=255)
                                                             gctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCTGGTGCCCAGCTCGGTGCCGa                                                                                   >  6:129788/1‑149 (MQ=255)
                                                              ctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGat                                                                                  <  5:73625/149‑1 (MQ=255)
                                                                     ggTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGgcgc                                                                           <  5:169120/149‑1 (MQ=255)
                                                                      gTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCa                                                                          <  1:123669/149‑1 (MQ=255)
                                                                           ggTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTg                                                                     <  5:94119/149‑1 (MQ=255)
                                                                            gTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGc                                                                                                        <  5:252531/113‑1 (MQ=255)
                                                                            gTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGc                                                                                                        >  6:252531/1‑113 (MQ=255)
                                                                             tttCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGAt                                                                   >  5:167064/1‑149 (MQ=255)
                                                                                 cGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCgg                                                               <  3:164766/149‑1 (MQ=255)
                                                                                      gTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCAt                                                           <  1:294633/148‑1 (MQ=255)
                                                                                       tGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCg                                                         <  2:17444/149‑1 (MQ=255)
                                                                                         cTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgc                                                       >  1:239058/1‑149 (MQ=255)
                                                                                          tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTGGGGCGAACGGGCGAGGAAATCGCTGTGCACCAAAGCCTTGGAGGTTGACATGCCGGTGGCCAGCTCGGTGCCGATATGgc                                                                             <  5:168361/126‑1 (MQ=255)
                                                                                          tCAAGTGCAGCTCTGGCCATTCGGTGACTGGAGTCTTCAGCTCGTCGGCCGAGCGGCCGAGGAAATCGCTGACGCCCAAAGCCTGGGACGTGGAGATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCGGGCTGATGCGGCCCTCGgcg                                                      >  8:175972/1‑149 (MQ=255)
                                                                                            aaTTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGCAAATCGCTGACCACCAAAGCCTGGGAGGTGGa                                                                                                             <  1:39010/93‑1 (MQ=255)
                                                                                            aaTTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGa                                                                                                             >  2:39010/1‑92 (MQ=255)
                                                                                              ttGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTg                                                   >  5:217555/1‑148 (MQ=255)
                                                                                                 cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCTGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACt                                               <  5:129788/149‑1 (MQ=255)
                                                                                                    ctctGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGAGCTGAACTCGa                                            >  7:301407/1‑149 (MQ=255)
                                                                                                      ctGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGaca                                          <  7:122439/149‑1 (MQ=255)
                                                                                                          ccATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGCGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGGCACTgg                                      >  4:213258/1‑149 (MQ=255)
                                                                                                                 gTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCgg                               <  2:99652/149‑1 (MQ=255)
                                                                                                                    aCTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCt                                                    >  6:198541/1‑125 (MQ=255)
                                                                                                                      tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTgg                                      >  3:218452/1‑137 (MQ=255)
                                                                                                                      tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTgg                                      <  4:218452/137‑1 (MQ=255)
                                                                                                                       gCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATg                                                                               <  5:8824/95‑1 (MQ=255)
                                                                                                                       gCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATg                                                                               >  6:8824/1‑95 (MQ=255)
                                                                                                                        cAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCg                        >  6:314476/1‑149 (MQ=255)
                                                                                                                           tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              <  7:299804/124‑1 (MQ=255)
                                                                                                                           tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              >  8:299804/1‑124 (MQ=255)
                                                                                                                           tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGaa                     <  2:132442/149‑1 (MQ=255)
                                                                                                                           tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGaa                     >  8:256897/1‑149 (MQ=255)
                                                                                                                           tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGaa                     >  6:35817/1‑149 (MQ=255)
                                                                                                                             ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAgg                                                                        <  1:252201/96‑1 (MQ=255)
                                                                                                                             ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAgg                                                                        >  2:252201/1‑96 (MQ=255)
                                                                                                                             ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAcc                   >  7:73904/1‑149 (MQ=255)
                                                                                                                             ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAcc                   >  1:124981/1‑149 (MQ=255)
                                                                                                                               cACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTg                 >  2:235178/1‑149 (MQ=255)
                                                                                                                                aCCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGt                <  3:270395/149‑1 (MQ=255)
                                                                                                                                aCCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGAAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTg                 >  3:103295/1‑148 (MQ=255)
                                                                                                                                 ccTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTa               >  4:21980/1‑149 (MQ=255)
                                                                                                                                  cTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTAc              <  6:217555/149‑1 (MQ=255)
                                                                                                                                  cTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTAc              <  2:239058/149‑1 (MQ=255)
                                                                                                                                     gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTTGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCg            <  5:37595/148‑1 (MQ=255)
                                                                                                                                      tcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCg            >  5:313016/1‑147 (MQ=255)
                                                                                                                                          ccGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTcc      >  8:150378/1‑149 (MQ=255)
                                                                                                                                          ccGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTcc      >  5:190634/1‑149 (MQ=255)
                                                                                                                                           cGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCt     >  1:255152/1‑149 (MQ=255)
                                                                                                                                             aaCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTc   <  5:314476/149‑1 (MQ=255)
                                                                                                                                             aaCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTc   <  7:150378/149‑1 (MQ=255)
                                                                                                                                              aCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCt  <  6:167064/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
ATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCT  >  NC_002947/4805974‑4806264

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: