Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,972,296 G→T 93.7% 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,2960GT93.7% 579.4 / 28.9 206E112D (GAG→GATwzyEputative ECA polysaccharide chain elongation protein
Reads supporting (aligned to +/- strand):  ref base G (6/7);  new base T (84/109);  total (90/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 207 total aligned reads displayed.

GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCT  >  CP009273/3972151‑3972441
                                                                                                                                                 |                                                                                                                                                 
gtgtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATAc                                                                                                                                                 <  1:2923631/148‑1 (MQ=255)
gtgtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATAc                                                                                                                                                 <  1:1275467/148‑1 (MQ=255)
  gtCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCa                                                                                                                                               <  1:2385579/148‑1 (MQ=255)
   tCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCaa                                                                                                                                              >  2:1981785/1‑148 (MQ=255)
    cgcgCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAAt                                                                                                                                             <  2:1288748/148‑1 (MQ=255)
    cgcgCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAAt                                                                                                                                             <  2:2161830/148‑1 (MQ=255)
       gccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCtt                                                                                                                                          <  2:19063/148‑1 (MQ=255)
        ccgccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTa                                                                                                                                         >  1:113907/1‑148 (MQ=255)
          gccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACg                                                                                                                                       <  1:699272/148‑1 (MQ=255)
          gccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCCCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACg                                                                                                                                       <  1:1953007/148‑1 (MQ=255)
           ccAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGt                                                                                                                                      >  1:3146375/1‑148 (MQ=255)
              gAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTggg                                                                                                                                   <  1:2006407/148‑1 (MQ=255)
               aaaTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGt                                                                                                                                  >  2:379740/1‑148 (MQ=255)
               aaaTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCCCCTAAGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCAAGAACCGCGTGGATACCAATCTTACGTGGGt                                                                                                                                  >  1:1226066/1‑148 (MQ=255)
                 aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGa                                                                                                                                <  2:2803289/148‑1 (MQ=255)
                 aTCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACGATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGa                                                                                                                                <  2:1356393/148‑1 (MQ=255)
                   cTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATc                                                                                                                              <  1:1676210/148‑1 (MQ=255)
                     tgttgCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCt                                                                                                                            >  2:2238139/1‑148 (MQ=255)
                        tgCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAAAGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTAAGTGGGTGATCCTGAt                                                                                                                         >  1:1494891/1‑148 (MQ=255)
                         gCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATg                                                                                                                        <  1:2656094/148‑1 (MQ=255)
                         gCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATg                                                                                                                        <  1:1002433/148‑1 (MQ=255)
                            ggCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGt                                                                                                                     >  1:987555/1‑148 (MQ=255)
                            ggCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGt                                                                                                                     <  1:884518/148‑1 (MQ=255)
                              cGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTAt                                                                                                                   <  1:812280/148‑1 (MQ=255)
                              cGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTAt                                                                                                                   <  1:1609640/148‑1 (MQ=255)
                                ttGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcg                                                                                                                 <  2:1050577/148‑1 (MQ=255)
                                 tGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgc                                                                                                                >  1:67917/1‑148 (MQ=255)
                                 tGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgc                                                                                                                >  2:59954/1‑148 (MQ=255)
                                 tGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATcgc                                                                                                                >  2:892321/1‑148 (MQ=255)
                                    tttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCt                                                                                                             >  1:658476/1‑148 (MQ=255)
                                    tttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCt                                                                                                             <  1:22479/148‑1 (MQ=255)
                                     ttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCTGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTg                                                                                                            <  2:462760/148‑1 (MQ=255)
                                     ttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTg                                                                                                            <  1:109641/148‑1 (MQ=255)
                                     ttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTg                                                                                                            >  2:3143084/1‑148 (MQ=255)
                                     ttCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTg                                                                                                            <  2:47343/148‑1 (MQ=255)
                                       cTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGt                                                                                                          >  1:185729/1‑148 (MQ=255)
                                        tGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTa                                                                                                         <  2:634315/148‑1 (MQ=255)
                                         gCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTaa                                                                                                        <  1:1820842/148‑1 (MQ=255)
                                            ggCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                     >  2:3066872/1‑148 (MQ=255)
                                                gctTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCgg                                                                                                 >  1:3080165/1‑148 (MQ=255)
                                                  ttCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCa                                                                                               <  2:1387160/148‑1 (MQ=255)
                                                   tCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAt                                                                                              >  2:581542/1‑148 (MQ=255)
                                                    cTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtc                                                                                             <  1:781445/148‑1 (MQ=255)
                                                      aCGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtctt                                                                                           >  2:1845108/1‑148 (MQ=255)
                                                      aCGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtctt                                                                                           >  2:175221/1‑148 (MQ=255)
                                                      aCGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGAGGGGTATCGCGCTGGTAAGCGTCGGCAtctt                                                                                           >  2:1644550/1‑148 (MQ=255)
                                                       cgcgGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtcttc                                                                                          <  2:3149099/148‑1 (MQ=255)
                                                         cgGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATcttctt                                                                                        >  1:486711/1‑148 (MQ=255)
                                                          ggTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCttcttc                                                                                       <  1:2804705/148‑1 (MQ=255)
                                                           gTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCa                                                                                      <  2:1995240/148‑1 (MQ=255)
                                                            tttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCAt                                                                                     >  1:2501284/1‑148 (MQ=255)
                                                                 tATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGcaca                                                                                >  1:624155/1‑148 (MQ=255)
                                                                    gTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACg                                                                             >  1:2881879/1‑148 (MQ=255)
                                                                     tCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACgg                                                                            >  2:832615/1‑148 (MQ=255)
                                                                          tACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCtttt                                                                       >  2:1454678/1‑148 (MQ=255)
                                                                           aCAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGTACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCGGGTCAGCGTCGGCATCTTCTTCAGGGACACCGGCttttt                                                                      >  2:1414464/1‑148 (MQ=255)
                                                                             aaaaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGc                                                                    <  2:2975838/148‑1 (MQ=255)
                                                                              aaaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCt                                                                   >  1:660486/1‑148 (MQ=255)
                                                                               aaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTg                                                                  <  2:1524505/148‑1 (MQ=255)
                                                                               aaCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTg                                                                  >  2:1391150/1‑148 (MQ=255)
                                                                                  ccGCCTACGCAAACGCTTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTc                                                               <  2:2110936/148‑1 (MQ=255)
                                                                                   cGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTcc                                                              >  2:3217216/1‑148 (MQ=255)
                                                                                         cGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTa                                                        >  2:2096230/1‑148 (MQ=255)
                                                                                          gCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTaa                                                       <  1:2881742/148‑1 (MQ=255)
                                                                                             aaCGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTc                                                    <  1:260415/148‑1 (MQ=255)
                                                                                              aCGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTcc                                                   >  2:690741/1‑148 (MQ=255)
                                                                                               cgcgTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                  <  1:2384920/148‑1 (MQ=255)
                                                                                                  gTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACa                                               <  2:2578809/148‑1 (MQ=255)
                                                                                                   ttGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAg                                              >  1:3294910/1‑148 (MQ=255)
                                                                                                   ttGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAg                                              <  2:2506887/148‑1 (MQ=255)
                                                                                                    tGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGt                                             >  1:1120932/1‑148 (MQ=255)
                                                                                                    tGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGt                                             <  1:2238139/148‑1 (MQ=255)
                                                                                                        gATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGa                                         <  2:1734652/148‑1 (MQ=255)
                                                                                                         aTGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGAt                                        <  1:2427689/148‑1 (MQ=255)
                                                                                                         aTGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGAt                                        >  1:1352371/1‑148 (MQ=255)
                                                                                                         aTGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGAt                                        >  1:2413435/1‑148 (MQ=255)
                                                                                                           gTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGTGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCt                                      >  1:3232820/1‑148 (MQ=255)
                                                                                                             aCCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCttt                                    <  2:2688543/148‑1 (MQ=255)
                                                                                                             aCCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCttt                                    <  2:2689234/148‑1 (MQ=255)
                                                                                                              ccGCGCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGAGGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCtttt                                   >  1:1924398/1‑148 (MQ=255)
                                                                                                               cgcgCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTc                                  <  1:554570/148‑1 (MQ=255)
                                                                                                               cgcgCCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTc                                  <  1:1497399/148‑1 (MQ=255)
                                                                                                                  gcCGTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAg                               <  1:1981785/148‑1 (MQ=255)
                                                                                                                  gcCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAg                               <  2:2327532/148‑1 (MQ=255)
                                                                                                                     gTCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGa                            <  2:1553404/148‑1 (MQ=255)
                                                                                                                      tCCGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGaa                           >  2:1922220/1‑148 (MQ=255)
                                                                                                                        cGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGt                         >  1:2693825/1‑148 (MQ=255)
                                                                                                                        cGCTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGt                         <  1:2040603/148‑1 (MQ=255)
                                                                                                                          cTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGtct                       <  2:1851260/148‑1 (MQ=255)
                                                                                                                          cTGTTTACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGtct                       <  2:1289060/148‑1 (MQ=255)
                                                                                                                             tttACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCg                    >  1:2132923/1‑148 (MQ=255)
                                                                                                                             tttACCATGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCg                    >  1:547264/1‑148 (MQ=255)
                                                                                                                                   aTGAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtgg              <  2:1073992/148‑1 (MQ=255)
                                                                                                                                     gAACCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtggcg            >  2:182834/1‑148 (MQ=255)
                                                                                                                                     gAACCGAGTGGATACCAATATTACGTGGGTGATCCTGATGGGTATCGCTCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGTCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGTAGTCTCCggcgtgtcg            <  2:1168834/148‑1 (MQ=255)
                                                                                                                                      aaCCGCGTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATCCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtggcgt           <  1:1952809/148‑1 (MQ=255)
                                                                                                                                          gcgTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTaaa       <  1:284354/148‑1 (MQ=255)
                                                                                                                                          gcgTGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTaaa       <  1:1009095/148‑1 (MQ=255)
                                                                                                                                             tGGATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACg    <  2:1120932/148‑1 (MQ=255)
                                                                                                                                              ggATACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGc   <  2:1800731/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCT  >  CP009273/3972151‑3972441

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: