Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,680,642 C→A 9.6% E68D (GAG→GAT fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,680,6420CA 9.6% 803.2 / 39.3 215E68D (GAG→GATfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base C (102/92);  new base A (13/8);  total (115/100)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 223 total aligned reads displayed.

AGTAATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTAT  >  CP009273/1680493‑1680790
                                                                                                                                                     |                                                                                                                                                    
aGTAATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTtc                                                                                                                                                      >  2:3114951/1‑150 (MQ=255)
    aTTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCtt                                                                                                                                                  <  1:560305/150‑1 (MQ=255)
      tCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCa                                                                                                                                                >  2:1012271/1‑150 (MQ=255)
         cTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGac                                                                                                                                             <  2:2763969/150‑1 (MQ=255)
         cTGGCCCGGTTAGCCAGCACTTCGTTAATGTTCATGTTACTATGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCAGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGac                                                                                                                                             >  2:4015754/1‑150 (MQ=255)
          tGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGaca                                                                                                                                            <  2:2670765/150‑1 (MQ=255)
            gCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGacaac                                                                                                                                          <  1:1924401/150‑1 (MQ=255)
             cccGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGacaaca                                                                                                                                         <  1:3126847/150‑1 (MQ=255)
               cGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAg                                                                                                                                       <  1:3175517/150‑1 (MQ=255)
                ggTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGc                                                                                                                                      <  1:3406575/150‑1 (MQ=255)
                  ttAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCt                                                                                                                                    <  1:2438573/150‑1 (MQ=255)
                  ttAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCt                                                                                                                                    <  2:2457854/150‑1 (MQ=255)
                  ttAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCt                                                                                                                                    <  1:2052437/150‑1 (MQ=255)
                   tAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTa                                                                                                                                   <  1:4045628/150‑1 (MQ=255)
                   tAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTa                                                                                                                                   <  1:1946257/150‑1 (MQ=255)
                     gCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTaaa                                                                                                                                 <  1:1901685/150‑1 (MQ=255)
                     gCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTaaa                                                                                                                                 >  1:3155546/1‑150 (MQ=255)
                     gCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGACAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTaaa                                                                                                                                 <  1:4000553/150‑1 (MQ=255)
                      ccAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAAt                                                                                                                                <  1:853190/150‑1 (MQ=255)
                      ccAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAAt                                                                                                                                >  1:2728761/1‑150 (MQ=255)
                       cAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATc                                                                                                                               >  1:1800015/1‑150 (MQ=255)
                       cAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATc                                                                                                                               <  2:2042515/150‑1 (MQ=255)
                         gCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCtt                                                                                                                             >  1:2121612/1‑150 (MQ=255)
                             ttCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCAtt                                                                                                                         <  2:634496/150‑1 (MQ=255)
                               cGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTaa                                                                                                                       >  2:2550497/1‑150 (MQ=255)
                               cGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTaa                                                                                                                       >  2:3555956/1‑150 (MQ=255)
                                gttcatgttcatGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAAc                                                                                                                      >  1:926844/1‑150 (MQ=255)
                                 ttcatgttcatgTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACt                                                                                                                     <  1:1605365/150‑1 (MQ=255)
                                  tcatgttcatgtTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACtt                                                                                                                    >  1:2804720/1‑150 (MQ=255)
                                  tcatgttcatgtTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACtt                                                                                                                    <  2:1764133/150‑1 (MQ=255)
                                   catgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACttt                                                                                                                   <  1:1476199/150‑1 (MQ=255)
                                    atgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                                  <  1:2520287/150‑1 (MQ=255)
                                    atgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                                  <  1:421057/150‑1 (MQ=255)
                                    atgttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACtttt                                                                                                                  <  2:3935257/150‑1 (MQ=255)
                                      gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGc                                                                                                                >  2:3751393/1‑150 (MQ=255)
                                      gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGc                                                                                                                >  2:1653056/1‑150 (MQ=255)
                                      gttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGc                                                                                                                >  2:1041531/1‑150 (MQ=255)
                                       ttcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGcc                                                                                                               <  2:2614594/150‑1 (MQ=255)
                                        tcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCg                                                                                                              <  2:3866891/150‑1 (MQ=255)
                                        tcatgttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCg                                                                                                              >  1:2470283/1‑150 (MQ=255)
                                         catgttACTTTGCGTGCCGGAGACGGTATGCCAGATAGCCAGCGCTAATTCGTCGTCATGCTGTCATGCCAGTACTTCATCCGCATCCTGACGAATGGCGCTAGATTTCTCTTCAGACAACAAGCCTAAATCTTCATTAAATTATGcctc                                                                                                             >  1:1306939/1‑148 (MQ=255)
                                            gttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGc                                                                                                          >  1:770954/1‑150 (MQ=255)
                                             ttACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCAGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCa                                                                                                         >  2:3156886/1‑150 (MQ=255)
                                                      gTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGtt                                                                                                <  2:2971161/150‑1 (MQ=255)
                                                         ccGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAgc                                                                                             <  1:1060648/150‑1 (MQ=255)
                                                              gCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAg                                                                                        >  1:783366/1‑150 (MQ=255)
                                                                 ggTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAgcgc                                                                                     >  2:2813821/1‑150 (MQ=255)
                                                                      gCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGaa                                                                                <  1:1006949/150‑1 (MQ=255)
                                                                                cAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGt                                                                      <  2:690265/150‑1 (MQ=255)
                                                                                 aGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTg                                                                     <  1:3804585/150‑1 (MQ=255)
                                                                                  gCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTgg                                                                    <  1:3582986/150‑1 (MQ=255)
                                                                                    gggAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCCGCGCATGAATCAGTGAGGTGGGc                                                                  >  1:602661/1‑150 (MQ=255)
                                                                                    gggAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGc                                                                  <  2:3495401/150‑1 (MQ=255)
                                                                                       aaTTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAtt                                                               >  2:2471535/1‑150 (MQ=255)
                                                                                       aaTTCGTCGTCATGCTGTCCTGCCAGGACGTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAtt                                                               <  2:3703475/150‑1 (MQ=255)
                                                                                         ttCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCAtttt                                                             >  2:2788053/1‑150 (MQ=255)
                                                                                           cgtcgtCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtct                                                           >  1:1005023/1‑150 (MQ=255)
                                                                                           cgtcgtCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtct                                                           >  1:1005073/1‑150 (MQ=255)
                                                                                            gtcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                          >  2:2168668/1‑150 (MQ=255)
                                                                                            gtcgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTtctc                                                          >  1:3376927/1‑150 (MQ=255)
                                                                                              cgtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCCTTTTCTCCg                                                        >  2:2261946/1‑150 (MQ=255)
                                                                                               gtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                       <  2:1796602/150‑1 (MQ=255)
                                                                                               gtcATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGt                                                       <  2:3316501/150‑1 (MQ=255)
                                                                                                 cATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCg                                                     >  1:3744524/1‑150 (MQ=255)
                                                                                                 cATGCTGTCCTGCCAGTACTTCATCCACCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCg                                                     <  2:433463/150‑1 (MQ=255)
                                                                                                   tGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGaa                                                   >  2:2969248/1‑150 (MQ=255)
                                                                                                             ccAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAAt                                         <  2:1424513/150‑1 (MQ=255)
                                                                                                             ccAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAAt                                         <  1:2581424/150‑1 (MQ=255)
                                                                                                              cAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATg                                        >  2:3408378/1‑150 (MQ=255)
                                                                                                               aGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGc                                       >  2:2883384/1‑150 (MQ=255)
                                                                                                                gTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCt                                      <  2:1139037/150‑1 (MQ=255)
                                                                                                                   cTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCa                                   <  1:2558696/150‑1 (MQ=255)
                                                                                                                     tCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagc                                 >  1:2792257/1‑150 (MQ=255)
                                                                                                                     tCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagc                                 >  1:169047/1‑150 (MQ=255)
                                                                                                                        tCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCagcgag                              <  1:2411469/150‑1 (MQ=255)
                                                                                                                             cgccTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTg                         >  2:591596/1‑150 (MQ=255)
                                                                                                                             cgccTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTg                         >  1:2841920/1‑150 (MQ=255)
                                                                                                                             cgccTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTg                         >  2:591576/1‑150 (MQ=255)
                                                                                                                             cgccTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTg                         >  2:1961981/1‑150 (MQ=255)
                                                                                                                             cgccTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTg                         >  1:1992068/1‑150 (MQ=255)
                                                                                                                                  gACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGGTGAGttt                    >  2:446075/1‑150 (MQ=255)
                                                                                                                                   aCGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtg                   >  1:1350922/1‑150 (MQ=255)
                                                                                                                                    cGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtgt                  >  1:2978239/1‑150 (MQ=255)
                                                                                                                                    cGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtgt                  <  2:1371104/150‑1 (MQ=255)
                                                                                                                                    cGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtgt                  <  2:209946/150‑1 (MQ=255)
                                                                                                                                    cGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTtgt                  <  2:1213588/150‑1 (MQ=255)
                                                                                                                                        tGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTgcgc              >  2:3356827/1‑150 (MQ=255)
                                                                                                                                        tGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTgcgc              >  2:1337042/1‑150 (MQ=255)
                                                                                                                                         ggCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGcc             >  1:3636414/1‑150 (MQ=255)
                                                                                                                                           cgctcgctTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGcccc           <  1:1630659/150‑1 (MQ=255)
                                                                                                                                            gctcgctTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCca          >  2:3469253/1‑150 (MQ=255)
                                                                                                                                            gctcgctTTATCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCca          >  1:1053268/1‑150 (MQ=255)
                                                                                                                                              tcgctTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCcaca        >  1:799454/1‑150 (MQ=255)
                                                                                                                                                gctTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGc      >  2:2516365/1‑150 (MQ=255)
                                                                                                                                                gctTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTGTGTGCGCCCCACAGc      >  2:1139190/1‑150 (MQ=255)
                                                                                                                                                 ctTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGACATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCt     >  2:3859472/1‑150 (MQ=255)
                                                                                                                                                  tttCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCtt    <  1:2256126/150‑1 (MQ=255)
                                                                                                                                                   ttCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTa   <  2:2569501/150‑1 (MQ=255)
                                                                                                                                                    tctcTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTAt  <  1:2331866/150‑1 (MQ=255)
                                                                                                                                                    tctcTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCGTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTAt  >  2:3957989/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTAATTCACTGGCCCGGTTAGCCAGCACTTCGTTCATGTTCATGTTACTTTGCGTGCCGGAGCCGGTCTGCCAGATAGCCAGCGGGAATTCGTCGTCATGCTGTCCTGCCAGTACTTCATCCGCCGCCTGACGAATGGCGCTCGCTTTCTCTTCAGACAACAAGCCTAAATCTTCATTAACTTTTGCCGCTGCACGCTTGGTTAGCGCCAGCGCATGAATCAGTGAGGTGGGCATTTTCTCCGTCGAAATGCGGAAATGCTCCAGCGAGCGTTGAGTTTGTGCGCCCCACAGCTTAT  >  CP009273/1680493‑1680790

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: