Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,514,177 G→T 100% E298D (GAG→GAT)  curA → curcumin/dihydrocurcumin reductase, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,514,1770GT98.9% 859.2 / ‑0.0 270E298D (GAG→GAT) curAcurcumin/dihydrocurcumin reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (145/122);  total (148/122)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 270 total aligned reads displayed.

GTGAGCAGCTATAACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAG  >  CP009273/1514031‑1514323
                                                                                                                                                  |                                                                                                                                                  
 tGAGCAGCTATAACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATgg                                                                                                                                                >  1:1782063/1‑150 (MQ=255)
        cTATAACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATg                                                                                                                                         >  1:3333960/1‑150 (MQ=255)
         tataACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATgg                                                                                                                                        <  1:3093599/150‑1 (MQ=255)
         tataACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATgg                                                                                                                                        >  2:4880693/1‑150 (MQ=255)
          ataACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCACTggg                                                                                                                                       >  1:621704/1‑150 (MQ=255)
          ataACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATggg                                                                                                                                       <  1:2592982/150‑1 (MQ=255)
           taACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGt                                                                                                                                      >  1:2046229/1‑150 (MQ=255)
           taACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGt                                                                                                                                      <  2:5001212/150‑1 (MQ=255)
            aaCGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTg                                                                                                                                     <  2:3909789/150‑1 (MQ=255)
             aCGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGa                                                                                                                                    <  1:3762191/150‑1 (MQ=255)
               gCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGaaa                                                                                                                                  >  2:4164068/1‑150 (MQ=255)
                cTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAag                                                                                                                                 >  1:4754387/1‑150 (MQ=255)
                    tgagCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGa                                                                                                                             >  1:3618172/2‑150 (MQ=255)
                      agCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATa                                                                                                                           <  2:134201/150‑1 (MQ=255)
                       gCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATaa                                                                                                                          <  1:438833/150‑1 (MQ=255)
                         tACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATaaaa                                                                                                                        >  1:5058786/1‑150 (MQ=255)
                           ccaccCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGATGAAAGAGGATAAAATc                                                                                                                      <  2:2426476/150‑1 (MQ=255)
                            caccCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATcc                                                                                                                     >  2:1813610/1‑150 (MQ=255)
                              cccGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCAc                                                                                                                   >  2:3358304/1‑150 (MQ=255)
                               ccGGTCCGGATCTTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGTCAATGGGTGAAAGAGGATAAAATCCACt                                                                                                                  <  1:1014433/150‑1 (MQ=255)
                                 ggTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTAc                                                                                                                >  2:3591943/1‑150 (MQ=255)
                                  gTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTAcc                                                                                                               >  1:829918/1‑150 (MQ=255)
                                  gTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTAcc                                                                                                               >  1:829793/1‑150 (MQ=255)
                                    ccGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACcgc                                                                                                             <  1:4679068/150‑1 (MQ=255)
                                    ccGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACcgc                                                                                                             <  2:4749147/150‑1 (MQ=255)
                                     cGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACcgcg                                                                                                            <  1:1150798/150‑1 (MQ=255)
                                     cGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACcgcg                                                                                                            <  2:3533454/150‑1 (MQ=255)
                                       gATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCgaa                                                                                                          <  1:4279347/150‑1 (MQ=255)
                                        aTCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCgaag                                                                                                         <  2:3147472/150‑1 (MQ=255)
                                         tCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCgaaga                                                                                                        >  1:1878403/1‑150 (MQ=255)
                                           gTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGaaa                                                                                                      >  1:3383157/1‑150 (MQ=255)
                                             ttACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAAtt                                                                                                    >  1:4626954/1‑150 (MQ=255)
                                              tACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTa                                                                                                   >  2:639834/1‑150 (MQ=255)
                                                   ctGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGTATAAAATCCACTACCGCGAAGAAATTACTGAc                                                                                              >  1:3033902/1‑150 (MQ=255)
                                                   ctGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGAc                                                                                              >  1:3343854/1‑150 (MQ=255)
                                                    tGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACg                                                                                             >  2:2608494/1‑150 (MQ=255)
                                                      ttGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGt                                                                                           <  1:4434242/150‑1 (MQ=255)
                                                        gATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGttt                                                                                         <  1:3850349/150‑1 (MQ=255)
                                                         aTGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTa                                                                                        <  1:1643981/150‑1 (MQ=255)
                                                              tACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGaa                                                                                   <  1:4785859/150‑1 (MQ=255)
                                                                     cTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCca                                                                            >  2:3828238/1‑150 (MQ=255)
                                                                      tGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCcac                                                                           >  2:4048331/1‑150 (MQ=255)
                                                                      tGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCcac                                                                           >  1:1357499/1‑150 (MQ=255)
                                                                      tGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCcac                                                                           <  2:1415079/150‑1 (MQ=255)
                                                                        aaaaaaCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAg                                                                         >  1:2993962/1‑150 (MQ=255)
                                                                         aaaaaCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGa                                                                        >  1:3196732/1‑150 (MQ=255)
                                                                         aaaaaCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGa                                                                        >  1:148402/1‑150 (MQ=255)
                                                                          aaaaCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGAc                                                                       >  2:2573250/1‑150 (MQ=255)
                                                                          aaaaCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGAc                                                                       <  2:1628955/150‑1 (MQ=255)
                                                                            aaCGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAATATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGt                                                                     >  2:544191/1‑150 (MQ=255)
                                                                              cGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGttt                                                                   >  1:2836751/1‑150 (MQ=255)
                                                                               gTATTCTCTTGAAAGGTTTTATTATCGCTCTGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTa                                                                  <  2:2938242/150‑1 (MQ=255)
                                                                               gTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTa                                                                  <  1:3646928/150‑1 (MQ=255)
                                                                               gTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTa                                                                  >  2:2737601/1‑150 (MQ=255)
                                                                               gTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTa                                                                  >  1:1678407/1‑150 (MQ=255)
                                                                                 aTTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATc                                                                >  1:3243220/1‑150 (MQ=255)
                                                                                  ttCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCg                                                               >  2:3202939/1‑150 (MQ=255)
                                                                                  ttCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCg                                                               >  2:3081933/1‑150 (MQ=255)
                                                                                   tCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCgg                                                              >  1:2297377/1‑150 (MQ=255)
                                                                                   tCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCgg                                                              >  2:3173050/1‑150 (MQ=255)
                                                                                      cTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCct                                                           >  2:761988/1‑150 (MQ=255)
                                                                                        tGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgc                                                         <  2:4515043/150‑1 (MQ=255)
                                                                                        tGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgc                                                         <  2:1943981/150‑1 (MQ=255)
                                                                                         gCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgct                                                        <  2:2304601/150‑1 (MQ=255)
                                                                                          cAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgctg                                                       <  2:3641994/150‑1 (MQ=255)
                                                                                          cAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgctg                                                       >  2:1906754/1‑150 (MQ=255)
                                                                                          cAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgctg                                                       >  1:1263356/1‑150 (MQ=255)
                                                                                          cAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgctg                                                       >  2:3785541/1‑150 (MQ=255)
                                                                                          cAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCctgctg                                                       <  2:1553617/150‑1 (MQ=255)
                                                                                           aaGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGa                                                      >  1:4103511/1‑150 (MQ=255)
                                                                                             ggTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATCAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAg                                                    >  1:1366047/1‑150 (MQ=255)
                                                                                                       tCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACt                                          <  1:2801718/150‑1 (MQ=255)
                                                                                                       tCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACt                                          >  1:3738878/1‑150 (MQ=255)
                                                                                                       tCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACt                                          >  2:700727/1‑150 (MQ=255)
                                                                                                           tCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCgg                                      >  2:4260330/1‑150 (MQ=255)
                                                                                                           tCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCgg                                      <  2:4690741/150‑1 (MQ=255)
                                                                                                           tCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCgg                                      <  1:3916531/150‑1 (MQ=255)
                                                                                                            cAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGc                                     >  2:4299489/1‑150 (MQ=255)
                                                                                                            cAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGc                                     >  1:4608648/1‑150 (MQ=255)
                                                                                                             aGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCa                                    >  2:2299411/1‑150 (MQ=255)
                                                                                                             aGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCa                                    <  2:173289/150‑1 (MQ=255)
                                                                                                             aGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCa                                    <  2:501434/150‑1 (MQ=255)
                                                                                                                  tATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAgtg                               <  2:610324/150‑1 (MQ=255)
                                                                                                                  tATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAgtg                               <  1:1078423/150‑1 (MQ=255)
                                                                                                                  tATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAgtg                               >  1:4454429/1‑150 (MQ=255)
                                                                                                                  tATGGTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAgtg                               <  1:3972317/150‑1 (MQ=255)
                                                                                                                     ggTCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAgtggtg                            <  2:177246/150‑1 (MQ=255)
                                                                                                                       tCACCGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGAt                          >  2:4095139/1‑150 (MQ=255)
                                                                                                                           cGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCcgc                      <  1:5101320/150‑1 (MQ=255)
                                                                                                                           cGCATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCcgc                      >  1:2269401/1‑150 (MQ=255)
                                                                                                                             cATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGt                    >  1:3828772/1‑150 (MQ=255)
                                                                                                                             cATCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGt                    >  1:2914193/1‑150 (MQ=255)
                                                                                                                              aTCCATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTg                   <  1:2692864/150‑1 (MQ=255)
                                                                                                                                ccATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGc                 <  1:706323/150‑1 (MQ=255)
                                                                                                                                ccATGAGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGc                 >  1:3497056/1‑150 (MQ=255)
                                                                                                                                     aGTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGtg            <  2:208005/150‑1 (MQ=255)
                                                                                                                                      gTTTCAGAGGGATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGtga           <  1:3561225/150‑1 (MQ=255)
                                                                                                                                               ggATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAg  <  2:1039082/150‑1 (MQ=255)
                                                                                                                                               ggATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAg  <  1:580896/150‑1 (MQ=255)
                                                                                                                                               ggATATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAg  >  1:1163899/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTGAGCAGCTATAACGCTACAGAGCTACCACCCGGTCCGGATCGTTTACCTCTGTTGATGGCTACAGTGCTGAAAAAACGTATTCGCTTGCAAGGTTTTATTATCGCTCAGGATTATGGTCACCGCATCCATGAGTTTCAGAGGGAGATGGGGCAATGGGTGAAAGAGGATAAAATCCACTACCGCGAAGAAATTACTGACGGTTTAGAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAG  >  CP009273/1514031‑1514323

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: ‑