Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,808,082 A→T 5.3% S593T (TCG→ACG)  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,808,0820AT 5.3% 1124.1 / 25.7 331S593T (TCG→ACG) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (157/154);  new base T (5/14);  total (162/169)
Fisher's exact test for biased strand distribution p-value = 5.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

GCCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCT  >  NC_002947/4807934‑4808231
                                                                                                                                                    |                                                                                                                                                     
gCCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGaa                                                                                                                                                      <  2:5060552/150‑1 (MQ=255)
gCCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGaa                                                                                                                                                      <  1:4458813/150‑1 (MQ=255)
gCCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGaa                                                                                                                                                      >  1:3266468/1‑150 (MQ=255)
   gacgaTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACgg                                                                                                                                                   >  2:1588100/1‑150 (MQ=255)
     cgaTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCCTCGGCCGTTCGGcc                                                                                                                                                 <  1:4256840/150‑1 (MQ=255)
       aTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa                                                                                                                                               <  2:3448601/150‑1 (MQ=255)
       aTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa                                                                                                                                               <  1:131729/150‑1 (MQ=255)
       aTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa                                                                                                                                               <  1:6534614/150‑1 (MQ=255)
       aTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa                                                                                                                                               <  1:186640/150‑1 (MQ=255)
         ggCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCTCCTCGTCGGCCGAACGGCCGAgg                                                                                                                                             <  2:9057517/150‑1 (MQ=255)
         ggCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGacg                                                                                                                                             <  2:6180277/150‑3 (MQ=255)
          gCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacga                                                                                                                                            <  1:297175/150‑1 (MQ=255)
           cGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgag                                                                                                                                           <  2:2300568/150‑2 (MQ=255)
            gACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGaaa                                                                                                                                          >  1:282081/1‑150 (MQ=255)
                 cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGGCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCt                                                                                                                                     <  2:8313998/150‑1 (MQ=255)
                 cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCt                                                                                                                                     <  2:555289/150‑1 (MQ=255)
                 cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCt                                                                                                                                     <  2:4588448/150‑1 (MQ=255)
                 cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCt                                                                                                                                     >  2:1863624/1‑150 (MQ=255)
                 cAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCCTCGGCCGTTCGGCcgacgagg                                                                                                                                          >  2:6002954/1‑143 (MQ=255)
                   gCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGa                                                                                                                                   <  1:4648569/150‑1 (MQ=255)
                     ttGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGacc                                                                                                                                 >  2:566069/1‑150 (MQ=255)
                     ttGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGacc                                                                                                                                 <  1:1112331/150‑1 (MQ=255)
                         cGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGAccacca                                                                                                                             >  1:5245683/1‑150 (MQ=255)
                          gACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCaa                                                                                                                            >  1:4701730/1‑150 (MQ=255)
                          gACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCaa                                                                                                                            >  1:5010048/1‑150 (MQ=255)
                             ccAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCCATTGCAGCTCTGGCCATTCGGTGACTGCAGTCGTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGc                                                                                                                         >  2:7142169/1‑150 (MQ=255)
                                gcgcGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTg                                                                                                                      <  1:8155175/150‑1 (MQ=255)
                                 cgcgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                     <  1:5746499/150‑1 (MQ=255)
                                 cgcgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                     <  1:6270624/150‑1 (MQ=255)
                                   cgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGa                                                                                                                   <  1:877759/150‑1 (MQ=255)
                                   cgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGa                                                                                                                   <  1:7604611/150‑1 (MQ=255)
                                         gcgcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGa                                                                                                             <  2:4256089/150‑1 (MQ=255)
                                          cgcAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGAc                                                                                                            <  2:3242565/150‑1 (MQ=255)
                                             aGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATg                                                                                                         >  1:1730224/1‑150 (MQ=255)
                                               cgcgGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGTAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGcc                                                                                                       >  2:4008550/1‑150 (MQ=255)
                                               cgcgGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACagcga                                                                                                       <  1:5418777/150‑5 (MQ=255)
                                               cgcgGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGcc                                                                                                       >  2:7653182/1‑150 (MQ=255)
                                                gcggcgTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCg                                                                                                      <  1:7202157/150‑1 (MQ=255)
                                                gcggcgTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCg                                                                                                      <  2:8705340/150‑1 (MQ=255)
                                                  ggcgTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGt                                                                                                    <  1:7240303/150‑1 (MQ=255)
                                                     gTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGcc                                                                                                 <  2:2716550/150‑1 (MQ=255)
                                                       ccTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCa                                                                                               <  1:6956886/150‑1 (MQ=255)
                                                       ccTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCa                                                                                               <  1:7487768/150‑1 (MQ=255)
                                                         tGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGc                                                                                             >  2:4767747/1‑150 (MQ=255)
                                                           cTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTc                                                                                           >  2:6920858/1‑150 (MQ=255)
                                                            tCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCg                                                                                          >  2:8703274/1‑150 (MQ=255)
                                                             ccGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCgg                                                                                         <  1:5601235/150‑1 (MQ=255)
                                                                gctggctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTTACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGc                                                                                      >  1:2206096/1‑150 (MQ=255)
                                                                   ggctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGa                                                                                   <  2:5231769/150‑1 (MQ=255)
                                                                      ggATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGatat                                                                                <  2:1796281/149‑1 (MQ=255)
                                                                       gATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATg                                                                               >  2:6644928/1‑150 (MQ=255)
                                                                            ggTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCa                                                                          <  2:8682205/150‑1 (MQ=255)
                                                                             gTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAg                                                                         <  2:4341411/150‑1 (MQ=255)
                                                                                cGGGTGTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAACGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCt                                                                      <  1:1490221/150‑1 (MQ=255)
                                                                                  ggTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGa                                                                    >  1:102425/1‑150 (MQ=255)
                                                                                   gTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGAt                                                                   >  1:8402578/1‑150 (MQ=255)
                                                                                    tttCCGCTGGTGCTCAATTGCAGATCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATg                                                                  >  1:6703008/1‑150 (MQ=255)
                                                                                      tCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAATTCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCg                                                                <  1:2566948/150‑1 (MQ=255)
                                                                                        cGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGc                                                              <  1:8344239/150‑1 (MQ=255)
                                                                                           tGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCAt                                                           >  2:1933995/1‑150 (MQ=255)
                                                                                             gTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCg                                                         <  2:8027827/150‑1 (MQ=255)
                                                                                               gCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgc                                                       >  2:2374700/1‑150 (MQ=255)
                                                                                               gCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgc                                                       >  1:5477791/1‑150 (MQ=255)
                                                                                                 tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgcgc                                                     >  2:2214128/1‑150 (MQ=255)
                                                                                                 tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgcgc                                                     >  1:333409/1‑150 (MQ=255)
                                                                                                 tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgcgc                                                     >  1:3531742/1‑150 (MQ=255)
                                                                                                 tCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGgcgc                                                     >  1:6818058/1‑150 (MQ=255)
                                                                                                  cAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCt                                                    >  1:1090980/1‑150 (MQ=255)
                                                                                                   aaTTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTg                                                   >  2:806507/1‑150 (MQ=255)
                                                                                                     ttGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGaa                                                 >  2:7288401/1‑150 (MQ=255)
                                                                                                       gCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACt                                               >  2:4937248/1‑150 (MQ=255)
                                                                                                        cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              <  2:7836759/150‑1 (MQ=255)
                                                                                                        cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              >  2:3026416/1‑150 (MQ=255)
                                                                                                        cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              <  1:6041192/150‑1 (MQ=255)
                                                                                                        cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              <  1:787641/150‑1 (MQ=255)
                                                                                                        cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              <  2:5860367/150‑1 (MQ=255)
                                                                                                        cAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTc                                              <  2:1091311/150‑1 (MQ=255)
                                                                                                                gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTgg                                      >  2:8861393/1‑150 (MQ=255)
                                                                                                                ccccttCGGCGACCGCAGTCACCACCACGTCGGCCGCCCGGCCGAGGAAATCGCTGACCACCAACGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTgg                                      <  1:7368681/146‑1 (MQ=255)
                                                                                                                   aTTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCa                                   >  2:1724768/1‑150 (MQ=255)
                                                                                                                        gTGACTGCAGTCTTCACCTCGGCGGCCGAACGGCCGAGGAAATCGCTGACCACCAACGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGgc                              <  2:5349369/150‑1 (MQ=255)
                                                                                                                             tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTc                         >  2:4092764/1‑150 (MQ=255)
                                                                                                                             tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTc                         >  2:6347942/1‑150 (MQ=255)
                                                                                                                                  tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAc                    >  1:7790426/1‑150 (MQ=255)
                                                                                                                                   cTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAcc                   >  2:4883635/1‑150 (MQ=255)
                                                                                                                                    ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGGCATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt                  >  2:6413539/1‑150 (MQ=255)
                                                                                                                                    ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt                  >  2:8467023/1‑150 (MQ=255)
                                                                                                                                    ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt                  >  1:2024031/1‑150 (MQ=255)
                                                                                                                                    ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt                  >  1:2403919/1‑150 (MQ=255)
                                                                                                                                    ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt                  <  2:1443357/150‑1 (MQ=255)
                                                                                                                                      cACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGt                >  1:4924104/1‑150 (MQ=255)
                                                                                                                                        ccTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTAc              >  1:8619792/1‑150 (MQ=255)
                                                                                                                                        ccTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTAc              <  1:9004173/150‑1 (MQ=255)
                                                                                                                                            gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa          <  1:3962379/150‑1 (MQ=255)
                                                                                                                                            gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa          <  2:365519/150‑1 (MQ=255)
                                                                                                                                            gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa          <  2:282081/150‑1 (MQ=255)
                                                                                                                                            ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa          >  1:1437535/2‑150 (MQ=255)
                                                                                                                                            ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa          >  2:3329600/2‑150 (MQ=255)
                                                                                                                                                gCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTcc      >  2:7765354/1‑150 (MQ=255)
                                                                                                                                                    aaCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCt  <  2:2920575/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTCCTTCT  >  NC_002947/4807934‑4808231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: