Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,972,296 G→T 100% E112D (GAG→GAT wzyE → putative ECA polysaccharide chain elongation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,972,2960GT96.3% 921.7 / 19.8 308E112D (GAG→GATwzyEputative ECA polysaccharide chain elongation protein
Reads supporting (aligned to +/- strand):  ref base G (7/5);  new base T (131/165);  total (138/170)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 308 total aligned reads displayed.

TGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTT  >  CP009273/3972149‑3972442
                                                                                                                                                   |                                                                                                                                                  
  gtgtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCa                                                                                                                                                <  1:3686194/150‑1 (MQ=255)
  gtgtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCa                                                                                                                                                <  1:4438402/150‑1 (MQ=255)
      cgcgCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCt                                                                                                                                            >  1:1835429/1‑150 (MQ=255)
       gcgcCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                           <  2:3216140/150‑1 (MQ=255)
       gcgcCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                           >  1:4903216/1‑150 (MQ=255)
       gcgcCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                           >  1:213447/1‑150 (MQ=255)
        cgccgcCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                          >  2:3686773/1‑150 (MQ=255)
        cgccgcCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                          >  1:5050361/1‑150 (MQ=255)
               aGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGt                                                                                                                                   >  2:3183/1‑150 (MQ=255)
                 aaaTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGa                                                                                                                                 >  2:2982146/1‑150 (MQ=255)
                   aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATc                                                                                                                               <  1:1392022/150‑1 (MQ=255)
                   aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATc                                                                                                                               >  1:4082740/1‑150 (MQ=255)
                     cTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCt                                                                                                                             <  1:3051516/150‑1 (MQ=255)
                        gttgCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGAt                                                                                                                          <  1:1329046/150‑1 (MQ=255)
                            cAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGt                                                                                                                      <  1:913159/150‑1 (MQ=255)
                             aGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                     >  2:1165923/1‑150 (MQ=255)
                             aGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                     >  2:2061414/1‑150 (MQ=255)
                             aGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                     <  2:4427792/150‑1 (MQ=255)
                                cGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcg                                                                                                                  <  2:4103032/150‑1 (MQ=255)
                                cGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcg                                                                                                                  >  1:5123320/1‑150 (MQ=255)
                                 gTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgc                                                                                                                 >  1:2989683/1‑150 (MQ=255)
                                  ttGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgcg                                                                                                                <  2:3359289/150‑1 (MQ=255)
                                   tGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCgcgc                                                                                                               >  1:1653521/1‑150 (MQ=255)
                                        tCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTa                                                                                                          <  1:3781769/150‑1 (MQ=255)
                                        tCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTa                                                                                                          <  2:1827072/150‑1 (MQ=255)
                                         cTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTaa                                                                                                         <  1:802464/150‑1 (MQ=255)
                                         cTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTaa                                                                                                         >  2:511215/1‑150 (MQ=255)
                                         cTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTaa                                                                                                         <  1:4725798/150‑1 (MQ=255)
                                          tGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAg                                                                                                        <  2:1520323/150‑1 (MQ=255)
                                          tGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAg                                                                                                        <  1:2338756/150‑1 (MQ=255)
                                            cGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                      >  1:1231379/1‑150 (MQ=255)
                                               gctgctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCg                                                                                                   <  2:1279983/150‑1 (MQ=255)
                                               gctgctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCg                                                                                                   <  1:3361267/150‑1 (MQ=255)
                                               gctgctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCg                                                                                                   <  1:216445/150‑1 (MQ=255)
                                               gctgctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCg                                                                                                   <  1:3837072/150‑1 (MQ=255)
                                                 tgctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGc                                                                                                 >  1:4658125/1‑150 (MQ=255)
                                                 tgctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGc                                                                                                 >  1:4546003/1‑150 (MQ=255)
                                                  gctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCa                                                                                                >  2:2704008/1‑150 (MQ=255)
                                                    ttCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtc                                                                                              >  2:2882633/1‑150 (MQ=255)
                                                         cgcgGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATcttctt                                                                                         >  2:1871598/1‑150 (MQ=255)
                                                         cgcgGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATcttctt                                                                                         <  1:3331295/150‑1 (MQ=255)
                                                         cgcgGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATcttctt                                                                                         <  1:1933259/150‑1 (MQ=255)
                                                            ggTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCAt                                                                                      <  2:2803653/150‑1 (MQ=255)
                                                            ggTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCAt                                                                                      >  2:844276/1‑150 (MQ=255)
                                                             gTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATg                                                                                     >  2:1705886/1‑150 (MQ=255)
                                                               ttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGca                                                                                   >  2:306872/1‑150 (MQ=255)
                                                                  cTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACaa                                                                                <  2:4523150/150‑1 (MQ=255)
                                                                    aTGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACg                                                                              <  2:1381237/150‑1 (MQ=255)
                                                                     tGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACgg                                                                             <  2:3626653/150‑1 (MQ=255)
                                                                       tCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCt                                                                           >  1:2804061/1‑150 (MQ=255)
                                                                        cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtt                                                                          <  1:2475677/150‑1 (MQ=255)
                                                                          ccTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtttt                                                                        >  1:772851/1‑150 (MQ=255)
                                                                             aCAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGc                                                                     >  2:4939097/1‑150 (MQ=255)
                                                                                 aaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGtt                                                                 <  1:2606149/150‑1 (MQ=255)
                                                                                 aaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGtt                                                                 >  2:631020/1‑150 (MQ=255)
                                                                                  aCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTc                                                                <  2:4252207/150‑1 (MQ=255)
                                                                                  aCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTc                                                                >  2:350090/1‑150 (MQ=255)
                                                                                     cGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCgg                                                             <  1:4178863/150‑1 (MQ=255)
                                                                                      gCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGc                                                            >  2:4074653/1‑150 (MQ=255)
                                                                                        cTACACAAACGCGTTGCTGATGTCCCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCtt                                                          <  2:4762930/150‑1 (MQ=255)
                                                                                         tACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTa                                                         >  2:1622482/1‑150 (MQ=255)
                                                                                         tACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTa                                                         <  1:1551964/150‑1 (MQ=255)
                                                                                           cGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAAc                                                       <  2:3819216/150‑1 (MQ=255)
                                                                                           cGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAAc                                                       <  1:544276/150‑1 (MQ=255)
                                                                                              aaaCGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTcc                                                    >  2:2832295/1‑150 (MQ=255)
                                                                                                   cgTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAg                                               <  1:604455/150‑1 (MQ=255)
                                                                                                   cgTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAg                                               <  2:3146173/150‑1 (MQ=255)
                                                                                                     ttGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTc                                             <  1:3185647/150‑1 (MQ=255)
                                                                                                      tGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCa                                            <  1:3606658/150‑1 (MQ=255)
                                                                                                         tGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGAt                                         >  2:1206903/1‑150 (MQ=255)
                                                                                                            tGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCtt                                      >  2:485911/1‑150 (MQ=255)
                                                                                                             gTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCttt                                     <  1:3890293/150‑1 (MQ=255)
                                                                                                             gTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCttt                                     <  2:2733158/150‑1 (MQ=255)
                                                                                                              tACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCtttt                                    >  1:3522794/1‑150 (MQ=255)
                                                                                                              tACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCtttt                                    >  2:2344960/1‑150 (MQ=255)
                                                                                                               aCCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTNGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTc                                   <  1:2676237/150‑1 (MQ=255)
                                                                                                                  gcgcCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAg                                <  2:564827/150‑1 (MQ=255)
                                                                                                                   cgcCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGt                               <  2:3108097/150‑1 (MQ=255)
                                                                                                                      cGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGaa                            >  2:398377/1‑150 (MQ=255)
                                                                                                                       gTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAg                           <  2:3262155/150‑1 (MQ=255)
                                                                                                                        tCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCATACAGTCAGATCTTTTCCAGTGAAGt                          >  2:4952639/1‑150 (MQ=255)
                                                                                                                         ccGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGGGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGtc                         <  2:4485627/150‑1 (MQ=255)
                                                                                                                            atGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTcc                      <  1:398632/149‑1 (MQ=255)
                                                                                                                             tGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCg                     <  1:2544222/150‑1 (MQ=255)
                                                                                                                              gTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCgg                    <  1:2012884/150‑1 (MQ=255)
                                                                                                                              gTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCgg                    >  2:910453/1‑150 (MQ=255)
                                                                                                                               tttACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggc                   <  1:4302914/150‑1 (MQ=255)
                                                                                                                                  aCCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtg                >  1:2847667/1‑150 (MQ=255)
                                                                                                                                    cATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTTAAGTCTCCggcgtggc              <  2:1758416/150‑1 (MQ=255)
                                                                                                                                    cATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtggc              <  2:5140993/150‑1 (MQ=255)
                                                                                                                                       gAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACCGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGtt           <  2:5011480/150‑1 (MQ=255)
                                                                                                                                       gAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGtt           <  1:3424591/150‑1 (MQ=255)
                                                                                                                                       gAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGtt           <  1:2539322/150‑1 (MQ=255)
                                                                                                                                       gAACCGAGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACTACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGtt           <  1:1706763/150‑1 (MQ=255)
                                                                                                                                        aaCCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTa          <  1:309768/150‑1 (MQ=255)
                                                                                                                                        aaCCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTa          <  2:5136182/150‑1 (MQ=255)
                                                                                                                                        aaCCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTa          <  1:1108186/150‑1 (MQ=255)
                                                                                                                                             cgTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACg     <  2:3948571/150‑1 (MQ=255)
                                                                                                                                             cgTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACg     >  2:2437243/1‑150 (MQ=255)
                                                                                                                                               tGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGct   <  1:1185748/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTT  >  CP009273/3972149‑3972442

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: