Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,588,904 A→T 37.0% L287Q (CTG→CAG)  livH ← leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,588,9040AT37.0% 252.5 / 285.0 240L287Q (CTG→CAG) livHleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (66/85);  new base T (40/49);  total (106/134)
Fisher's exact test for biased strand distribution p-value = 8.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 247 total aligned reads displayed.

CAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCATCCCGGC  >  CP009273/3588757‑3589052
                                                                                                                                                   |                                                                                                                                                    
cAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTgg                                                                                                                                                    >  2:2934768/1‑150 (MQ=255)
cAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAgg                                                                                                                                                    <  1:2266556/150‑1 (MQ=255)
  gTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGc                                                                                                                                                  <  1:2189404/150‑1 (MQ=255)
   ttGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCg                                                                                                                                                 <  2:268290/150‑1 (MQ=255)
    tGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGa                                                                                                                                                <  1:2417775/150‑1 (MQ=255)
         gCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATgag                                                                                                                                           >  2:2970704/1‑150 (MQ=255)
         gCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATgag                                                                                                                                           <  2:2997974/150‑1 (MQ=255)
          cccATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATgaga                                                                                                                                          >  1:825032/1‑150 (MQ=255)
            cATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGAcc                                                                                                                                        >  2:3727379/1‑150 (MQ=255)
            cATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGAcc                                                                                                                                        >  2:470110/1‑150 (MQ=255)
             aTAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACca                                                                                                                                       <  1:3973006/150‑1 (MQ=255)
                 aGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATc                                                                                                                                   <  1:1858398/150‑1 (MQ=255)
                     gcccgccAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTta                                                                                                                               <  1:1339200/150‑1 (MQ=255)
                          ccAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTc                                                                                                                          <  2:902047/150‑1 (MQ=255)
                           cAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTcc                                                                                                                         >  2:879723/1‑150 (MQ=255)
                           cAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTcc                                                                                                                         >  2:3126244/1‑150 (MQ=255)
                           aaGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAAAACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTcc                                                                                                                         >  2:1035046/2‑150 (MQ=255)
                            aGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCATCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCg                                                                                                                        >  1:199667/1‑150 (MQ=255)
                             gCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGt                                                                                                                       <  1:868828/150‑1 (MQ=255)
                                caAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACt                                                                                                                    <  1:2078058/150‑1 (MQ=255)
                                     aaCATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGAt                                                                                                               >  1:1477820/1‑150 (MQ=255)
                                      aCATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATa                                                                                                              <  1:3559126/150‑1 (MQ=255)
                                       cATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAg                                                                                                             <  1:2177476/150‑1 (MQ=255)
                                        aTCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAgg                                                                                                            >  1:3500678/1‑150 (MQ=255)
                                          cgcgGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACAGCTTTATATTCCGTACTCAGATAGGCa                                                                                                          >  2:2111098/1‑150 (MQ=255)
                                              gCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCagaag                                                                                                      >  2:1254512/1‑150 (MQ=255)
                                                 gagagCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAagaga                                                                                                   >  1:719782/1‑150 (MQ=255)
                                                   gagCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAgc                                                                                                 <  1:2602335/150‑1 (MQ=255)
                                                    agcagcGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAgcg                                                                                                >  2:4090572/1‑150 (MQ=255)
                                                      cagcGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGcc                                                                                              <  1:987179/150‑1 (MQ=255)
                                                        gcgcCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTc                                                                                            <  1:2330102/150‑1 (MQ=255)
                                                         cgcCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTcc                                                                                           <  1:1711065/150‑1 (MQ=255)
                                                            cATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCTTACTCAGATAGGCAGAAGAGAGCGCCTCCGCa                                                                                        <  2:1472390/150‑1 (MQ=255)
                                                             aTTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCaa                                                                                       >  1:3235673/1‑150 (MQ=255)
                                                               tGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATc                                                                                     >  1:1164398/1‑150 (MQ=255)
                                                                  aaTATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATcccc                                                                                  >  2:1360981/1‑150 (MQ=255)
                                                                  aaTATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATcccc                                                                                  <  2:2442596/150‑1 (MQ=255)
                                                                    tatGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAg                                                                                <  1:694129/150‑1 (MQ=255)
                                                                    tatGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAg                                                                                <  1:4126690/150‑1 (MQ=255)
                                                                    tatGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAg                                                                                <  2:1071342/150‑1 (MQ=255)
                                                                     atGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGACTGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGa                                                                               >  2:400111/1‑150 (MQ=255)
                                                                     atGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGa                                                                               >  2:2695673/1‑150 (MQ=255)
                                                                     atGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGa                                                                               <  1:1000871/150‑1 (MQ=255)
                                                                      tGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGaa                                                                              <  1:2334753/150‑1 (MQ=255)
                                                                        cATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATc                                                                            <  2:178559/150‑1 (MQ=255)
                                                                         aTCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCa                                                                           <  2:2250926/150‑1 (MQ=255)
                                                                         aTCGGTTTCATACTTTCTCGACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCa                                                                           <  2:3142117/150‑1 (MQ=255)
                                                                          tCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAg                                                                          <  1:3557649/150‑1 (MQ=255)
                                                                            ggTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGgc                                                                        <  2:3745211/150‑1 (MQ=255)
                                                                            ggTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGgc                                                                        <  1:1558964/150‑1 (MQ=255)
                                                                             gTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTCTATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGgcc                                                                       >  2:2294950/1‑150 (MQ=255)
                                                                             gTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGgcc                                                                       <  1:3847291/150‑1 (MQ=255)
                                                                             gTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGgcc                                                                       >  2:2622048/1‑150 (MQ=255)
                                                                                 cATTCTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCTGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATTGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCa                                                                   <  2:2950743/150‑1 (MQ=255)
                                                                                 cATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCa                                                                   >  2:1546893/1‑150 (MQ=255)
                                                                                  aTACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCaa                                                                  >  1:4075263/1‑150 (MQ=255)
                                                                                      tttCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAatcat                                                              <  1:1828168/150‑1 (MQ=255)
                                                                                         ctctACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGc                                                           >  1:819183/1‑150 (MQ=255)
                                                                                         ctctACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAAAGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGc                                                           >  1:3340190/1‑150 (MQ=255)
                                                                                               cTCCGGGCGACCCAGCATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCgg                                                     <  1:4074454/150‑1 (MQ=255)
                                                                                                  cGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAAt                                                  <  1:1957156/150‑1 (MQ=255)
                                                                                                      cGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAAtgctg                                              >  1:2378798/1‑150 (MQ=255)
                                                                                                       gACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCACTCCCGGACtgctgc                                             >  2:4101553/1‑150 (MQ=255)
                                                                                                       gACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCACTCCCGGAAtgctgc                                             >  2:2835745/1‑150 (MQ=255)
                                                                                                        aCCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGcc                                            <  2:3176164/150‑1 (MQ=255)
                                                                                                        aCCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGcc                                            >  2:3734841/1‑150 (MQ=255)
                                                                                                         cccAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCa                                           <  1:3332422/150‑1 (MQ=255)
                                                                                                             gAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATcc                                       <  2:1688914/150‑1 (MQ=255)
                                                                                                             gAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATcc                                       <  2:2161539/150‑1 (MQ=255)
                                                                                                             gAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATcc                                       <  2:2161577/150‑1 (MQ=255)
                                                                                                                 aCCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCAcc                                   >  1:1267304/1‑150 (MQ=255)
                                                                                                                 aCCGGTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCAcc                                   >  1:890066/1‑150 (MQ=255)
                                                                                                                    ggTCGGCATCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAg                                >  2:3830734/1‑150 (MQ=255)
                                                                                                                      tCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCa                              <  2:12112/150‑1 (MQ=255)
                                                                                                                       cGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCTCACCGAGCAc                             <  2:3402853/150‑1 (MQ=255)
                                                                                                                        ggCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCAcc                            >  1:518169/1‑150 (MQ=255)
                                                                                                                          cATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCAccgc                          >  2:4149058/1‑150 (MQ=255)
                                                                                                                           aTCACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCAccgcc                         >  1:281524/1‑150 (MQ=255)
                                                                                                                           aTCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCAccgcc                         >  2:1394509/1‑150 (MQ=255)
                                                                                                                             cACCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACcgccgc                       >  2:3669327/1‑150 (MQ=255)
                                                                                                                              aCCAGCAGCACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCcgcg                      >  1:3262545/1‑150 (MQ=255)
                                                                                                                                 agcagcACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTa                   <  2:1961257/150‑1 (MQ=255)
                                                                                                                                 agcagcACCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTa                   <  2:1961078/150‑1 (MQ=255)
                                                                                                                                    agcaCCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAg                >  1:3253552/1‑150 (MQ=255)
                                                                                                                                     gcaCCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAgg               <  1:319758/150‑1 (MQ=255)
                                                                                                                                      caCCAGAATCAGCTGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGc              <  2:3523048/150‑1 (MQ=255)
                                                                                                                                       aCCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCt             >  1:3868102/1‑150 (MQ=255)
                                                                                                                                       aCCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCt             >  1:396708/1‑150 (MQ=255)
                                                                                                                                        ccAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCtt            <  1:1017843/150‑1 (MQ=255)
                                                                                                                                         cAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCttt           <  2:2775393/150‑1 (MQ=255)
                                                                                                                                          aGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTc          <  2:3911539/150‑1 (MQ=255)
                                                                                                                                          aGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTc          >  1:3099607/1‑150 (MQ=255)
                                                                                                                                          aGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTc          >  1:2950073/1‑150 (MQ=255)
                                                                                                                                          aGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTc          <  1:1041378/150‑1 (MQ=255)
                                                                                                                                            aaTCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCAt        <  1:2916333/150‑1 (MQ=255)
                                                                                                                                              tCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCATcc      <  1:1276712/150‑1 (MQ=255)
                                                                                                                                               cagcagGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCATccc     >  1:882597/1‑150 (MQ=255)
                                                                                                                                                agctgGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCATCCCg    >  1:938173/1‑150 (MQ=255)
                                                                                                                                                agcagGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCATCCCg    >  1:1734434/1‑150 (MQ=255)
                                                                                                                                                  cagGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCATCCCGGc  >  2:3469365/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAGTTGCACGCCCATAAAGACGCCCGCCAGCACAAAGAACATCGCGGCAGAGAGCAGCGCCATTGCAATATGCATCGGTTTCATACTTTCTCTACCTCCGGGCGACCCAGAATACCGGTCGGCATCACCAGCAGCACCAGAATCAGCAGGGCGAATGAGACCACATCTTTATATTCCGTACTCAGATAGGCAGAAGAGAGCGCCTCCGCAATCCCCAGAATCAGGCCGCCAATCATCGCTCCCGGAATGCTGCCAATCCCACCGAGCACCGCCGCGGTAAAGGCTTTCATCCCGGC  >  CP009273/3588757‑3589052

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: