Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,471,705 C→T 100% G51D (GGT→GAT)  slyD ← FKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,471,7050CT97.8% 1028.4 / 8.5 316G51D (GGT→GAT) slyDFKBP‑type peptidyl prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (6/1);  new base T (160/148);  total (167/149)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01
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 316 total aligned reads displayed.

GTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACT  >  CP009273/3471558‑3471851
                                                                                                                                                   |                                                                                                                                                  
gTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCt                                                                                                                                                  >  2:893778/1‑150 (MQ=255)
  ttCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTc                                                                                                                                                <  1:2545399/150‑1 (MQ=255)
  ttCAGACATGAAACGCATAACTAACTTCAGTTCATCAACTCCCATAAATACGTATTTAGGAACACGTTGAACCAGGTTTTCTTCGTACTGACCGTAAGCGTCGTTAGCTCCAACAGCGACATCAAATTTGTAGCCAACTTCATGATCTTc                                                                                                                                                <  2:3451745/150‑1 (MQ=255)
   tCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTcc                                                                                                                                               <  2:1017430/150‑1 (MQ=255)
    cAGCCAGGAATCGCATACCTACCTGCAGTTCATAAACGCCCATAAATACGTCTTTTTGAACACGTTGAACCATTTTTTCGTCGTAATGACCGTAAGCGTAGTTCTCGCCTACAGCGACATCTAATTTGGCGCCAACTTCATGATCTTCCa                                                                                                                                              <  1:2389914/150‑1 (MQ=255)
    cAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCa                                                                                                                                              <  2:4390611/150‑1 (MQ=255)
    cAGCCAGGAAACGCATACCTACCTGAAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTTCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCa                                                                                                                                              <  1:2501835/150‑1 (MQ=255)
     aGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAg                                                                                                                                             <  2:2944255/150‑1 (MQ=255)
       ccAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAgcg                                                                                                                                           <  1:538175/150‑1 (MQ=255)
         aGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGcc                                                                                                                                         <  1:641734/150‑1 (MQ=255)
         aGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGcc                                                                                                                                         >  1:4292694/1‑150 (MQ=255)
          ggAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCg                                                                                                                                        >  2:2561286/1‑150 (MQ=255)
            ataCGCATACCTACCTGCATTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGtt                                                                                                                                      <  1:413907/148‑1 (MQ=255)
            aaaCGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGtt                                                                                                                                      >  1:3653873/1‑150 (MQ=255)
               cGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTcc                                                                                                                                   <  2:216909/150‑1 (MQ=255)
               cGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGATGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACATCATGATCTTCCAGCGCCGGTTcc                                                                                                                                   >  2:827668/1‑150 (MQ=255)
                  aTACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAgg                                                                                                                                >  1:4419822/1‑150 (MQ=255)
                        acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCagag                                                                                                                          >  1:2700748/1‑150 (MQ=255)
                        acctGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCagag                                                                                                                          <  2:4361197/150‑1 (MQ=255)
                             cAGTTCATCAACGCCCATAACTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAg                                                                                                                     <  1:2635548/150‑1 (MQ=255)
                             cAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAg                                                                                                                     <  1:745155/150‑1 (MQ=255)
                              aGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAgg                                                                                                                    >  1:3581665/1‑150 (MQ=255)
                              aGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAgg                                                                                                                    <  1:1401197/150‑1 (MQ=255)
                               gTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAggg                                                                                                                   >  1:5152243/1‑150 (MQ=255)
                                ttCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGa                                                                                                                  <  1:4798806/150‑1 (MQ=255)
                                 tcatcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGGTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGaa                                                                                                                 <  2:3254104/150‑1 (MQ=255)
                                  catcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAAc                                                                                                                >  2:1184282/1‑150 (MQ=255)
                                   atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAAcc                                                                                                               <  2:1912629/150‑1 (MQ=255)
                                   atcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAAcc                                                                                                               <  2:1969300/150‑1 (MQ=255)
                                    tcaACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCg                                                                                                              <  1:1225968/150‑1 (MQ=255)
                                       aCGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGa                                                                                                           <  2:4150763/150‑1 (MQ=255)
                                        cGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCATCGCCGTTTCAAGGCCAGATATCAGGGAACCGTGAc                                                                                                          >  1:1945831/1‑150 (MQ=255)
                                         gCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGAcc                                                                                                         <  1:3957604/150‑1 (MQ=255)
                                         gCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGAcc                                                                                                         <  1:4071650/150‑1 (MQ=255)
                                           ccATAAATCCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAAGTTGTCGCCAACTTCATGATCTTCCATCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCAt                                                                                                       <  2:1168644/150‑1 (MQ=255)
                                           ccATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCAt                                                                                                       <  2:331325/150‑1 (MQ=255)
                                            cATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATg                                                                                                      >  2:3113402/1‑150 (MQ=255)
                                               aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAg                                                                                                   <  2:3672439/150‑1 (MQ=255)
                                               aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAg                                                                                                   >  2:4917000/1‑150 (MQ=255)
                                               aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAg                                                                                                   <  1:2351718/150‑1 (MQ=255)
                                               aaaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAg                                                                                                   >  1:1599222/1‑150 (MQ=255)
                                                aaTACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAgg                                                                                                  >  2:2386152/1‑150 (MQ=255)
                                                   aCGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCCTGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAg                                                                                               <  1:3791006/150‑1 (MQ=255)
                                                    cGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCCGAGATCAGGGAACCGTGACCATGCAGGTAGt                                                                                              >  1:2455609/1‑150 (MQ=255)
                                                     gTCTTGAGGAACACGTTGCACCAGGTGTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTc                                                                                             <  2:1556501/150‑1 (MQ=255)
                                                         ttAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAgc                                                                                         <  1:2025248/150‑1 (MQ=255)
                                                            ggAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCCTGCAGGTAGTCCAgcggc                                                                                      >  2:4443743/1‑150 (MQ=255)
                                                                 acGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTAATGATATTACAGCGCCGTTTCCCGGCCAGAGATCAGGGAACCGTGACCAGGAAGGTAGTCCAGCGGCGCACt                                                                                 >  2:1431932/1‑150 (MQ=255)
                                                                  cGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTGTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTc                                                                                >  2:1651032/1‑150 (MQ=255)
                                                                   gTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCa                                                                               <  1:582165/150‑1 (MQ=255)
                                                                   gTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCa                                                                               >  1:311679/1‑150 (MQ=255)
                                                                     tGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCAcc                                                                             <  1:2271956/150‑1 (MQ=255)
                                                                        aCCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAACTTTGTCGCCAACTTCATGATATTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGTTAGTCCAGCGGCGCACTCACCGga                                                                          <  2:3835324/150‑1 (MQ=255)
                                                                        aCCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGga                                                                          <  1:1036766/150‑1 (MQ=255)
                                                                           aGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGATATCAGGGAACCGTGACCATGCAGGTAGTCCATCGGCGCACTCACCGGAGAc                                                                       >  1:3209436/1‑150 (MQ=255)
                                                                           aGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGAc                                                                       >  1:2882754/1‑150 (MQ=255)
                                                                             gTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAACTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtc                                                                     <  2:5139730/150‑1 (MQ=255)
                                                                             gTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtc                                                                     >  2:3878909/1‑150 (MQ=255)
                                                                             gTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtc                                                                     >  1:2938462/1‑150 (MQ=255)
                                                                               tttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcat                                                                   <  2:2606508/150‑1 (MQ=255)
                                                                               tttCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACtcat                                                                   <  1:1177223/150‑1 (MQ=255)
                                                                                    tcgtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaacc                                                              <  1:3512840/150‑1 (MQ=255)
                                                                                      gtACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATcaaccaa                                                            >  1:5026637/1‑150 (MQ=255)
                                                                                         cTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacac                                                         <  1:3034884/150‑1 (MQ=255)
                                                                                         cTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacac                                                         <  1:3868476/150‑1 (MQ=255)
                                                                                          tGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAacaca                                                        >  2:5245894/1‑150 (MQ=255)
                                                                                            aCCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACAcc                                                      >  2:3151460/1‑150 (MQ=255)
                                                                                             ccGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCg                                                     <  1:1741968/150‑1 (MQ=255)
                                                                                                 aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                                 <  1:280552/150‑1 (MQ=255)
                                                                                                 aaGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtctt                                                 <  1:4378116/150‑1 (MQ=255)
                                                                                                  aGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtcttc                                                >  2:3563966/1‑150 (MQ=255)
                                                                                                   gCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGtcttct                                               <  2:2700748/150‑1 (MQ=255)
                                                                                                    cgtcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTg                                              <  2:2316214/150‑1 (MQ=255)
                                                                                                    cgtcgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATAAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTg                                              >  1:1819493/1‑150 (MQ=255)
                                                                                                       cgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCACCACACCGTCTTCTgtac                                           <  1:884405/150‑1 (MQ=255)
                                                                                                       cgtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtac                                           >  2:1242436/1‑150 (MQ=255)
                                                                                                        gtTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacg                                          >  2:650834/1‑150 (MQ=255)
                                                                                                           cgcgCCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTgtacgtac                                       >  1:328255/1‑150 (MQ=255)
                                                                                                             cgcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCt                                     >  1:540517/1‑150 (MQ=255)
                                                                                                              gcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTAccag                                    >  1:2296171/1‑148 (MQ=255)
                                                                                                              gcCAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTg                                    <  1:5254249/150‑1 (MQ=255)
                                                                                                                cAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                  >  1:4888901/1‑150 (MQ=255)
                                                                                                                cAACAGCGACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGAt                                  >  2:4311388/1‑150 (MQ=255)
                                                                                                                       gACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAg                           <  2:4997804/150‑1 (MQ=255)
                                                                                                                       gACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAg                           <  2:2867947/150‑1 (MQ=255)
                                                                                                                       gACATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAg                           <  2:4635022/150‑1 (MQ=255)
                                                                                                                        aCATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAgg                          >  2:1712334/1‑150 (MQ=255)
                                                                                                                         cATCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGc                         >  1:3001243/1‑150 (MQ=255)
                                                                                                                          aTCAAATTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCt                        >  1:5096858/1‑150 (MQ=255)
                                                                                                                             aaaTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGac                     <  1:4566490/150‑1 (MQ=255)
                                                                                                                              aaTTTGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCCTCTTCAGTACGTACCTGATAGGCCAGGCTGacc                    >  1:5185034/1‑150 (MQ=255)
                                                                                                                                tttGTCGCCAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGaccac                  >  2:3789514/1‑150 (MQ=255)
                                                                                                                                       ccAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCt           <  2:4105567/150‑1 (MQ=255)
                                                                                                                                       ccAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCt           <  1:2151143/150‑1 (MQ=255)
                                                                                                                                       ccAACTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCt           <  2:3996626/150‑1 (MQ=255)
                                                                                                                                           cTTCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTg       <  1:1675885/150‑1 (MQ=255)
                                                                                                                                            ttCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGc      <  1:3987663/150‑1 (MQ=255)
                                                                                                                                             tCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCt     <  2:2505355/150‑1 (MQ=255)
                                                                                                                                             tCATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCt     <  2:1086256/150‑1 (MQ=255)
                                                                                                                                              cATGATCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGTGAACCGTGACCATGCATGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAAAACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTa    >  1:3413007/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTGTCGCCAACTTCATGACCTTCCAGCGCCGTTTCCAGGCCAGAGATCAGGGAACCGTGACCATGCAGGTAGTCCAGCGGCGCACTCACCGGAGACTCATCAACCAACACACCGTCTTCTGTACGTACCTGATAGGCCAGGCTGACCACCAGGTCTTTTGCTACT  >  CP009273/3471558‑3471851

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑