Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,358,855 G→T 100% 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,8550GT95.6% 602.9 / 17.7 207G123C (GGC→TGC) puuBgamma‑glutamylputrescine oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base T (86/112);  total (90/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 217 total aligned reads displayed.

GAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAC  >  CP009273/1358710‑1358997
                                                                                                                                                 |                                                                                                                                              
gAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCctgc                                                                                                                                              >  1:2461188/1‑148 (MQ=255)
  aCTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCctgcgg                                                                                                                                            <  1:3160096/148‑1 (MQ=255)
  aCTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCctgcgg                                                                                                                                            >  2:2326872/1‑148 (MQ=255)
   cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCtgcggc                                                                                                                                           <  1:3302398/148‑1 (MQ=255)
   cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCtgcggc                                                                                                                                           <  1:3379087/148‑1 (MQ=255)
   cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCtgcggc                                                                                                                                           <  2:3133841/148‑1 (MQ=255)
       tACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGt                                                                                                                                       <  2:177804/148‑1 (MQ=255)
        aCAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGtt                                                                                                                                      <  1:1025219/148‑1 (MQ=255)
         cAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGttt                                                                                                                                     <  1:400822/148‑1 (MQ=255)
         cAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGttt                                                                                                                                     <  2:2825926/148‑1 (MQ=255)
             cgcgACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTgg                                                                                                                                 >  1:3528990/1‑148 (MQ=255)
               cgACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCTAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCg                                                                                                                               >  1:2170787/1‑148 (MQ=255)
                gACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGa                                                                                                                              <  2:418164/148‑1 (MQ=255)
                 aCATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGat                                                                                                                             >  1:1653945/1‑148 (MQ=255)
                 aCATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGat                                                                                                                             >  2:3102353/1‑148 (MQ=255)
                   aTCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatga                                                                                                                           >  2:1374345/1‑148 (MQ=255)
                      gACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatgaatg                                                                                                                        <  2:179575/148‑1 (MQ=255)
                        cGTGATCTAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                      <  1:17573/148‑1 (MQ=255)
                        cGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                      <  1:1961408/148‑1 (MQ=255)
                          tGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                    <  2:1314946/148‑1 (MQ=255)
                          tGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGAGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                    <  2:2946489/148‑1 (MQ=255)
                          tGATCGAAAAAAGCTACGGCATGGACACCGCCCGGATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGTACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                    <  2:3360390/148‑1 (MQ=255)
                           gATCGAGAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaaa                                                                                                                   <  2:3180702/148‑1 (MQ=255)
                           gATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaaa                                                                                                                   <  2:3180210/148‑1 (MQ=255)
                           gATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaaa                                                                                                                   <  2:1503411/148‑1 (MQ=255)
                           gATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTGTGTGGCGATGAATGATaaa                                                                                                                   >  2:3277231/1‑148 (MQ=255)
                            aTCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAAc                                                                                                                  <  2:191259/148‑1 (MQ=255)
                                   aaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGc                                                                                                           >  1:1166829/1‑148 (MQ=255)
                                    aaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGcc                                                                                                          <  2:1285883/148‑1 (MQ=255)
                                      gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                        >  1:3060110/1‑148 (MQ=255)
                                       cTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcaca                                                                                                       >  2:2053188/1‑148 (MQ=255)
                                        tACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCacac                                                                                                      <  1:1404051/148‑1 (MQ=255)
                                         aCGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                     >  2:1718233/1‑148 (MQ=255)
                                             cATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGaa                                                                                                 >  2:1864616/1‑148 (MQ=255)
                                             cATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGaa                                                                                                 >  1:750014/1‑148 (MQ=255)
                                              aTGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAag                                                                                                <  1:2517529/148‑1 (MQ=255)
                                                ggACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAagag                                                                                              <  2:3069195/148‑1 (MQ=255)
                                                ggACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAagag                                                                                              >  2:762860/1‑148 (MQ=255)
                                                 gACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGc                                                                                             <  2:3104812/148‑1 (MQ=255)
                                                 gACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGc                                                                                             >  2:32233/1‑148 (MQ=255)
                                                 gACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGc                                                                                             <  1:1718233/148‑1 (MQ=255)
                                                   cacCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAg                                                                                           >  2:1992929/1‑148 (MQ=255)
                                                    acCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGa                                                                                          <  2:928394/148‑1 (MQ=255)
                                                     ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaa                                                                                         <  2:1145174/148‑1 (MQ=255)
                                                       gCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAag                                                                                       <  1:929528/148‑1 (MQ=255)
                                                        cccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                      <  1:1462987/148‑1 (MQ=255)
                                                         ccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAagag                                                                                     <  1:958371/148‑1 (MQ=255)
                                                         ccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAagag                                                                                     <  1:2059887/148‑1 (MQ=255)
                                                          cGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                    <  1:1866667/148‑1 (MQ=255)
                                                          cGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                    <  1:509593/148‑1 (MQ=255)
                                                          cGTATGCTCGGCAGCAAGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                    <  2:2524880/148‑1 (MQ=255)
                                                            tATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAAc                                                                                  >  2:3048049/1‑148 (MQ=255)
                                                             aTGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                 <  1:2173964/148‑1 (MQ=255)
                                                             aTGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                 <  1:2173874/148‑1 (MQ=255)
                                                             aTGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                 >  1:852700/1‑148 (MQ=255)
                                                              tGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTg                                                                                <  2:852700/148‑1 (MQ=255)
                                                               gCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTgg                                                                               <  1:784827/148‑1 (MQ=255)
                                                                 tCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGa                                                                             <  2:2404084/148‑1 (MQ=255)
                                                                   ggCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAAc                                                                           <  1:2671759/148‑1 (MQ=255)
                                                                            atgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCa                                                                  >  1:105671/1‑148 (MQ=255)
                                                                            atgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCa                                                                  >  1:30743/1‑148 (MQ=255)
                                                                                   aaGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACa                                                           >  2:753043/1‑148 (MQ=255)
                                                                                     ggtggtGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGc                                                         >  2:2066563/1‑148 (MQ=255)
                                                                                     ggtggtGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGc                                                         <  1:3223716/148‑1 (MQ=255)
                                                                                      gtggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCt                                                        <  2:46883/148‑1 (MQ=255)
                                                                                          tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGaa                                                    <  2:360804/148‑1 (MQ=255)
                                                                                          tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGaa                                                    <  1:2811971/148‑1 (MQ=255)
                                                                                          tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGaa                                                    <  2:1629016/148‑1 (MQ=255)
                                                                                              atcatcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGc                                                <  2:3178284/148‑1 (MQ=255)
                                                                                                 atcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTgg                                             <  2:3539635/148‑1 (MQ=255)
                                                                                                  tcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGa                                            >  1:3436042/1‑148 (MQ=255)
                                                                                                  tcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGa                                            >  2:1760874/1‑148 (MQ=255)
                                                                                                  tcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGa                                            <  1:1302485/148‑1 (MQ=255)
                                                                                                  tcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACATGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGATGCAATTGCTGGa                                            >  1:3289228/1‑148 (MQ=255)
                                                                                                    cgcgAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacg                                          <  2:1035308/148‑1 (MQ=255)
                                                                                                    cgcgAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacg                                          <  2:1478260/148‑1 (MQ=255)
                                                                                                      cgAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgcc                                        >  1:3480215/1‑148 (MQ=255)
                                                                                                       gAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgcca                                       <  2:653258/148‑1 (MQ=255)
                                                                                                       gAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgcca                                       <  1:2814729/148‑1 (MQ=255)
                                                                                                       gAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgcca                                       <  1:3546866/148‑1 (MQ=255)
                                                                                                           ttATCAAAAGTTATCTGATTGACTGCGTCTACCGCCCCTGCGGCATTTTTGTGGAGATGAATTATAAACATCTCGACACACTTGAATAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacgc                                   <  1:659589/147‑1 (MQ=255)
                                                                                                             aTCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacgcca                                 >  1:2890557/1‑148 (MQ=255)
                                                                                                              tCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAt                                >  2:1042748/1‑148 (MQ=255)
                                                                                                              tCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGCCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAt                                >  2:2822301/1‑148 (MQ=255)
                                                                                                                  aCGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGc                            >  2:1787105/1‑148 (MQ=255)
                                                                                                                      tATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTg                        >  2:2083533/1‑148 (MQ=255)
                                                                                                                        tCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGACTTGCTGGACGCCAACGCCATTCGCCGTGaa                      >  2:1544790/1‑148 (MQ=255)
                                                                                                                        tCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGaa                      >  2:1778262/1‑148 (MQ=255)
                                                                                                                         cAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                     <  2:2428365/148‑1 (MQ=255)
                                                                                                                         cAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                     <  1:2645382/148‑1 (MQ=255)
                                                                                                                         cAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                     <  2:439938/148‑1 (MQ=255)
                                                                                                                         cAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                     >  2:2656103/1‑148 (MQ=255)
                                                                                                                         cAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAg                     >  2:906184/1‑148 (MQ=255)
                                                                                                                           gATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTa                   >  2:1723499/1‑148 (MQ=255)
                                                                                                                               gACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCa               >  1:2882038/1‑148 (MQ=255)
                                                                                                                                   gCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGa           <  2:230468/148‑1 (MQ=255)
                                                                                                                                      aCTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCg        >  1:2357860/1‑148 (MQ=255)
                                                                                                                                          ccGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtat    >  2:2633200/1‑148 (MQ=255)
                                                                                                                                           cGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtata   <  1:753043/148‑1 (MQ=255)
                                                                                                                                            gCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAACCAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAc  <  2:1153172/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAC  >  CP009273/1358710‑1358997

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: