Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,501,571 C→A 93.5% 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,5710CA93.5% 880.4 / 45.6 278G127C (GGC→TGC) phpphosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (7/11);  new base A (137/123);  total (144/134)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 279 total aligned reads displayed.

CCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCA  >  CP009273/3501422‑3501717
                                                                                                                                                     |                                                                                                                                                  
cccGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGcc                                                                                                                                                    >  1:1985354/1‑150 (MQ=255)
  cGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGa                                                                                                                                                  <  2:1740949/150‑1 (MQ=255)
     gggCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGAtct                                                                                                                                               >  2:2565823/1‑150 (MQ=255)
        cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAGTACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCcg                                                                                                                                            >  2:804858/1‑150 (MQ=255)
        cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTCATCTTTCCTTCGCGGGTGCAGCTCTCcg                                                                                                                                            >  2:1204222/1‑150 (MQ=255)
        cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCcg                                                                                                                                            >  2:3208647/1‑150 (MQ=255)
        cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCcg                                                                                                                                            <  1:3117734/150‑1 (MQ=255)
        cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCcg                                                                                                                                            <  2:3355254/150‑1 (MQ=255)
        cTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCcg                                                                                                                                            >  1:365410/1‑150 (MQ=255)
         ttGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCcgc                                                                                                                                           >  1:188191/1‑150 (MQ=255)
            tAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCgat                                                                                                                                        >  2:4194469/1‑150 (MQ=255)
             aGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCgatg                                                                                                                                       >  1:4133169/1‑150 (MQ=255)
               cAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCgatgat                                                                                                                                     >  1:671045/1‑150 (MQ=255)
               cAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCgatgat                                                                                                                                     <  1:2497044/150‑1 (MQ=255)
               cAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCgatgat                                                                                                                                     >  1:4052562/1‑150 (MQ=255)
                aaCGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATc                                                                                                                                    >  2:1059811/1‑150 (MQ=255)
                  cGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATccc                                                                                                                                  <  2:2085302/150‑1 (MQ=255)
                     cAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGc                                                                                                                               >  2:3887571/1‑150 (MQ=255)
                     cAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGc                                                                                                                               >  2:1519281/1‑150 (MQ=255)
                      aGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCt                                                                                                                              <  1:2336972/150‑1 (MQ=255)
                        ttGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCttt                                                                                                                            <  1:655616/150‑1 (MQ=255)
                         tGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTc                                                                                                                           <  1:3440437/150‑1 (MQ=255)
                         tGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCACTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTc                                                                                                                           >  2:850190/1‑150 (MQ=255)
                          gCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCa                                                                                                                          >  2:3874902/1‑150 (MQ=255)
                           cTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAg                                                                                                                         >  1:1010012/1‑150 (MQ=255)
                               aGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTcc                                                                                                                     >  1:2758162/1‑150 (MQ=255)
                                 ccccATCGTGCTGACCGACGTCGGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATCAAGACCTTCTCTTCCAGCGGCGTAATCTGTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGt                                                                                                                   <  2:3955081/150‑1 (MQ=255)
                                 ccccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGt                                                                                                                   <  1:2461061/150‑1 (MQ=255)
                                  cccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTg                                                                                                                  >  2:3827833/1‑150 (MQ=255)
                                  cccATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTg                                                                                                                  >  2:4335323/1‑150 (MQ=255)
                                     aTCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCa                                                                                                               <  1:1452723/150‑1 (MQ=255)
                                       cGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATc                                                                                                             <  2:56072/150‑1 (MQ=255)
                                       cGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATc                                                                                                             >  2:3475112/1‑150 (MQ=255)
                                       cGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATc                                                                                                             <  2:4020948/150‑1 (MQ=255)
                                        gTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCg                                                                                                            <  2:3309790/150‑1 (MQ=255)
                                         tGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGa                                                                                                           <  2:2822215/150‑1 (MQ=255)
                                           cTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATa                                                                                                         >  1:3377659/1‑150 (MQ=255)
                                             gAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAcc                                                                                                       <  1:145527/150‑1 (MQ=255)
                                             gAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAcc                                                                                                       <  1:3265849/150‑1 (MQ=255)
                                              aaCGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAccc                                                                                                      >  2:488352/1‑150 (MQ=255)
                                               acgacgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCt                                                                                                     <  2:2025387/150‑1 (MQ=255)
                                                  acgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGtt                                                                                                  <  1:2936357/150‑1 (MQ=255)
                                                  acgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGtt                                                                                                  <  1:1844838/150‑1 (MQ=255)
                                                  acgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGtt                                                                                                  <  1:2529099/150‑1 (MQ=255)
                                                   cgTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATCAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGAGGATCCCGGCTTGCAGCTCCGTGCCAGCGATACCCTGTTc                                                                                                 >  2:2141646/1‑150 (MQ=255)
                                                    gTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCa                                                                                                >  2:2626004/1‑150 (MQ=255)
                                                    gTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCATAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCa                                                                                                <  2:146375/150‑1 (MQ=255)
                                                      aTGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAAt                                                                                              <  2:954265/150‑1 (MQ=255)
                                                       tGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAAtt                                                                                             <  1:1466157/150‑1 (MQ=255)
                                                       tGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATAACCTGTTCAAtt                                                                                             <  2:1621706/150‑1 (MQ=255)
                                                        gCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAAttt                                                                                            <  1:3326378/150‑1 (MQ=255)
                                                        gCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAAttt                                                                                            <  2:3016816/150‑1 (MQ=255)
                                                         cGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTc                                                                                           >  2:1380536/1‑150 (MQ=255)
                                                              agaTCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGa                                                                                      <  2:2728790/150‑1 (MQ=255)
                                                                   gggCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCTATCTCCGCGATGATCCCGGCTTTCATCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCAtc                                                                                 >  2:638484/1‑150 (MQ=255)
                                                                                  ttATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTcc                                                                  >  1:469462/1‑150 (MQ=255)
                                                                                    aTGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTg                                                                >  2:1272990/1‑150 (MQ=255)
                                                                                          cagcgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCt                                                          <  2:3726226/150‑1 (MQ=255)
                                                                                           agcgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTg                                                         >  1:1390654/1‑150 (MQ=255)
                                                                                             cgcagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCg                                                       <  1:3879925/150‑1 (MQ=255)
                                                                                               cagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCggg                                                     >  2:894476/1‑150 (MQ=255)
                                                                                               cagcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCggg                                                     >  2:1582495/1‑150 (MQ=255)
                                                                                                agcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGt                                                    >  2:3935829/1‑150 (MQ=255)
                                                                                                 gcgGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTc                                                   <  1:850190/150‑1 (MQ=255)
                                                                                                    gCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGcc                                                <  2:2662432/150‑1 (MQ=255)
                                                                                                     cAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCca                                               >  2:3578365/1‑150 (MQ=255)
                                                                                                       ataaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCcaca                                             >  2:125516/1‑150 (MQ=255)
                                                                                                       ataaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCcaca                                             >  2:2302197/1‑150 (MQ=255)
                                                                                                        taaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACAt                                            <  1:3653033/150‑1 (MQ=255)
                                                                                                         aattaCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATg                                           <  1:804858/147‑1 (MQ=255)
                                                                                                         aaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCCTCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATg                                           <  1:481368/150‑1 (MQ=255)
                                                                                                         aaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATg                                           <  1:1585478/150‑1 (MQ=255)
                                                                                                         aaataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATg                                           <  1:2978927/150‑1 (MQ=255)
                                                                                                          aataCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGt                                          >  2:2325042/1‑150 (MQ=255)
                                                                                                            taCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTc                                        <  2:3160382/150‑1 (MQ=255)
                                                                                                             aCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTcc                                       <  1:436855/150‑1 (MQ=255)
                                                                                                             aCCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTcc                                       >  2:1878141/1‑150 (MQ=255)
                                                                                                              ccTTCTCTTCCTGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCg                                      <  1:903547/150‑1 (MQ=255)
                                                                                                              ccTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCg                                      >  1:1782706/1‑150 (MQ=255)
                                                                                                                  ctctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaa                                  <  1:3486085/150‑1 (MQ=255)
                                                                                                                   tctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaa                                 >  2:2065517/1‑150 (MQ=255)
                                                                                                                    ctTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGaaaa                                <  2:73956/150‑1 (MQ=255)
                                                                                                                       ccAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAACAAcg                             >  2:4263560/1‑150 (MQ=255)
                                                                                                                       ccAGCGGCGTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAcg                             >  1:3769374/1‑150 (MQ=255)
                                                                                                                            ggcgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTcc                        >  1:2633276/1‑150 (MQ=255)
                                                                                                                            ggcgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAAAGCGTcc                        >  1:626615/1‑150 (MQ=255)
                                                                                                                             gcgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCt                       <  2:2456878/150‑1 (MQ=255)
                                                                                                                              cgTAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTg                      >  2:664015/1‑150 (MQ=255)
                                                                                                                                tAATCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGt                    >  1:1048402/1‑150 (MQ=255)
                                                                                                                                 aaTCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGta                   <  2:2963620/150‑1 (MQ=255)
                                                                                                                                 aaTCTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGta                   >  1:1206277/1‑150 (MQ=255)
                                                                                                                                    cTTTCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGtaata                >  2:538585/1‑150 (MQ=255)
                                                                                                                                       tCCTTCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTGCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAAcc             >  1:3191651/1‑150 (MQ=255)
                                                                                                                                        ccTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCg            <  2:1207968/150‑1 (MQ=255)
                                                                                                                                           tCGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTa         >  1:3048381/1‑150 (MQ=255)
                                                                                                                                            cGCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGa                                   <  1:4214484/123‑1 (MQ=255)
                                                                                                                                             gCTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACa       >  2:4252825/1‑150 (MQ=255)
                                                                                                                                              cTGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAg      <  2:826455/150‑1 (MQ=255)
                                                                                                                                               tGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAgg     <  2:3164303/150‑1 (MQ=255)
                                                                                                                                               tGGTGCAGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAgg     <  2:1846815/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCCGTGGGCTTGTAGCAACGCCAGTTGCTCCAGCCCCATCGTGCTGAACGACGTATGCGTGGAGATCGGGCGTCCGGTCTGGTTATGCGCCAGCGCAGCGGCAATAAATACCTTCTCTTCCAGCGGCGTAATCTTTCCTTCGCTGGTGCCGATCTCCGCGATGATCCCGGCTTTCAGCTCCGTGCCATCGATACCCTGTTCAATTTCATCGACCATCTCCTGCGCCAGTTCCTGCACGCTGCGGGTCGCCACATGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCA  >  CP009273/3501422‑3501717

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: