Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 598,618 G→T 5.3% E402* (GAA→TAA)  pheP → phenylalanine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273598,6180GT 5.3% 1099.9 / 22.1 282E402* (GAA→TAA) phePphenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base G (120/146);  new base T (6/10);  total (126/156)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

CTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCA  >  CP009273/598469‑598767
                                                                                                                                                     |                                                                                                                                                     
    gCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAac                                                                                                                                                   <  1:2819925/150‑1 (MQ=255)
      tggtggTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAacac                                                                                                                                                 >  1:2160515/1‑150 (MQ=255)
       tgtgttgTTAATCAAATATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTATCAACGCTGATTTTGAACTGGATTATGATCTGTCTGGCGCATCTGCTTTTTAGTTCAGCGATTAGACGTCAGGGGCGTTAAAcaca                                                                                                                                                <  2:2431325/149‑1 (MQ=255)
       ggtggtGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAAcaca                                                                                                                                                >  1:2963111/1‑150 (MQ=255)
       ggtggtGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAAcaca                                                                                                                                                <  2:693373/150‑1 (MQ=255)
       ggtggtGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAAcaca                                                                                                                                                >  1:2726171/1‑150 (MQ=255)
        gtggtgTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAg                                                                                                                                               >  2:4382770/1‑150 (MQ=255)
           gtgtTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGttt                                                                                                                                            >  1:2404083/1‑150 (MQ=255)
             gtTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTaa                                                                                                                                          <  1:1694342/150‑1 (MQ=255)
             gtTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTaa                                                                                                                                          <  1:1924324/150‑1 (MQ=255)
              ttAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAg                                                                                                                                         <  1:4639428/150‑1 (MQ=255)
                   cAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCgct                                                                                                                                    <  2:306898/150‑1 (MQ=255)
                   cAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCgct                                                                                                                                    <  2:12884/150‑1 (MQ=255)
                      cTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCgctgct                                                                                                                                 >  2:4407744/1‑150 (MQ=255)
                        aTCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGctct                                                                                                                               <  2:848134/150‑1 (MQ=255)
                        aTCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGctct                                                                                                                               >  2:1115896/1‑150 (MQ=255)
                        aTCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGctct                                                                                                                               <  2:2931646/150‑1 (MQ=255)
                        aTCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGctct                                                                                                                               >  1:3231496/1‑150 (MQ=255)
                        aTCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGAGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAAAACAGTTTAAGGCGCTGctct                                                                                                                               <  2:3865259/150‑1 (MQ=255)
                          ctgctgCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTAt                                                                                                                             >  1:1282760/1‑150 (MQ=255)
                           tgctgcCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATc                                                                                                                            <  1:4240270/150‑1 (MQ=255)
                           tgctgcCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATc                                                                                                                            <  1:4206892/150‑1 (MQ=255)
                             ctgcCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCg                                                                                                                          <  2:4230507/150‑1 (MQ=255)
                                ccGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTc                                                                                                                       >  2:1218750/1‑150 (MQ=255)
                                  gCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCgg                                                                                                                     <  1:1676522/150‑1 (MQ=255)
                                  gCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCgg                                                                                                                     <  2:2073031/150‑1 (MQ=255)
                                     aaaaGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCaa                                                                                                                  >  1:2661777/1‑150 (MQ=255)
                                        aGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTa                                                                                                               >  2:1822083/1‑150 (MQ=255)
                                        aGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTa                                                                                                               >  2:780511/1‑150 (MQ=255)
                                           gTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTAtct                                                                                                            <  1:4103055/150‑1 (MQ=255)
                                             ttGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATctct                                                                                                          <  2:875276/150‑1 (MQ=255)
                                             ttGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATctct                                                                                                          <  2:3145415/150‑1 (MQ=255)
                                               ggTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGc                                                                                                        <  2:3666368/150‑1 (MQ=255)
                                                gTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCa                                                                                                       <  1:3967017/150‑1 (MQ=255)
                                                  ctgctgATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCAtt                                                                                                     >  1:2581396/1‑150 (MQ=255)
                                                     ctgATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGcc                                                                                                  <  1:1459898/150‑1 (MQ=255)
                                                       gATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCtt                                                                                                <  1:318704/150‑1 (MQ=255)
                                                       gATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCtt                                                                                                >  1:2033007/1‑150 (MQ=255)
                                                               tggtggTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCa                                                                                        <  2:3904429/150‑1 (MQ=255)
                                                                 gtggtAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATg                                                                                      <  1:1473078/150‑1 (MQ=255)
                                                                    gtAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGAtt                                                                                   <  2:5749/150‑1 (MQ=255)
                                                                     tAGCATCGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGAttt                                                                                  <  1:1001165/150‑1 (MQ=255)
                                                                      aGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGAtttt                                                                                 <  1:616334/150‑1 (MQ=255)
                                                                         aaCGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTtgct                                                                              <  1:2319618/150‑1 (MQ=255)
                                                                            gctgctGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTgctgct                                                                           >  1:3340390/1‑150 (MQ=255)
                                                                              tgctgTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGa                                                                         >  1:3822901/1‑150 (MQ=255)
                                                                               gctgTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGAt                                                                        <  1:1876155/150‑1 (MQ=255)
                                                                                    tGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCa                                                                   <  2:1542604/150‑1 (MQ=255)
                                                                                    tGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCa                                                                   <  1:2066531/150‑1 (MQ=255)
                                                                                     gAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCAc                                                                  >  1:2198735/1‑150 (MQ=255)
                                                                                      aaCTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACg                                                                 <  2:241340/150‑1 (MQ=255)
                                                                                      aaCTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACg                                                                 >  2:2370871/1‑150 (MQ=255)
                                                                                        cTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACgat                                                               >  2:2627479/1‑150 (MQ=255)
                                                                                        cTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACgat                                                               <  1:3538750/150‑1 (MQ=255)
                                                                                            aTTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACgatggat                                                           >  2:2185375/1‑150 (MQ=255)
                                                                                                 gATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGatat                                                      >  2:2491705/1‑150 (MQ=255)
                                                                                                 gATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGatat                                                      >  2:302394/1‑150 (MQ=255)
                                                                                                    ctgtctgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATgcg                                                   >  2:1334620/1‑150 (MQ=255)
                                                                                                      gtctgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCt                                                 >  1:4559571/1‑150 (MQ=255)
                                                                                                        ctgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTg                                               <  1:3370768/150‑1 (MQ=255)
                                                                                                         tgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGt                                              >  1:763327/1‑150 (MQ=255)
                                                                                                         tgGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGt                                              <  1:1264476/150‑1 (MQ=255)
                                                                                                          ggCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTc                                             <  1:2220752/150‑1 (MQ=255)
                                                                                                          ggCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTc                                             >  1:1890728/1‑150 (MQ=255)
                                                                                                           gcgcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCa                                            <  2:1398973/150‑1 (MQ=255)
                                                                                                           gcgcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCa                                            <  2:4595173/150‑1 (MQ=255)
                                                                                                            cgcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAg                                           <  1:1206048/150‑1 (MQ=255)
                                                                                                            cgcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAg                                           <  2:489627/150‑1 (MQ=255)
                                                                                                             gcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGc                                          >  2:3937718/1‑150 (MQ=255)
                                                                                                             gcATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGc                                          >  2:3745086/1‑150 (MQ=255)
                                                                                                             gcATCTGCGTTTTCGTGCAGCGATGCGACGACAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGc                                          >  1:4290759/1‑150 (MQ=255)
                                                                                                               aTCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGa                                        <  1:934126/150‑1 (MQ=255)
                                                                                                                tCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGAt                                       <  2:3436858/150‑1 (MQ=255)
                                                                                                                   gCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCct                                    <  1:3617573/150‑1 (MQ=255)
                                                                                                                      ttttCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCctgct                                 <  2:1907793/150‑1 (MQ=255)
                                                                                                                         tCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGcc                              <  2:3231496/150‑1 (MQ=255)
                                                                                                                            tGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGgt                           >  1:1974666/1‑150 (MQ=255)
                                                                                                                              cAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCAgtgt                         <  1:4421046/150‑1 (MQ=255)
                                                                                                                              cAGCGATGCGACGTCAGGGGCGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCAgtgt                         <  1:3082353/150‑1 (MQ=255)
                                                                                                                              cAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGgtgt                         >  1:3277795/1‑150 (MQ=255)
                                                                                                                               aGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCGTTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGtgtg                        <  2:4012187/150‑1 (MQ=255)
                                                                                                                               aGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGtgtg                        >  1:1575032/1‑150 (MQ=255)
                                                                                                                                 cGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGa                      <  2:2726171/150‑1 (MQ=255)
                                                                                                                                  gATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGAt                     >  1:4120519/1‑150 (MQ=255)
                                                                                                                                  gATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGAt                     >  1:1160047/1‑150 (MQ=255)
                                                                                                                                   aTGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGAtt                    >  1:2496475/1‑150 (MQ=255)
                                                                                                                                      cGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTa                 <  1:1035767/150‑1 (MQ=255)
                                                                                                                                      cGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTa                 <  1:96907/150‑1 (MQ=255)
                                                                                                                                      cGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTa                 <  1:2453788/150‑1 (MQ=255)
                                                                                                                                        aCGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTAtt               <  2:2586232/150‑1 (MQ=255)
                                                                                                                                          gTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTcc             >  1:3770607/1‑150 (MQ=255)
                                                                                                                                           tCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCt            <  1:2113629/150‑1 (MQ=255)
                                                                                                                                            cAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTg           >  2:1465766/1‑150 (MQ=255)
                                                                                                                                            cAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTg           >  1:4480471/1‑150 (MQ=255)
                                                                                                                                                ggCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTa       <  2:1766291/150‑1 (MQ=255)
                                                                                                                                                  cGTTAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATg     >  1:2861283/1‑150 (MQ=255)
                                                                                                                                                  cGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATg     <  1:2893819/150‑1 (MQ=255)
                                                                                                                                                  cGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATg     <  1:2839175/150‑1 (MQ=255)
                                                                                                                                                   gTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATgg    <  1:560887/150‑1 (MQ=255)
                                                                                                                                                     gAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCa  >  1:498949/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTCGCTGGTGGTGTTAATCAACTATCTGCTGCCGCAAAAAGCGTTTGGTCTGCTGATGGCGCTGGTGGTAGCAACGCTGCTGTTGAACTGGATTATGATCTGTCTGGCGCATCTGCGTTTTCGTGCAGCGATGCGACGTCAGGGGCGTGAAACACAGTTTAAGGCGCTGCTCTATCCGTTCGGCAACTATCTCTGCATTGCCTTCCTCGGCATGATTTTGCTGCTGATGTGCACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCA  >  CP009273/598469‑598767

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: