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,8550GT97.6% 1045.9 / 12.1 329G123C (GGC→TGC) puuBgamma‑glutamylputrescine oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base T (187/134);  total (192/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 340 total aligned reads displayed.

GTGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGT  >  CP009273/1358708‑1359001
                                                                                                                                                   |                                                                                                                                                  
    aCTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGcc                                                                                                                                              >  2:1800110/1‑150 (MQ=255)
     cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCt                                                                                                                                             <  2:4272366/150‑1 (MQ=255)
     cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCt                                                                                                                                             <  2:4382856/150‑1 (MQ=255)
     cTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCt                                                                                                                                             <  2:5149308/150‑1 (MQ=255)
      tCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTg                                                                                                                                            <  1:1940516/150‑1 (MQ=255)
       ccTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGt                                                                                                                                           >  1:682935/1‑150 (MQ=255)
        cTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGtt                                                                                                                                          >  2:3260503/1‑150 (MQ=255)
          aCAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTTTGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTtg                                                                                                                                        >  2:377870/1‑150 (MQ=255)
           cAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTtgt                                                                                                                                       <  1:4152116/150‑1 (MQ=255)
             gCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTgg                                                                                                                                     <  1:133082/150‑1 (MQ=255)
              ccGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGc                                                                                                                                    <  1:1205526/150‑1 (MQ=255)
              ccGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGc                                                                                                                                    >  1:4464547/1‑150 (MQ=255)
                 cgACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGat                                                                                                                                 >  2:1342835/1‑150 (MQ=255)
                  gACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGatg                                                                                                                                >  1:139018/1‑150 (MQ=255)
                    cATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGctgac                                                                                                                              >  1:3702667/1‑149 (MQ=255)
                        gACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                          >  1:4740205/1‑150 (MQ=255)
                        gACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGAt                                                                                                                          >  2:1119648/1‑150 (MQ=255)
                          cGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaa                                                                                                                        >  2:2588287/1‑150 (MQ=255)
                           gTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATaaa                                                                                                                       >  2:3079399/1‑150 (MQ=255)
                               tCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGc                                                                                                                   <  1:3311089/150‑1 (MQ=255)
                               tCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGc                                                                                                                   >  1:2077647/1‑150 (MQ=255)
                                 gAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTc                                                                                                                 <  1:1864215/150‑1 (MQ=255)
                                     aaaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCca                                                                                                             >  1:2781126/1‑150 (MQ=255)
                                      aaGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcac                                                                                                            <  2:1956165/150‑1 (MQ=255)
                                       aGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcaca                                                                                                           <  1:4443053/150‑1 (MQ=255)
                                       aGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcaca                                                                                                           >  1:1001695/1‑150 (MQ=255)
                                       aGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCcaca                                                                                                           >  1:1825608/1‑150 (MQ=255)
                                        gCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCacac                                                                                                          >  1:921010/1‑150 (MQ=255)
                                         cTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACt                                                                                                         >  2:4101125/1‑150 (MQ=255)
                                           aCGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTg                                                                                                       >  2:2694947/1‑150 (MQ=255)
                                                aTGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAagag                                                                                                  >  1:3758571/1‑150 (MQ=255)
                                                 tGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGc                                                                                                 <  2:4710251/150‑1 (MQ=255)
                                                 tGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGc                                                                                                 >  1:2233322/1‑150 (MQ=255)
                                                 tGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGc                                                                                                 <  2:2053655/150‑1 (MQ=255)
                                                  ggACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCa                                                                                                >  1:3004724/1‑150 (MQ=255)
                                                   gACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAg                                                                                               <  1:1016319/150‑1 (MQ=255)
                                                   gACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAg                                                                                               <  1:3237249/150‑1 (MQ=255)
                                                    acacCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGa                                                                                              >  2:2232267/1‑150 (MQ=255)
                                                     cacCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGaa                                                                                             <  2:1889699/150‑1 (MQ=255)
                                                       ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAag                                                                                           >  2:2219289/1‑150 (MQ=255)
                                                       ccGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAag                                                                                           >  1:589993/1‑150 (MQ=255)
                                                        cGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                          >  2:2380000/1‑150 (MQ=255)
                                                        cGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAaga                                                                                          >  1:1537199/1‑150 (MQ=255)
                                                          cccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                        >  2:1594342/1‑150 (MQ=255)
                                                          cccGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAgaga                                                                                        >  1:2629798/1‑150 (MQ=255)
                                                             gTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACt                                                                                     <  2:2563313/150‑1 (MQ=255)
                                                                tGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTggg                                                                                  <  1:4614105/150‑1 (MQ=255)
                                                                   tCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAAc                                                                               >  2:2901595/1‑150 (MQ=255)
                                                                     ggCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGc                                                                             <  1:5096886/150‑1 (MQ=255)
                                                                      gcagcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGTTAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCt                                                                            <  1:4306189/150‑1 (MQ=255)
                                                                      gcagcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCt                                                                            >  2:4881507/1‑150 (MQ=255)
                                                                         gcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACg                                                                         >  1:2846091/1‑150 (MQ=255)
                                                                         gcaTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACg                                                                         <  1:2322843/150‑1 (MQ=255)
                                                                          caTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGACCGCTActg                                                                        >  2:4804618/1‑148 (MQ=255)
                                                                          caTGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACgg                                                                        <  1:1456789/150‑1 (MQ=255)
                                                                           atgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGc                                                                       >  1:1686583/1‑150 (MQ=255)
                                                                           atgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGc                                                                       >  2:4144897/1‑150 (MQ=255)
                                                                           atgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGc                                                                       >  1:1116363/1‑150 (MQ=255)
                                                                            tgatgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCa                                                                      <  2:922711/150‑1 (MQ=255)
                                                                              atgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAAt                                                                    <  2:822706/150‑1 (MQ=255)
                                                                               tgTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGACCGCTACGGCAATa                                                                   >  2:165861/1‑150 (MQ=255)
                                                                                 ttCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATaaa                                                                 >  2:3229892/1‑150 (MQ=255)
                                                                                     aaGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGc                                                             >  1:3574297/1‑150 (MQ=255)
                                                                                      aGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCt                                                            >  2:1076993/1‑150 (MQ=255)
                                                                                      aGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCt                                                            >  1:3785228/1‑150 (MQ=255)
                                                                                       ggtggtGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTg                                                           <  1:4180242/150‑1 (MQ=255)
                                                                                       ggtggtGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTg                                                           <  1:4925387/150‑1 (MQ=255)
                                                                                        gtggtgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTgg                                                          <  2:3577293/150‑1 (MQ=255)
                                                                                            tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAAtt                                                      <  2:2370855/150‑1 (MQ=255)
                                                                                            tgAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTAGGAACGCTACGGCAATAAACAGCTGGAAtt                                                      <  2:617733/150‑1 (MQ=255)
                                                                                             gagaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTg                                                     >  1:2887133/1‑150 (MQ=255)
                                                                                              agaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCACTAAACAGCTGGAATTGc                                                    >  1:679728/1‑150 (MQ=255)
                                                                                               gaTCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCt                                                   >  2:3101901/1‑150 (MQ=255)
                                                                                                   atcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGac                                               <  1:3941386/150‑1 (MQ=255)
                                                                                                   atcCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGac                                               >  2:3919648/1‑150 (MQ=255)
                                                                                                      cgcgAACGTATCAAACGTTATCAGATTGACTGCGACNACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgcc                                            >  1:2674819/1‑150 (MQ=255)
                                                                                                        cgAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGacgccaa                                          >  2:3482856/1‑150 (MQ=255)
                                                                                                               aTCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCAtt                                   >  1:3192630/1‑150 (MQ=255)
                                                                                                                tCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTc                                  <  2:568091/150‑1 (MQ=255)
                                                                                                                tCAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTc                                  >  1:913475/1‑150 (MQ=255)
                                                                                                                 cAAACGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCg                                 <  1:2541685/150‑1 (MQ=255)
                                                                                                                    aCGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCg                              >  1:2201106/1‑150 (MQ=255)
                                                                                                                     cGTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGt                             >  2:2390674/1‑150 (MQ=255)
                                                                                                                      gTTATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTg                            >  2:1231650/1‑150 (MQ=255)
                                                                                                                       ttATCAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGa                           >  2:4498768/1‑150 (MQ=255)
                                                                                                                        tATCAGAGTGACTGCGACGACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGaa                          <  2:507975/150‑1 (MQ=255)
                                                                                                                           cAGATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTa                       >  2:4979573/1‑150 (MQ=255)
                                                                                                                             gATTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGc                     >  2:3300924/1‑150 (MQ=255)
                                                                                                                              aTTGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGcc                    <  2:2284328/150‑1 (MQ=255)
                                                                                                                               ttGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCa                   <  2:1704322/150‑1 (MQ=255)
                                                                                                                                tGACTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAg                  >  1:1772045/1‑150 (MQ=255)
                                                                                                                                  aCTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCg                >  1:552561/1‑150 (MQ=255)
                                                                                                                                   cTGCGACTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGa               >  2:384529/1‑150 (MQ=255)
                                                                                                                                        aCTACCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCt          <  2:1399448/150‑1 (MQ=255)
                                                                                                                                          tACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtat        <  1:5115861/150‑1 (MQ=255)
                                                                                                                                           aCCGCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCtata       <  2:2781126/150‑1 (MQ=255)
                                                                                                                                             cGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATAcc     >  2:4644882/1‑150 (MQ=255)
                                                                                                                                              gCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCg    >  1:725187/1‑150 (MQ=255)
                                                                                                                                              gCCCCTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCg    <  1:3072481/150‑1 (MQ=255)
                                                                                                                                                cccTGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCTGt  >  2:4665243/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTGAACTCCTACAGCCGCGACATCGACGTGATCGAAAAAAGCTACGGCATGGACACCGCCCGTATGCTCGGCAGCATGATGTTCGAAGGTGGTGAGATCATCCGCGAACGTATCAAACGTTATCAGATTGACTGCGACTACCGCCCCGGCGGCCTGTTTGTGGCGATGAATGATAAACAGCTCGCCACACTTGAAGAGCAGAAAGAGAACTGGGAACGCTACGGCAATAAACAGCTGGAATTGCTGGACGCCAACGCCATTCGCCGTGAAGTAGCCAGCGATCGCTATACCGGT  >  CP009273/1358708‑1359001

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: