Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,358,855 G→T 7.1% G123C (GGC→TGC)  puuB → gamma‑glutamylputrescine oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,358,8550GT 7.1% 1141.9 / 36.4 300G123C (GGC→TGC) puuBgamma‑glutamylputrescine oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (133/145);  new base T (12/10);  total (145/155)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 314 total aligned reads displayed.

TGTGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCG  >  CP009273/1358707‑1359004
                                                                                                                                                    |                                                                                                                                                     
tgtgAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCcgg                                                                                                                                                      >  1:3929551/1‑150 (MQ=255)
  tgAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCcggcg                                                                                                                                                    <  1:1381240/150‑1 (MQ=255)
  tgAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCcggcg                                                                                                                                                    <  2:2705783/150‑1 (MQ=255)
    aaCTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCggcggc                                                                                                                                                  >  2:1399096/1‑150 (MQ=255)
    aaCTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCggcggc                                                                                                                                                  >  1:706215/1‑150 (MQ=255)
       tCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTg                                                                                                                                               >  2:2024798/1‑150 (MQ=255)
        ccTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGt                                                                                                                                              >  2:2394945/1‑150 (MQ=255)
         cTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGtt                                                                                                                                             <  2:1436104/150‑1 (MQ=255)
           aCAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCATATTGACTGCGACTACCGCCCCGGCGGCCTGTTtg                                                                                                                                           <  1:3223109/150‑1 (MQ=255)
            cAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTtgt                                                                                                                                          <  1:651189/150‑1 (MQ=255)
                cgcgACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCg                                                                                                                                      >  1:3316507/1‑150 (MQ=255)
                  cgACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGat                                                                                                                                    <  1:27377/150‑1 (MQ=255)
                    aCATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGatga                                                                                                                                  <  2:1691115/150‑1 (MQ=255)
                      aTCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGatgaat                                                                                                                                <  2:4224873/150‑1 (MQ=255)
                         gACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                             <  1:2920381/150‑1 (MQ=255)
                            gTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATaaa                                                                                                                          <  1:19471/150‑1 (MQ=255)
                            gTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATaaa                                                                                                                          <  1:533123/150‑1 (MQ=255)
                              gATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACa                                                                                                                        >  2:2552453/1‑150 (MQ=255)
                              gATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACa                                                                                                                        >  2:3634586/1‑150 (MQ=255)
                               aTCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAg                                                                                                                       >  1:933930/1‑150 (MQ=255)
                                   aaaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCg                                                                                                                   >  2:3795477/1‑150 (MQ=255)
                                   aaaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCg                                                                                                                   >  2:341095/1‑150 (MQ=255)
                                    aaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGc                                                                                                                  >  2:4605421/1‑150 (MQ=255)
                                      aaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCca                                                                                                                >  2:1935153/1‑150 (MQ=255)
                                      aaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCca                                                                                                                <  2:3251030/150‑1 (MQ=255)
                                         gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCacac                                                                                                             <  2:2965271/150‑1 (MQ=255)
                                          cTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                            <  1:4384459/150‑1 (MQ=255)
                                           tACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACtt                                                                                                           >  1:1522090/1‑150 (MQ=255)
                                                   ggACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCa                                                                                                   >  2:2496189/1‑150 (MQ=255)
                                                    gACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAg                                                                                                  <  1:2349650/150‑1 (MQ=255)
                                                      cacCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaa                                                                                                <  2:2194501/150‑1 (MQ=255)
                                                       acCGCCCGTATGCTCGGCAGCATGATGTTCGACGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaaa                                                                                               <  2:4294825/150‑1 (MQ=255)
                                                       acCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaaa                                                                                               >  1:2782324/1‑150 (MQ=255)
                                                        ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAag                                                                                              <  2:2944405/150‑1 (MQ=255)
                                                           cccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                           <  2:600994/150‑1 (MQ=255)
                                                            ccGTATGCTCGTCAGAATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGaa                                                                                          <  2:4018938/150‑1 (MQ=255)
                                                            ccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGaa                                                                                          <  1:4130750/150‑1 (MQ=255)
                                                            ccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGaa                                                                                          <  2:1111987/150‑1 (MQ=255)
                                                            ccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGaa                                                                                          <  2:461988/150‑1 (MQ=255)
                                                             cGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAAc                                                                                         >  2:550377/1‑150 (MQ=255)
                                                                 tGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTggg                                                                                     <  1:1165407/150‑1 (MQ=255)
                                                                     cGGCAGCATGTTGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACg                                                                                 <  1:349719/150‑1 (MQ=255)
                                                                     cGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACg                                                                                 <  1:411076/150‑1 (MQ=255)
                                                                      ggCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGc                                                                                >  2:3951359/1‑150 (MQ=255)
                                                                       gcagcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCt                                                                               <  2:3807904/150‑1 (MQ=255)
                                                                         agcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTAc                                                                             <  1:1008083/150‑1 (MQ=255)
                                                                           caTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACgg                                                                           >  2:1425584/1‑150 (MQ=255)
                                                                            atgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGc                                                                          >  2:3543462/1‑150 (MQ=255)
                                                                            atgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGc                                                                          >  1:3787104/1‑150 (MQ=255)
                                                                               atgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAAt                                                                       <  2:712538/150‑1 (MQ=255)
                                                                               atgTTCGAAGGTGGTGAGATCATCCGCGAACGCATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAAt                                                                       <  2:3883583/150‑1 (MQ=255)
                                                                                tgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATa                                                                      <  1:3541469/150‑1 (MQ=255)
                                                                                  ttCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATaaa                                                                    >  2:3354381/1‑150 (MQ=255)
                                                                                   tCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAAc                                                                   <  1:4268404/150‑1 (MQ=255)
                                                                                    cGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACa                                                                  >  1:4464001/1‑150 (MQ=255)
                                                                                     gAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAg                                                                 >  2:4189909/1‑150 (MQ=255)
                                                                                     gAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAg                                                                 <  1:419040/150‑1 (MQ=255)
                                                                                        ggtggtGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTg                                                              >  1:3286909/1‑150 (MQ=255)
                                                                                        ggtggtGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTg                                                              >  1:1415109/1‑150 (MQ=255)
                                                                                         gtggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTgg                                                             >  1:1706661/1‑150 (MQ=255)
                                                                                         gtggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTgg                                                             <  1:1894923/150‑1 (MQ=255)
                                                                                          tggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGa                                                            >  1:1864210/1‑150 (MQ=255)
                                                                                            gtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAt                                                          <  2:4369565/150‑1 (MQ=255)
                                                                                             tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                         <  1:917584/150‑1 (MQ=255)
                                                                                             tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                         <  1:2165301/150‑1 (MQ=255)
                                                                                               agaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGc                                                       >  2:2900721/1‑150 (MQ=255)
                                                                                                 atcatcCGCGAACGTATCTAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTg                                                     <  2:2157076/150‑1 (MQ=255)
                                                                                                 atcatcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTg                                                     <  1:17993/150‑1 (MQ=255)
                                                                                                   catcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGa                                                   >  1:2263133/1‑150 (MQ=255)
                                                                                                   catcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGa                                                   >  2:2327901/1‑150 (MQ=255)
                                                                                                    atcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGac                                                  >  2:3708031/1‑150 (MQ=255)
                                                                                                             cGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacgcc                                         <  1:1941704/150‑1 (MQ=255)
                                                                                                               tATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAt                                       <  2:4113051/150‑1 (MQ=255)
                                                                                                               tATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAt                                       >  2:2978238/1‑150 (MQ=255)
                                                                                                               tATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAt                                       <  1:3276692/150‑1 (MQ=255)
                                                                                                                aTCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAtt                                      <  2:760477/150‑1 (MQ=255)
                                                                                                                aTCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAtt                                      >  2:3165961/1‑150 (MQ=255)
                                                                                                                 tCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTc                                     <  2:2122275/150‑1 (MQ=255)
                                                                                                                 tCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTc                                     >  2:3992235/1‑150 (MQ=255)
                                                                                                                 tCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTc                                     <  2:3918993/150‑1 (MQ=255)
                                                                                                                   aaaCGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGc                                   >  2:4472874/1‑150 (MQ=255)
                                                                                                                    aaCGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGcc                                  <  1:2958444/150‑1 (MQ=255)
                                                                                                                        ttATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGa                              <  2:1287083/150‑1 (MQ=255)
                                                                                                                         tATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGaa                             <  1:2492483/150‑1 (MQ=255)
                                                                                                                         tATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGaa                             >  1:3690048/1‑150 (MQ=255)
                                                                                                                          aTCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                            >  2:2825365/1‑150 (MQ=255)
                                                                                                                          aTCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                            >  2:3137439/1‑150 (MQ=255)
                                                                                                                             aGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAg                         <  2:1888858/150‑1 (MQ=255)
                                                                                                                              gATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGc                        <  1:2647781/150‑1 (MQ=255)
                                                                                                                               aTTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGcc                       >  1:1092815/1‑150 (MQ=255)
                                                                                                                               aTTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGcc                       >  1:3743774/1‑150 (MQ=255)
                                                                                                                                 tGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAg                     <  1:3112836/150‑1 (MQ=255)
                                                                                                                                        gACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGc              >  1:4351553/1‑150 (MQ=255)
                                                                                                                                           tACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtat           >  1:556263/1‑150 (MQ=255)
                                                                                                                                           tACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtat           <  2:187840/150‑1 (MQ=255)
                                                                                                                                             ccGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAc         <  2:1528458/150‑1 (MQ=255)
                                                                                                                                                 cccGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGt     >  1:1018146/1‑150 (MQ=255)
                                                                                                                                                  ccGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTAGACCGGTg    <  2:525088/150‑1 (MQ=255)
                                                                                                                                                   cggcggCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTgc   >  1:3164514/1‑150 (MQ=255)
                                                                                                                                                   cggcggCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTgc   <  1:747180/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGCG  >  CP009273/1358707‑1359004

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: