Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,358,855 G→T 86.7% 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,8550GT86.7% 718.1 / 114.9 299G123C (GGC→TGC) puuBgamma‑glutamylputrescine oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (16/24);  new base T (140/119);  total (156/143)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

TGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGC  >  CP009273/1358709‑1359003
                                                                                                                                                  |                                                                                                                                                    
    cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGcc                                                                                                                                                <  1:1450026/149‑1 (MQ=255)
    cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGcc                                                                                                                                                <  2:1007254/149‑1 (MQ=255)
    cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGcc                                                                                                                                                <  1:406650/149‑1 (MQ=255)
       cTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGt                                                                                                                                             <  2:3727214/149‑1 (MQ=255)
        tACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGtt                                                                                                                                            >  1:3709265/1‑149 (MQ=255)
        tACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGtt                                                                                                                                            >  2:3005625/1‑149 (MQ=255)
         aCAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGttt                                                                                                                                           <  2:3335751/149‑1 (MQ=255)
           aGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGTGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTtgt                                                                                                                                         <  2:812191/149‑1 (MQ=255)
            gCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTtgtg                                                                                                                                        >  2:3842274/1‑149 (MQ=255)
               gcgACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCg                                                                                                                                     >  1:2812470/1‑149 (MQ=255)
                cgACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGa                                                                                                                                    <  1:105700/149‑1 (MQ=255)
                 gACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGat                                                                                                                                   >  1:1727455/1‑149 (MQ=255)
                    aTCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatgaa                                                                                                                                <  2:4526032/149‑1 (MQ=255)
                     tCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatgaat                                                                                                                               >  2:3242654/1‑149 (MQ=255)
                       gACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGatgaatga                                                                                                                             >  2:4333464/1‑149 (MQ=255)
                        aCGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                            <  2:1599765/149‑1 (MQ=255)
                        aCGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                            <  2:1350305/149‑1 (MQ=255)
                        aCGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                            <  1:3239536/149‑1 (MQ=255)
                        aCGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACATATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                            >  2:2692642/1‑149 (MQ=255)
                          gTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                          >  1:2195058/1‑149 (MQ=255)
                                gAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCt                                                                                                                    <  2:3437587/149‑1 (MQ=255)
                                gAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCt                                                                                                                    <  1:3872276/149‑1 (MQ=255)
                                 aaaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTc                                                                                                                   >  1:3135753/1‑149 (MQ=255)
                                 aaaaaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTc                                                                                                                   >  2:1161887/1‑149 (MQ=255)
                                      aGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                              <  2:2893479/149‑1 (MQ=255)
                                        cTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCacac                                                                                                            <  2:4739900/149‑1 (MQ=255)
                                         tACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                           >  2:622447/1‑149 (MQ=255)
                                         tACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                           <  2:1475175/149‑1 (MQ=255)
                                         tACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                           <  2:1858019/149‑1 (MQ=255)
                                          aCGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACtt                                                                                                          >  1:3959639/1‑149 (MQ=255)
                                            ggCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGa                                                                                                        <  1:3802081/149‑1 (MQ=255)
                                              cATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAag                                                                                                      <  1:69350/149‑1 (MQ=255)
                                               aTGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAaga                                                                                                     >  2:4450582/1‑149 (MQ=255)
                                                    cacCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGa                                                                                                >  1:3621949/1‑149 (MQ=255)
                                                     acCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaa                                                                                               <  1:3255009/149‑1 (MQ=255)
                                                     acCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaa                                                                                               >  2:4184727/1‑149 (MQ=255)
                                                      ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaaa                                                                                              <  1:3308616/149‑1 (MQ=255)
                                                      ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaaa                                                                                              <  1:3064190/149‑1 (MQ=255)
                                                       cGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAag                                                                                             <  2:758971/149‑1 (MQ=255)
                                                        gCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                            <  1:1290046/149‑1 (MQ=255)
                                                        gCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                            <  1:4709678/149‑1 (MQ=255)
                                                        gCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                            >  2:1320208/1‑149 (MQ=255)
                                                        gCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                            >  1:3294161/1‑149 (MQ=255)
                                                            gTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAAc                                                                                        <  2:1935548/149‑1 (MQ=255)
                                                             tATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                       >  2:177141/1‑149 (MQ=255)
                                                             tATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                       <  1:1835379/149‑1 (MQ=255)
                                                             tATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                       <  1:537380/149‑1 (MQ=255)
                                                              aTGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTg                                                                                      >  1:4349874/1‑149 (MQ=255)
                                                                 cTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGa                                                                                   <  1:1867999/149‑1 (MQ=255)
                                                                  tCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGaa                                                                                  <  2:1404646/149‑1 (MQ=255)
                                                                  tCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGaa                                                                                  <  1:508044/149‑1 (MQ=255)
                                                                     gcagcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGc                                                                               >  2:247025/1‑149 (MQ=255)
                                                                     gcagcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGc                                                                               <  1:575853/149‑1 (MQ=255)
                                                                        gcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTAc                                                                            >  1:2337231/1‑149 (MQ=255)
                                                                         caTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACg                                                                           >  2:3391379/1‑149 (MQ=255)
                                                                           tgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGc                                                                         >  1:2365090/1‑149 (MQ=255)
                                                                               gTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATa                                                                     >  2:3750031/1‑149 (MQ=255)
                                                                               gTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATa                                                                     >  1:4615394/1‑149 (MQ=255)
                                                                                 tCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATaaa                                                                   <  1:3022006/149‑1 (MQ=255)
                                                                                   gAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACa                                                                 >  2:4716723/1‑149 (MQ=255)
                                                                                   gAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACa                                                                 >  2:2647980/1‑149 (MQ=255)
                                                                                    aaGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAg                                                                <  2:4387031/149‑1 (MQ=255)
                                                                                      ggtggtGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCt                                                              <  1:4381784/149‑1 (MQ=255)
                                                                                       gtggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTg                                                             >  2:3322541/1‑149 (MQ=255)
                                                                                       gtggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTg                                                             <  1:205137/149‑1 (MQ=255)
                                                                                         ggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGa                                                           >  1:1773407/1‑149 (MQ=255)
                                                                                           tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAt                                                         <  2:1076025/149‑1 (MQ=255)
                                                                                            gagaTCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                        >  1:1230539/1‑149 (MQ=255)
                                                                                            gagaTCATCCGCGAACGTATCAAACGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                        >  2:4019210/1‑149 (MQ=255)
                                                                                             agaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTg                                                       >  2:527323/1‑149 (MQ=255)
                                                                                                tcatcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTg                                                    >  1:2466675/1‑149 (MQ=255)
                                                                                                tcatcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCAATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTg                                                    >  1:1502777/1‑149 (MQ=255)
                                                                                                 catcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTgg                                                   >  1:4701536/1‑149 (MQ=255)
                                                                                                   tcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGac                                                 >  2:3158685/1‑149 (MQ=255)
                                                                                                   gccGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGac                                                 <  2:4237468/148‑1 (MQ=255)
                                                                                                         aaCGTATCAAACGTTATCAGAGTGACTGCGACTACCGCCCCGGCGGCCTGTGTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaac                                           <  2:3940021/149‑1 (MQ=255)
                                                                                                          aCGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacg                                          >  1:3409554/1‑149 (MQ=255)
                                                                                                          aCGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacg                                          >  1:781461/1‑149 (MQ=255)
                                                                                                           cGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacgc                                         <  2:2819624/149‑1 (MQ=255)
                                                                                                             tATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaacgcca                                       >  1:2134958/1‑149 (MQ=255)
                                                                                                                 aaaCGTTATCAGATTTACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCg                                   >  1:2732152/1‑149 (MQ=255)
                                                                                                                  aaCGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGc                                  >  2:3574297/1‑149 (MQ=255)
                                                                                                                  aaCGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGc                                  >  1:4280391/1‑149 (MQ=255)
                                                                                                                  aaCGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGc                                  >  1:823349/1‑149 (MQ=255)
                                                                                                                   aCGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGcc                                 >  2:1747876/1‑149 (MQ=255)
                                                                                                                        aTCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGaa                            <  2:2165570/149‑1 (MQ=255)
                                                                                                                            gATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAg                        >  2:2114046/1‑149 (MQ=255)
                                                                                                                             aTTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGc                       >  2:462665/1‑149 (MQ=255)
                                                                                                                              ttGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGcc                      >  2:2026256/1‑149 (MQ=255)
                                                                                                                              ttGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGcc                      >  2:2965316/1‑149 (MQ=255)
                                                                                                                                 aCTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGc                   <  2:2319427/149‑1 (MQ=255)
                                                                                                                                 aCTGCGACGACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGc                   <  1:3607088/149‑1 (MQ=255)
                                                                                                                                   tGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGa                 <  1:4204194/149‑1 (MQ=255)
                                                                                                                                       aCTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGc             >  1:266807/1‑149 (MQ=255)
                                                                                                                                          aCCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtat          <  1:2347279/149‑1 (MQ=255)
                                                                                                                                            cGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAc        <  2:3058021/149‑1 (MQ=255)
                                                                                                                                             gCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAcc       <  2:68377/149‑1 (MQ=255)
                                                                                                                                              ccccTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCg      >  2:4435602/1‑149 (MQ=255)
                                                                                                                                              ccccTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCg      >  2:1195756/1‑149 (MQ=255)
                                                                                                                                                  ggcggcCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTgc  >  2:2749926/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGTGC  >  CP009273/1358709‑1359003

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: