Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,806,116 A→T 7.5% 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.5% 182.2 / 3.0 67R592R (CGT→CGA) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base T (2/3);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCT  >  NC_002947/4805972‑4806264
                                                                                                                                                |                                                                                                                                                    
cgaTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGc                                                                                                                                                  <  1:184698/149‑1 (MQ=255)
 gaTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGcc                                                                                                                                                 >  4:180308/1‑149 (MQ=255)
    ggCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAg                                                                                                                                              >  4:210615/1‑149 (MQ=255)
        aCGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgagg                                                                                                                                          <  4:322144/149‑3 (MQ=255)
        aCGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGaaa                                                                                                                                          >  2:108926/1‑149 (MQ=255)
            cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGcc                                                                                                                                                 >  8:34997/1‑138 (MQ=255)
            cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGcc                                                                                                                                                 <  7:34997/138‑1 (MQ=255)
            cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGc                                                                                                                                      <  6:369503/149‑1 (MQ=255)
              gCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCCAACGTCCGAGTAAAACGCt                                                                                                                                     >  6:266467/1‑148 (MQ=255)
                  gCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgag                                                                                                                                           >  1:141541/1‑137 (MQ=255)
                  gCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgag                                                                                                                                           <  2:141541/138‑2 (MQ=255)
                  gCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGac                                                                                                                                  <  4:14173/147‑1 (MQ=255)
                       cccTGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTTTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCaaa                                                                                                                           <  1:276669/149‑1 (MQ=255)
                             gcgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                     <  5:29313/149‑1 (MQ=255)
                             gcgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                     <  8:88529/149‑1 (MQ=255)
                                  ttGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                                <  6:218184/149‑1 (MQ=255)
                                    gcgcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCGATTCGGTGACTGCAGGCTTCACCTCGTCGGCCGAGCGGCCGAGGCAATCGCTGACCACCAAAGCCTGGGAGGTgg                                                                                                              >  8:681/1‑149 (MQ=255)
                                    gcgcAGCGCGGCGGCCTGCTCCGGCGGGCTGATCAGGTACGGGGATCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGATGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTgg                                                                                                              <  8:326087/149‑1 (MQ=255)
                                     cgcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGa                                                                                                             <  5:266467/149‑1 (MQ=255)
                                                 tCCTGCTCCGGCTGGCTGATCAGGTACGGGGTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGcc                                                                                                 <  6:153285/149‑1 (MQ=255)
                                                   cTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGccc                                                                                                <  1:255424/148‑1 (MQ=255)
                                                         cGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCgg                                                                                         <  1:132043/149‑1 (MQ=255)
                                                           ggtggctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGGTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAATGCCTGGGAGGTGGACCTGCCGGTGCCCAGCTAGGTg                                                                                       >  5:10553/3‑149 (MQ=255)
                                                            cggtctgCTCAGTGACGGGGCTCCGCGGGCGCGCAATGGCCGCTCTGGCCATTCGGTGACTGCAGGCTGCACCTCGGCGGCCGCCCGGCCGAGGAAATCGCGGACCACCAAAGCCGGGGAGGTGGACCTGCCGGTGCCCAGCTCGGTGc                                                                                      <  5:306804/145‑1 (MQ=255)
                                                                 tgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGata                                                                                 >  5:219072/1‑149 (MQ=255)
                                                                   aTCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGatat                                                                                <  3:180308/148‑1 (MQ=255)
                                                                   aTCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATg                                                                               >  5:110703/1‑149 (MQ=255)
                                                                        gTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGgcg                                                                            >  4:226113/1‑147 (MQ=255)
                                                                            gggTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCt                                                                      >  3:265424/1‑149 (MQ=255)
                                                                                   cgcgGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCgg                                                               <  6:228211/149‑1 (MQ=255)
                                                                                        gTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATc                                                          >  1:190524/1‑149 (MQ=255)
                                                                                           cTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgc                                                       >  7:295356/1‑149 (MQ=255)
                                                                                            tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgcg                                                      >  5:13167/1‑149 (MQ=255)
                                                                                             cAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgcgc                                                     >  6:11659/1‑149 (MQ=255)
                                                                                              aaTTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgcgc                                                     <  5:11659/148‑1 (MQ=255)
                                                                                              aaTTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCt                                                    <  2:190524/149‑1 (MQ=255)
                                                                                              aaTTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCt                                                    >  5:151618/1‑149 (MQ=255)
                                                                                                  gCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGACGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACAAGCCGGTGCCCAGCTCGGTGCCGATATCGCGCAGGCTGATGCGGCCATCGGCGCTGAAc                                                >  1:281807/1‑149 (MQ=255)
                                                                                                     gCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCg                                             <  1:108926/149‑1 (MQ=255)
                                                                                                      ctctGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGa                                            >  8:188419/1‑149 (MQ=255)
                                                                                                      ctctGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCGGAACTCGa                                            >  7:298116/1‑149 (MQ=255)
                                                                                                         tGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGacac                                         >  4:376042/1‑149 (MQ=255)
                                                                                                         tGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGac                                           <  7:188419/147‑1 (MQ=255)
                                                                                                           gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGa                                            >  6:245793/1‑144 (MQ=255)
                                                                                                           gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGa                                            <  5:245793/144‑1 (MQ=255)
                                                                                                           gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTg                                       <  3:210615/149‑1 (MQ=255)
                                                                                                           gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTg                                       >  4:80340/1‑149 (MQ=255)
                                                                                                           gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTg                                       >  2:353291/1‑149 (MQ=255)
                                                                                                             cATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGc                                     >  6:60099/1‑149 (MQ=255)
                                                                                                               ttCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCa                                   <  6:219072/149‑1 (MQ=255)
                                                                                                                   gTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCgg                               >  3:142448/1‑149 (MQ=255)
                                                                                                                      aCTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGgcgc                            >  2:62212/1‑149 (MQ=255)
                                                                                                                        tGCAGTGTTCAGCTCGTCGGCCGAAGGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCCATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTg                                       <  8:56683/136‑1 (MQ=255)
                                                                                                                        tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCg                                                                                          >  5:124848/1‑85 (MQ=255)
                                                                                                                        tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCg                                                                                          <  6:124848/85‑1 (MQ=255)
                                                                                                                        tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGACCTCGACACTGGc                                     >  3:240969/1‑138 (MQ=255)
                                                                                                                        tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTg                                       <  8:58721/136‑1 (MQ=255)
                                                                                                                        tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTg                                       >  7:58721/1‑136 (MQ=255)
                                                                                                                        tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGACGTGGACATGCCGGTGCCCAGGTGGGTGCCGATATGGGGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTg                                       >  7:56683/1‑136 (MQ=255)
                                                                                                                           aggcTTCAGCTGGTCGGCCGATCGGCCGAGGAAATCGCTGCCCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGGTGATGCGGGCATCGGCGGTGAACTCGTCTCTGGCCATCGGGGCCTCgg                       <  2:281807/146‑1 (MQ=255)
                                                                                                                             tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGaa                     <  1:353291/149‑1 (MQ=255)
                                                                                                                                  aCCTCGTCGGCCGAACGGCCGAGGAAATCGCTGCCCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGACATCGGCGATGAACTCGACACTGGCCATCGGCGCCTCGGAACCTg                 >  6:149885/1‑148 (MQ=255)
                                                                                                                                    cTCGTCGGCCGAACGCCCGAGGAAATCGCTGACCACCACAGCCTGGCAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTAc              <  6:10553/149‑1 (MQ=255)
                                                                                                                                       gtcgGCCGAACGGGCGAGGAAATCGCTGACCACGAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGACCTCGACACTGGCCATCGGCGGCTCGGAACCTGTatcat           >  6:150227/1‑145 (MQ=255)
                                                                                                                                       gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCg            >  7:51588/1‑148 (MQ=255)
                                                                                                                                       ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGCCACTGGCCATCGGCGCCTCGGAACCTGTACCg            >  6:296879/2‑148 (MQ=255)
                                                                                                                                          ggCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGt        >  4:111022/1‑149 (MQ=255)
                                                                                                                                                aCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCt  <  8:295356/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCT  >  NC_002947/4805972‑4806264

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