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.2% 623.9 / 14.6 208E112D (GAG→GATwzyEputative ECA polysaccharide chain elongation protein
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base T (77/123);  total (80/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 208 total aligned reads displayed.

GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACG  >  CP009273/3972151‑3972439
                                                                                                                                                 |                                                                                                                                               
 tgtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCGACAAAACCCGCCTACGCAAACGCGTTGCTGATGGACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATAcc                                                                                                                                              >  1:2401763/1‑148 (MQ=255)
  gtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCa                                                                                                                                             <  1:2611326/148‑1 (MQ=255)
  gtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCa                                                                                                                                             <  1:2612305/148‑1 (MQ=255)
    cgcgCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAAt                                                                                                                                           <  2:1112051/148‑1 (MQ=255)
    cgcgCCGCCAGAAATCTTGGTGGAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAAt                                                                                                                                           <  2:3020465/148‑1 (MQ=255)
     gcgcCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCACTc                                                                                                                                          <  2:2844067/148‑1 (MQ=255)
      cgccgcCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCt                                                                                                                                         <  1:2188722/148‑1 (MQ=255)
      cgccgcCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCt                                                                                                                                         >  2:422241/1‑148 (MQ=255)
       gccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                        <  1:3489689/148‑1 (MQ=255)
       gccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                        <  2:717953/148‑1 (MQ=255)
        ccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                       <  1:3263367/148‑1 (MQ=255)
        ccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                       <  1:2825540/148‑1 (MQ=255)
        ccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                       <  2:3423642/148‑1 (MQ=255)
        ccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                       <  2:3423598/148‑1 (MQ=255)
            cAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTg                                                                                                                                   <  2:1698722/148‑1 (MQ=255)
               aaaTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGt                                                                                                                                >  2:361766/1‑148 (MQ=255)
               aaaTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGt                                                                                                                                >  2:2466621/1‑148 (MQ=255)
                aaTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTg                                                                                                                               >  1:425337/1‑148 (MQ=255)
                 aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGa                                                                                                                              >  1:3229591/1‑148 (MQ=255)
                 aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGa                                                                                                                              >  1:879583/1‑148 (MQ=255)
                     tgttgCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCt                                                                                                                          >  1:696480/1‑148 (MQ=255)
                         gCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATg                                                                                                                      <  1:2470180/148‑1 (MQ=255)
                         gCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGAGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATg                                                                                                                      <  1:106932/148‑1 (MQ=255)
                          cAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATgg                                                                                                                     <  1:3245013/148‑1 (MQ=255)
                          cAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATgg                                                                                                                     <  1:238351/148‑1 (MQ=255)
                           aGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATggg                                                                                                                    <  1:2122672/148‑1 (MQ=255)
                             gCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                  <  2:1653892/148‑1 (MQ=255)
                             gCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                  >  1:1433346/1‑148 (MQ=255)
                               gTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATc                                                                                                                >  2:2390012/1‑148 (MQ=255)
                                ttGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcg                                                                                                               >  1:2063368/1‑148 (MQ=255)
                                  gCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgcg                                                                                                             >  1:1234256/1‑148 (MQ=255)
                                   cTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCgcgc                                                                                                            >  1:596937/1‑148 (MQ=255)
                                    tttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCt                                                                                                           >  2:1602150/1‑148 (MQ=255)
                                       cTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGt                                                                                                        <  1:1300291/148‑1 (MQ=255)
                                          cGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAg                                                                                                     <  1:958740/148‑1 (MQ=255)
                                          cGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAg                                                                                                     <  2:1775208/148‑1 (MQ=255)
                                             gctgctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGt                                                                                                  <  2:1792051/148‑1 (MQ=255)
                                                   tCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAt                                                                                            <  1:425451/148‑1 (MQ=255)
                                                   tCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAt                                                                                            <  2:562013/148‑1 (MQ=255)
                                                    cTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtc                                                                                           <  1:3742204/148‑1 (MQ=255)
                                                    cTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtc                                                                                           <  2:541577/148‑1 (MQ=255)
                                                          ggTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCttcttc                                                                                     <  1:3687762/148‑1 (MQ=255)
                                                           gTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCa                                                                                    <  1:521920/148‑1 (MQ=255)
                                                           gTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCa                                                                                    <  2:1219253/148‑1 (MQ=255)
                                                             ttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATg                                                                                  >  2:2031462/1‑148 (MQ=255)
                                                             ttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATg                                                                                  <  1:1590000/148‑1 (MQ=255)
                                                                 tATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGcaca                                                                              <  2:629655/148‑1 (MQ=255)
                                                                  aTGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACaa                                                                             <  2:2322864/148‑1 (MQ=255)
                                                                  aTGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACaa                                                                             >  2:367529/1‑148 (MQ=255)
                                                                   tGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAAc                                                                            >  2:3055326/1‑148 (MQ=255)
                                                                   tGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAAc                                                                            >  2:1149667/1‑148 (MQ=255)
                                                                    gTCACCTACAAAACCCGCCTACGAAAACGCGTTGCTGATGTACCGCGCCGTCCGCTCTTTACCATGAACCGCGAGGCTACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTAATGAAGAACg                                                                           >  1:1975462/1‑148 (MQ=255)
                                                                      cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGc                                                                         >  2:2562576/1‑148 (MQ=255)
                                                                          tACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtttt                                                                     >  2:1306439/1‑148 (MQ=255)
                                                                          tACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtttt                                                                     <  2:3522965/148‑1 (MQ=255)
                                                                            cAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTg                                                                   <  1:3797774/148‑1 (MQ=255)
                                                                            cAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTg                                                                   <  1:3385179/148‑1 (MQ=255)
                                                                            cAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGAGTCGGCATCTTCTTCATGCACAACGGCTTTTTg                                                                   <  2:1802555/148‑1 (MQ=255)
                                                                             aaaaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGc                                                                  <  2:1024869/148‑1 (MQ=255)
                                                                              aaaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCt                                                                 >  2:667341/1‑148 (MQ=255)
                                                                               aaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTg                                                                >  1:2027694/1‑148 (MQ=255)
                                                                                 cccGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGtt                                                              >  2:3758542/1‑148 (MQ=255)
                                                                                  ccGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTc                                                             <  2:1071993/148‑1 (MQ=255)
                                                                                   cGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTcc                                                            >  2:3052951/1‑148 (MQ=255)
                                                                                      cTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGc                                                         <  1:2466621/148‑1 (MQ=255)
                                                                                       tACGCAACCGCGGTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCt                                                        <  2:3687499/148‑1 (MQ=255)
                                                                                       tACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCt                                                        <  1:2665202/148‑1 (MQ=255)
                                                                                       tACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCt                                                        >  2:1356985/1‑148 (MQ=255)
                                                                                       tACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCt                                                        <  1:42136/148‑1 (MQ=255)
                                                                                       tACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCt                                                        >  2:196135/1‑148 (MQ=255)
                                                                                            aaaCGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACt                                                   <  1:2129792/148‑1 (MQ=255)
                                                                                            aaaCGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACt                                                   <  1:2531856/148‑1 (MQ=255)
                                                                                              aCGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTcc                                                 >  1:2590412/1‑148 (MQ=255)
                                                                                               cgcgTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                >  1:1180698/1‑148 (MQ=255)
                                                                                               cgcgTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                <  1:3053438/148‑1 (MQ=255)
                                                                                                 cgTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTAc                                              <  2:956012/148‑1 (MQ=255)
                                                                                                    tGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGt                                           >  1:3584243/1‑148 (MQ=255)
                                                                                                    tGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGt                                           <  2:1229066/148‑1 (MQ=255)
                                                                                                       tGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAg                                        <  1:1201068/148‑1 (MQ=255)
                                                                                                       tGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAg                                        <  1:807027/148‑1 (MQ=255)
                                                                                                          tGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATc                                     >  2:508468/1‑148 (MQ=255)
                                                                                                           gTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCt                                    >  2:2669903/1‑148 (MQ=255)
                                                                                                           gTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCt                                    <  2:425338/148‑1 (MQ=255)
                                                                                                              ccGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCtttt                                 <  1:1004580/148‑1 (MQ=255)
                                                                                                              ccGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCtttt                                 >  2:262805/1‑148 (MQ=255)
                                                                                                                     gTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGa                          <  1:2949153/148‑1 (MQ=255)
                                                                                                                     gTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGa                          <  1:2949083/148‑1 (MQ=255)
                                                                                                                       ccGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAg                        >  1:741361/1‑148 (MQ=255)
                                                                                                                        cGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGt                       >  1:1668402/1‑148 (MQ=255)
                                                                                                                        cGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGt                       >  1:1957923/1‑148 (MQ=255)
                                                                                                                          cTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGtct                     <  1:2792458/148‑1 (MQ=255)
                                                                                                                            gTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTcc                   >  1:1017037/1‑148 (MQ=255)
                                                                                                                                 ccATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGACGTCTCCggcgt              >  1:1100044/1‑148 (MQ=255)
                                                                                                                                 ccATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgt              <  1:1276094/148‑1 (MQ=255)
                                                                                                                                     gAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtggcg          <  2:2361150/148‑1 (MQ=255)
                                                                                                                                     gAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtggcg          >  1:1356444/1‑148 (MQ=255)
                                                                                                                                        ccGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTa       >  2:2604654/1‑148 (MQ=255)
                                                                                                                                         cgcgTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTaa      <  2:2079086/148‑1 (MQ=255)
                                                                                                                                             tGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACg  <  1:2342910/148‑1 (MQ=255)
                                                                                                                                             tGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACg  <  1:367490/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACG  >  CP009273/3972151‑3972439

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: