Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,501,571 C→A 48.1% G127C (GGC→TGC)  php ← phosphotriesterase homology protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,501,5710CA48.1% 131.4 / inf 322G127C (GGC→TGC) phpphosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (92/75);  new base A (83/72);  total (175/147)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 322 total aligned reads displayed.

CGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCAC  >  CP009273/3501424‑3501718
                                                                                                                                                   |                                                                                                                                                   
cGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGa                                                                                                                                                   <  1:2537392/150‑1 (MQ=255)
cGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGa                                                                                                                                                   >  1:1060376/1‑150 (MQ=255)
cGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGa                                                                                                                                                   <  2:1456463/150‑1 (MQ=255)
 gTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGAt                                                                                                                                                  <  2:1266819/150‑1 (MQ=255)
 gTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGAt                                                                                                                                                  <  2:4347411/150‑1 (MQ=255)
  tGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGAtc                                                                                                                                                 >  2:2831398/1‑150 (MQ=255)
    ggCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGAtctc                                                                                                                                               >  1:2347332/1‑150 (MQ=255)
      cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCcg                                                                                                                                             >  1:1334371/1‑150 (MQ=255)
      cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCcg                                                                                                                                             >  1:2263797/1‑150 (MQ=255)
      cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCcg                                                                                                                                             <  1:1836200/150‑1 (MQ=255)
         gTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCga                                                                                                                                          <  1:1799173/150‑1 (MQ=255)
           aGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCgatg                                                                                                                                        <  2:3731754/150‑1 (MQ=255)
           aGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCgatg                                                                                                                                        <  1:1874024/150‑1 (MQ=255)
             cAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCgatgat                                                                                                                                      <  2:3709412/150‑1 (MQ=255)
              aaCGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATc                                                                                                                                     <  2:3729344/150‑1 (MQ=255)
                   cAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGc                                                                                                                                >  1:1006083/1‑150 (MQ=255)
                   cAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGc                                                                                                                                >  1:910752/1‑150 (MQ=255)
                       tGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTc                                                                                                                            <  1:2534190/150‑1 (MQ=255)
                       tGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCGTCTCTTCCAGCGTCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTc                                                                                                                            <  1:1232232/150‑1 (MQ=255)
                             aGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTcc                                                                                                                      >  1:1598952/1‑150 (MQ=255)
                              gCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCg                                                                                                                     >  2:890519/1‑150 (MQ=255)
                               ccccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGt                                                                                                                    >  1:1905393/1‑150 (MQ=255)
                                cccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTg                                                                                                                   >  1:579378/1‑150 (MQ=255)
                                cccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTg                                                                                                                   >  2:335464/1‑150 (MQ=255)
                                 ccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGc                                                                                                                  >  1:3673284/1‑150 (MQ=255)
                                 ccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGc                                                                                                                  >  2:402462/1‑150 (MQ=255)
                                 ccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGc                                                                                                                  >  2:2779851/1‑150 (MQ=255)
                                         cTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATa                                                                                                          <  2:1528219/150‑1 (MQ=255)
                                          tGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAc                                                                                                         <  1:767363/150‑1 (MQ=255)
                                          tGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAc                                                                                                         <  1:3137103/150‑1 (MQ=255)
                                          tGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAc                                                                                                         >  1:2403424/1‑150 (MQ=255)
                                           gAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAcc                                                                                                        >  2:792786/1‑150 (MQ=255)
                                           gAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAcc                                                                                                        >  2:389873/1‑150 (MQ=255)
                                            aaCGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAccc                                                                                                       >  1:2614592/1‑150 (MQ=255)
                                               gacgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGt                                                                                                    >  1:649574/1‑150 (MQ=255)
                                                 cgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTAGGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTc                                                                                                  >  1:4483977/1‑150 (MQ=255)
                                                  gTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCa                                                                                                 <  1:1636378/150‑1 (MQ=255)
                                                  gTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCa                                                                                                 <  1:3935370/150‑1 (MQ=255)
                                                  gTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCa                                                                                                 <  1:1546072/150‑1 (MQ=255)
                                                   tATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCaa                                                                                                >  2:445024/1‑150 (MQ=255)
                                                   tATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCaa                                                                                                >  1:3451230/1‑150 (MQ=255)
                                                     tGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAAtt                                                                                              >  2:1897852/1‑150 (MQ=255)
                                                     tGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAAtt                                                                                              <  1:1106895/150‑1 (MQ=255)
                                                       cGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTc                                                                                            <  1:3604743/150‑1 (MQ=255)
                                                          ggAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATc                                                                                         >  2:2866602/1‑150 (MQ=255)
                                                          ggAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATc                                                                                         >  2:3145677/1‑150 (MQ=255)
                                                              aTCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGAcc                                                                                     >  2:3190664/1‑150 (MQ=255)
                                                               tCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCa                                                                                    <  1:2829866/150‑1 (MQ=255)
                                                                cGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCTATGATCCCGGCTTTCAGCTCCGTGCCATCTATACCCTGTTCAATTTCATCGACCAt                                                                                   <  1:1274526/150‑1 (MQ=255)
                                                                cGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAt                                                                                   <  2:429365/150‑1 (MQ=255)
                                                                 gggCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAtc                                                                                  >  1:61888/1‑150 (MQ=255)
                                                                 gggCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAtc                                                                                  >  2:2665738/1‑150 (MQ=255)
                                                                 gggCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAtc                                                                                  >  2:4210235/1‑150 (MQ=255)
                                                                     gTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCt                                                                              <  2:4570870/150‑1 (MQ=255)
                                                                        cGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTgcg                                                                           <  2:1163224/150‑1 (MQ=255)
                                                                         ggTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTgcgc                                                                          >  2:2213702/1‑150 (MQ=255)
                                                                         ggTCTGCTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTgcgc                                                                          >  1:4177837/1‑150 (MQ=255)
                                                                            cTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAg                                                                       >  2:1042635/1‑150 (MQ=255)
                                                                            cTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAg                                                                       <  2:2946067/150‑1 (MQ=255)
                                                                            cTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCGTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAg                                                                       <  2:2429570/150‑1 (MQ=255)
                                                                              ggTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGtt                                                                     <  2:1060692/150‑1 (MQ=255)
                                                                               gTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTc                                                                    >  2:207883/1‑150 (MQ=255)
                                                                                ttATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTcc                                                                   >  2:3914707/1‑150 (MQ=255)
                                                                                 tATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCt                                                                  <  2:137607/150‑1 (MQ=255)
                                                                                  aTGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTg                                                                 <  1:335464/150‑1 (MQ=255)
                                                                                       ccAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGc                                                            >  2:2616909/1‑150 (MQ=255)
                                                                                         agcgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTg                                                          >  2:1970593/1‑150 (MQ=255)
                                                                                         agcgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCAGCGACCATCTCCTGCGCCAGTTCGTGCACGCTg                                                          >  1:4377373/1‑150 (MQ=255)
                                                                                           cgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCg                                                        <  2:3952851/150‑1 (MQ=255)
                                                                                           cgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCg                                                        <  2:903717/150‑1 (MQ=255)
                                                                                                cgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCg                                                   <  2:2347332/150‑1 (MQ=255)
                                                                                                cgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCg                                                   <  1:3899626/150‑1 (MQ=255)
                                                                                                cgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCg                                                   <  2:4727717/150‑1 (MQ=255)
                                                                                                 ggCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGc                                                  >  1:2155302/1‑150 (MQ=255)
                                                                                                    aataaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCcac                                               >  1:487178/1‑150 (MQ=255)
                                                                                                       aaataCCTTCTCTTCCCGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATg                                            >  1:2221439/1‑150 (MQ=255)
                                                                                                       aaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATg                                            >  2:1698610/1‑150 (MQ=255)
                                                                                                         ataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGtt                                          <  2:3315007/150‑1 (MQ=255)
                                                                                                         ataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGtt                                          >  1:874675/1‑150 (MQ=255)
                                                                                                              ttCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCggg                                     <  2:3029871/150‑1 (MQ=255)
                                                                                                                ctctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaa                                   >  2:4589474/1‑150 (MQ=255)
                                                                                                                 tctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaa                                  >  1:4130256/1‑150 (MQ=255)
                                                                                                                  ctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaa                                 >  1:3399380/1‑150 (MQ=255)
                                                                                                                    tCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAc                               <  1:2747337/150‑1 (MQ=255)
                                                                                                                     ccAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcg                              <  2:350068/150‑1 (MQ=255)
                                                                                                                       aGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcgcg                            >  1:2987239/1‑150 (MQ=255)
                                                                                                                         cggcgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGACAAACGCGTc                          >  2:2585087/1‑150 (MQ=255)
                                                                                                                         cggcgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTc                          >  1:3020824/1‑150 (MQ=255)
                                                                                                                            cgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCTCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTg                       <  1:1000833/150‑1 (MQ=255)
                                                                                                                            cgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTg                       >  1:1428185/1‑150 (MQ=255)
                                                                                                                              tAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGt                     <  1:4350872/150‑1 (MQ=255)
                                                                                                                              tAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGt                     <  1:1970593/150‑1 (MQ=255)
                                                                                                                              tAATCTTTCCTGCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGt                     <  1:2024036/150‑1 (MQ=255)
                                                                                                                               aaTCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGta                    <  2:649574/150‑1 (MQ=255)
                                                                                                                                  cTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaata                 >  1:3816413/1‑150 (MQ=255)
                                                                                                                                  cTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaata                 <  1:890519/150‑1 (MQ=255)
                                                                                                                                      ccTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCg             <  2:1598952/150‑1 (MQ=255)
                                                                                                                                       cTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGACCAACGCGTCCTGGTAATAACCgg            >  2:2121825/1‑150 (MQ=255)
                                                                                                                                            cTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAg       <  1:2901786/150‑1 (MQ=255)
                                                                                                                                             tGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAgg      >  1:4591357/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCAC  >  CP009273/3501424‑3501718

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: