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,2960GT99.0% 642.8 / ‑1.8 201E112D (GAG→GATwzyEputative ECA polysaccharide chain elongation protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (94/104);  total (95/106)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 201 total aligned reads displayed.

GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGC  >  CP009273/3972151‑3972440
                                                                                                                                                 |                                                                                                                                                
  gtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCa                                                                                                                                              <  1:667647/148‑1 (MQ=255)
   tCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCaa                                                                                                                                             <  2:1783568/148‑1 (MQ=255)
    cgcgCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAAt                                                                                                                                            <  1:2719422/148‑1 (MQ=255)
    cgcgCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAAt                                                                                                                                            <  1:942121/148‑1 (MQ=255)
     gcgcCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATc                                                                                                                                           <  2:301149/148‑1 (MQ=255)
     gcgcCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATc                                                                                                                                           <  2:2525627/148‑1 (MQ=255)
       gccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                         <  2:2495677/148‑1 (MQ=255)
       gccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                         <  1:902631/148‑1 (MQ=255)
       gccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                         <  1:903085/148‑1 (MQ=255)
        ccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                        >  1:576834/1‑148 (MQ=255)
           ccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGt                                                                                                                                     <  2:123411/148‑1 (MQ=255)
             aGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTgg                                                                                                                                   <  2:2143492/148‑1 (MQ=255)
              gAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTggg                                                                                                                                  >  2:3175935/1‑148 (MQ=255)
               aaaTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGt                                                                                                                                 >  1:1825652/1‑148 (MQ=255)
                 aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGa                                                                                                                               >  1:200744/1‑148 (MQ=255)
                 aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGa                                                                                                                               >  2:2182321/1‑148 (MQ=255)
                 aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCAAGAACCGCGTGGATACCAATCTTAAGTGGGTGa                                                                                                                               >  1:1947854/1‑148 (MQ=255)
                     tgttgCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCt                                                                                                                           >  2:2628736/1‑148 (MQ=255)
                      gttgCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTg                                                                                                                          <  1:58091/148‑1 (MQ=255)
                           aGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATggg                                                                                                                     <  2:139209/148‑1 (MQ=255)
                           aGGCGGTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATggg                                                                                                                     <  2:2029174/148‑1 (MQ=255)
                            ggCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGt                                                                                                                    <  2:2406852/148‑1 (MQ=255)
                             gCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                   >  1:989839/1‑148 (MQ=255)
                             gCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                   <  1:1055690/148‑1 (MQ=255)
                              cGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCAACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATTAAACGCGTGGAAACCAATCTTACGTGGGTGATCCTGATGGGTAt                                                                                                                  >  1:753553/1‑148 (MQ=255)
                                ttGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcg                                                                                                                >  1:87485/1‑148 (MQ=255)
                                 tGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgc                                                                                                               >  1:282777/1‑148 (MQ=255)
                                  gCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgcg                                                                                                              >  1:3243916/1‑148 (MQ=255)
                                    tttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCt                                                                                                            <  2:2372272/148‑1 (MQ=255)
                                    tttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCt                                                                                                            >  1:3342422/1‑148 (MQ=255)
                                    tttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCt                                                                                                            <  2:3293089/148‑1 (MQ=255)
                                     ttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTTCAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCTATCTTACGTGGGTGATCCTGATGGGTATCGCGCTg                                                                                                           <  1:1382106/148‑1 (MQ=255)
                                       cTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGAGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGt                                                                                                         <  1:2112546/148‑1 (MQ=255)
                                         gcgggcggCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTaa                                                                                                       <  1:1809367/148‑1 (MQ=255)
                                            ggCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                    <  2:1623400/148‑1 (MQ=255)
                                            ggCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                    <  1:2967273/148‑1 (MQ=255)
                                                    cTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtc                                                                                            >  1:2305520/1‑148 (MQ=255)
                                                     tACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtct                                                                                           <  1:424526/148‑1 (MQ=255)
                                                      aCGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtctt                                                                                          >  1:3419063/1‑148 (MQ=255)
                                                       cgcgGTTTACTATGTCACCGACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtcttc                                                                                         <  2:1908966/148‑1 (MQ=255)
                                                          ggTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCttcttc                                                                                      >  1:1188670/1‑148 (MQ=255)
                                                            tttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCAt                                                                                    >  2:1723607/1‑148 (MQ=255)
                                                            tttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCAt                                                                                    >  2:474730/1‑148 (MQ=255)
                                                            tttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCACCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCAt                                                                                    >  2:3136988/1‑148 (MQ=255)
                                                             ttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATg                                                                                   <  2:2692086/148‑1 (MQ=255)
                                                             ttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATg                                                                                   >  2:1713810/1‑148 (MQ=255)
                                                                 tATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGcaca                                                                               <  1:1156031/148‑1 (MQ=255)
                                                                 tATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGcaca                                                                               <  1:2789267/148‑1 (MQ=255)
                                                                  aTGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACaa                                                                              <  1:2112100/148‑1 (MQ=255)
                                                                  aTGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACaa                                                                              >  2:1974387/1‑148 (MQ=255)
                                                                    gTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGAACCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACg                                                                            >  1:2280394/1‑148 (MQ=255)
                                                                     tCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACgg                                                                           >  1:323665/1‑148 (MQ=255)
                                                                     tCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACgg                                                                           >  1:325235/1‑148 (MQ=255)
                                                                      cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGc                                                                          >  2:92856/1‑148 (MQ=255)
                                                                        ccTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtt                                                                        <  1:169015/148‑1 (MQ=255)
                                                                        ccTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtt                                                                        <  1:1598549/148‑1 (MQ=255)
                                                                         cTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCttt                                                                       <  2:807234/148‑1 (MQ=255)
                                                                         cTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCttt                                                                       <  2:131776/148‑1 (MQ=255)
                                                                         cTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCttt                                                                       <  2:1444722/148‑1 (MQ=255)
                                                                         cTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCttt                                                                       >  2:2909477/1‑148 (MQ=255)
                                                                          tACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtttt                                                                      >  2:3438168/1‑148 (MQ=255)
                                                                           aCAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCttttt                                                                     >  1:2969534/1‑148 (MQ=255)
                                                                            cAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTg                                                                    <  2:1014822/148‑1 (MQ=255)
                                                                               aaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTg                                                                 >  2:441290/1‑148 (MQ=255)
                                                                               aaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTg                                                                 <  1:1198519/148‑1 (MQ=255)
                                                                                aCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGt                                                                >  2:3329251/1‑148 (MQ=255)
                                                                                 cccGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGtt                                                               <  1:2496798/148‑1 (MQ=255)
                                                                                     ccTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCgg                                                           <  2:91894/148‑1 (MQ=255)
                                                                                       tACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCt                                                         <  2:465186/148‑1 (MQ=255)
                                                                                        aCGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCtt                                                        >  1:1281089/1‑148 (MQ=255)
                                                                                         cGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTa                                                       <  2:638768/148‑1 (MQ=255)
                                                                                         cGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTa                                                       >  2:2088585/1‑148 (MQ=255)
                                                                                         cGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTa                                                       >  1:2141692/1‑148 (MQ=255)
                                                                                          gCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTaa                                                      <  1:2678948/148‑1 (MQ=255)
                                                                                              aCGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTcc                                                  <  2:282777/148‑1 (MQ=255)
                                                                                               cgcgTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                 <  1:1284938/148‑1 (MQ=255)
                                                                                                    tGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGt                                            >  1:944410/1‑148 (MQ=255)
                                                                                                       tGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAg                                         >  1:3535513/1‑148 (MQ=255)
                                                                                                        gATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGa                                        >  2:1051690/1‑148 (MQ=255)
                                                                                                         aTGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGAt                                       >  2:1795278/1‑148 (MQ=255)
                                                                                                           gTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCt                                     <  1:3077297/148‑1 (MQ=255)
                                                                                                                gcgcCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTcc                                <  1:2628736/148‑1 (MQ=255)
                                                                                                                 cgcCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCa                               <  1:1799555/148‑1 (MQ=255)
                                                                                                                  gcCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAg                              <  2:2450052/148‑1 (MQ=255)
                                                                                                                  gcCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAg                              <  2:327428/148‑1 (MQ=255)
                                                                                                                    cGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTg                            <  1:1985045/148‑1 (MQ=255)
                                                                                                                     gTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGa                           <  2:3484201/148‑1 (MQ=255)
                                                                                                                       ccGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACCGTCAGATCTTTTCCAGTGAAg                         >  2:74963/1‑148 (MQ=255)
                                                                                                                             tttACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCg                   >  1:967994/1‑148 (MQ=255)
                                                                                                                             tttACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCg                   >  1:819936/1‑148 (MQ=255)
                                                                                                                             tttACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCg                   >  1:3332263/1‑148 (MQ=255)
                                                                                                                              ttACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCgg                  <  2:466538/148‑1 (MQ=255)
                                                                                                                                aCCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcg                <  2:989783/148‑1 (MQ=255)
                                                                                                                                 ccATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggagt               >  2:3061674/1‑148 (MQ=255)
                                                                                                                                  cATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtg              <  1:695801/148‑1 (MQ=255)
                                                                                                                                  cATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtg              >  1:2542472/1‑148 (MQ=255)
                                                                                                                                     gAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtggcg           <  2:1947854/148‑1 (MQ=255)
                                                                                                                                        ccGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTa        >  2:1315379/1‑148 (MQ=255)
                                                                                                                                          gcgTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTaaa      <  2:2517615/148‑1 (MQ=255)
                                                                                                                                          gcgTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTaaa      <  1:1124676/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGC  >  CP009273/3972151‑3972440

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: