Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,586,676 G→T 7.4% H152N (CAC→AAC)  ypfH ← palmitoyl‑CoA esterase activity, uncertain physiological substrate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,586,6760GT 7.4% 1166.4 / 42.9 304H152N (CAC→AAC) ypfHpalmitoyl‑CoA esterase activity, uncertain physiological substrate
Reads supporting (aligned to +/- strand):  ref base G (138/143);  new base T (11/12);  total (149/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

GGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAGCGCCGTGGCGTTAG  >  CP009273/2586527‑2586825
                                                                                                                                                     |                                                                                                                                                     
 gCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCGCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCAGCAATCAGGTGa                                                                                                                                                      >  1:2748563/1‑150 (MQ=255)
   aaaTTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGact                                                                                                                                                    >  2:3394754/1‑148 (MQ=255)
    aaTTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAAtt                                                                                                                                                   <  2:433794/150‑1 (MQ=255)
    aaTTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAAtt                                                                                                                                                   <  2:2905374/150‑1 (MQ=255)
    aaTTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAAtt                                                                                                                                                   <  2:3855969/150‑1 (MQ=255)
       tGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTc                                                                                                                                                <  2:1062812/150‑1 (MQ=255)
           tACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTgg                                                                                                                                            >  2:2297450/1‑150 (MQ=255)
             cTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCt                                                                                                                                          >  1:17286/1‑150 (MQ=255)
              tGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTg                                                                                                                                         <  2:760571/150‑1 (MQ=255)
               gCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGt                                                                                                                                        <  1:4484457/150‑1 (MQ=255)
               gCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGt                                                                                                                                        >  2:1468349/1‑150 (MQ=255)
                cGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTa                                                                                                                                       >  1:4193794/1‑150 (MQ=255)
                 gATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAg                                                                                                                                      >  2:2128819/1‑150 (MQ=255)
                 gATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAg                                                                                                                                      >  1:1478659/1‑150 (MQ=255)
                 gATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAg                                                                                                                                      >  1:2862484/1‑150 (MQ=255)
                        gATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGt                                                                                                                               >  1:4663541/1‑150 (MQ=255)
                        gATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGt                                                                                                                               >  2:326792/1‑150 (MQ=255)
                        gATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGt                                                                                                                               <  1:2257455/150‑1 (MQ=255)
                            gCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTcc                                                                                                                           >  2:2346412/1‑150 (MQ=255)
                            gCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTcc                                                                                                                           <  2:1407307/150‑1 (MQ=255)
                                ggCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCa                                                                                                                       <  1:3480879/150‑1 (MQ=255)
                                ggCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCa                                                                                                                       <  1:4028884/150‑1 (MQ=255)
                                  ccTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAgg                                                                                                                     <  1:4111223/150‑1 (MQ=255)
                                  ccTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATAAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAAACATGTGAATAGTCGTGGCTGTAGACGCAGTTTCCGGCAgg                                                                                                                     >  1:804013/1‑150 (MQ=255)
                                    tAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCTCTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCt                                                                                                                   <  1:4878095/150‑1 (MQ=255)
                                     aaaTCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTc                                                                                                                  <  2:2441943/150‑1 (MQ=255)
                                     aaaTCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTc                                                                                                                  <  1:2346580/150‑1 (MQ=255)
                                     aaaTCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTc                                                                                                                  <  1:4153117/150‑1 (MQ=255)
                                     aaaTCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTc                                                                                                                  <  1:2624953/150‑1 (MQ=255)
                                      aaTCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCg                                                                                                                 >  1:1201084/1‑150 (MQ=255)
                                       aTCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTTAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGc                                                                                                                <  2:1052890/150‑1 (MQ=255)
                                        tcctccACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCa                                                                                                               <  2:128175/150‑1 (MQ=255)
                                           tccACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCGGTAGACGCAGTTTCCGGCAGGCTCGCATAg                                                                                                            <  1:1538816/150‑1 (MQ=255)
                                               cGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTTAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTc                                                                                                        <  1:4715505/150‑1 (MQ=255)
                                                gATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTcc                                                                                                       >  2:629433/1‑150 (MQ=255)
                                                 atatCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCg                                                                                                      >  2:79168/1‑150 (MQ=255)
                                                      cAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTaaa                                                                                                 <  1:3747034/150‑1 (MQ=255)
                                                       aGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTaaaa                                                                                                >  1:1787934/1‑150 (MQ=255)
                                                         cgtcacgtcaCCACCGGCACTGATTAACGCTTCCTGCTCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTAACCGCCATGAATCAGGTTAATTGTCGTGGCTGTATACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGc                                                                                              >  1:3822117/1‑150 (MQ=255)
                                                          gtcacgtcacCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCa                                                                                             >  1:2966538/1‑150 (MQ=255)
                                                                tcacCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGAcg                                                                                       <  2:2963669/150‑1 (MQ=255)
                                                                   ccaccGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGt                                                                                    <  1:851030/150‑1 (MQ=255)
                                                                      ccGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTTAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGaa                                                                                 <  2:2079220/150‑1 (MQ=255)
                                                                          cACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCaa                                                                             >  1:5121293/1‑150 (MQ=255)
                                                                           aCTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAg                                                                            <  2:2517280/150‑1 (MQ=255)
                                                                            cTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAgg                                                                           <  2:2698508/150‑1 (MQ=255)
                                                                            cTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAgg                                                                           >  1:3213025/1‑150 (MQ=255)
                                                                            cTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGAGTGAAGCAAgg                                                                           >  1:1682233/1‑150 (MQ=255)
                                                                               aTTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGccc                                                                        >  1:1038097/1‑150 (MQ=255)
                                                                                ttAACGCGTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCg                                                                       <  1:1936581/150‑1 (MQ=255)
                                                                                     gCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTTAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCg                                                                  >  1:3942542/1‑150 (MQ=255)
                                                                                      cTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGGCCGGTTCgg                                                                 >  2:1641701/1‑150 (MQ=255)
                                                                                        tCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCt                                                               <  1:2529534/150‑1 (MQ=255)
                                                                                         ccTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCtt                                                              <  1:3027470/150‑1 (MQ=255)
                                                                                          cTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCttt                                                             <  2:2751925/150‑1 (MQ=255)
                                                                                           tGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTa                                                            >  2:4054256/1‑150 (MQ=255)
                                                                                                 gcAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTTAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGctc                                                      >  1:3566364/1‑150 (MQ=255)
                                                                                                  cAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGctct                                                     >  1:824462/1‑150 (MQ=255)
                                                                                                  cAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGctct                                                     >  1:3400042/1‑150 (MQ=255)
                                                                                                   aaCCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCtctc                                                    >  1:5063021/1‑150 (MQ=255)
                                                                                                   aaCCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCtctc                                                    <  1:929861/150‑1 (MQ=255)
                                                                                                     ccGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCa                                                  <  1:3877678/150‑1 (MQ=255)
                                                                                                     ccGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCa                                                  <  1:1548665/150‑1 (MQ=255)
                                                                                                      cGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAg                                                 <  2:4830809/150‑1 (MQ=255)
                                                                                                       gCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTTAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGc                                                <  1:4228254/150‑1 (MQ=255)
                                                                                                       gCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCAGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGc                                                >  1:3354703/1‑150 (MQ=255)
                                                                                                           gcgcCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCa                                            <  2:3061926/150‑1 (MQ=255)
                                                                                                            cgcCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCAt                                           >  2:4727140/1‑150 (MQ=255)
                                                                                                            cgcCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCAt                                           >  1:1092627/1‑150 (MQ=255)
                                                                                                             gcCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa                                          >  2:1568890/1‑150 (MQ=255)
                                                                                                             gcCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATa                                          <  1:500441/150‑1 (MQ=255)
                                                                                                               cAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAAt                                        <  2:902102/150‑1 (MQ=255)
                                                                                                               cAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAAt                                        <  1:4517989/150‑1 (MQ=255)
                                                                                                                aGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATc                                       >  2:1812264/1‑150 (MQ=255)
                                                                                                                aGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATc                                       >  1:1760567/1‑150 (MQ=255)
                                                                                                                aGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATc                                       <  1:4412102/150‑1 (MQ=255)
                                                                                                                     gATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTcc                                  >  1:3095982/1‑150 (MQ=255)
                                                                                                                          ccGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTgag                             <  1:2719290/150‑1 (MQ=255)
                                                                                                                           cGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTCCTTgaga                            <  2:3445657/150‑1 (MQ=255)
                                                                                                                            ggATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTgagag                           >  2:4458735/1‑150 (MQ=255)
                                                                                                                            ggATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTGgagag                           >  2:2980395/1‑150 (MQ=255)
                                                                                                                             gATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTgagaga                          >  1:3299490/1‑150 (MQ=255)
                                                                                                                                cgtCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAAc                       <  2:3262750/148‑1 (MQ=255)
                                                                                                                                cTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAAc                       <  1:3140251/150‑1 (MQ=255)
                                                                                                                                 ttCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTCCTTGAGAGAAAcc                      <  1:599371/150‑1 (MQ=255)
                                                                                                                                 ttCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAAcc                      >  1:4849519/1‑150 (MQ=255)
                                                                                                                                 ttCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAAcc                      <  2:4189193/150‑1 (MQ=255)
                                                                                                                                  tCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTCCTTGAGAGAAACCg                     >  1:1397487/1‑150 (MQ=255)
                                                                                                                                    aCCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTCCTTGAGAGAAACCGAt                   <  1:4280123/150‑1 (MQ=255)
                                                                                                                                    aCCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGAt                   >  2:4884399/1‑150 (MQ=255)
                                                                                                                                     ccGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTCCTTGAGAGAAACCGATg                  >  2:3357207/1‑150 (MQ=255)
                                                                                                                                     ccGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTCCTTGAGAGAAACCGATg                  >  2:5141879/1‑150 (MQ=255)
                                                                                                                                       gCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAg                >  2:3166782/1‑150 (MQ=255)
                                                                                                                                          ttGATTCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAgcgc             <  2:3405128/149‑1 (MQ=255)
                                                                                                                                          aTGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAgcgc             >  1:4973523/1‑150 (MQ=255)
                                                                                                                                              aTCAGGTTAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAGCGCCGTg         >  2:4919227/1‑150 (MQ=255)
                                                                                                                                                  ggTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAGCGCCGTGGCGt     <  1:4070682/150‑1 (MQ=255)
                                                                                                                                                   gTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAGCGCCGTGGCGtt    <  2:4643788/150‑1 (MQ=255)
                                                                                                                                                    tGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCTCCTTGAGAGAAACCGATGAGCGCCGTGGCGTTa   >  1:4608874/1‑150 (MQ=255)
                                                                                                                                                     gAAGTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAGCGCCGTGGCGTTAg  <  2:3607999/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCAAATTGCATACTGCGATTGTCGATCGCATGGCCTAAATCCTCCACGATATCCAGCGTCACGTCACCACCGGCACTGATTAACGCTTCCTGCGCCGCAACCGCATGCGCCAGATCGATAACCGGATCTTCACCGCCATGAATCAGGTGAATTGTCGTGGCTGTAGACGCAGTTTCCGGCAGGCTCGCATAGCGTCCGTTAAAAGCAATGACGCGTGAAGCAAGGCCCGGTTCGGCTTTAATGCTCTCCAGCACCATAATCGCGCCTTGAGAGAAACCGATGAGCGCCGTGGCGTTAG  >  CP009273/2586527‑2586825

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: