Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,124,159 C→A 25.4% A67A (GCG→GCT wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,124,1590CA25.4% 520.3 / 155.8 230A67A (GCG→GCTwcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base C (89/82);  new base A (23/36);  total (112/118)
Fisher's exact test for biased strand distribution p-value = 9.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

AAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAG  >  CP009273/2124010‑2124308
                                                                                                                                                     |                                                                                                                                                     
aaaCGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAc                                                                                                                                                       >  1:1053751/1‑150 (MQ=255)
aaaCGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAc                                                                                                                                                       <  1:3261306/150‑1 (MQ=255)
aaaCGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAc                                                                                                                                                       >  2:361963/1‑150 (MQ=255)
   cGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGcc                                                                                                                                                    >  2:164122/1‑150 (MQ=255)
    gTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCg                                                                                                                                                   <  2:2914816/150‑1 (MQ=255)
      aaTATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTg                                                                                                                                                 >  2:3077474/1‑150 (MQ=255)
      aaTATCTTGAGATGGTTGTGTGGCGATATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTg                                                                                                                                                 <  1:4623110/150‑1 (MQ=255)
        tatCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAg                                                                                                                                               <  1:2839707/150‑1 (MQ=255)
        tatCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAg                                                                                                                                               <  1:1336844/150‑1 (MQ=255)
            ttGCGTTGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTccc                                                                                                                                           <  2:4945654/150‑1 (MQ=255)
             tGCGTTGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAACCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCcg                                                                                                                                          <  2:2402245/150‑1 (MQ=255)
             tGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCcg                                                                                                                                          >  1:870046/1‑150 (MQ=255)
             tGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCcg                                                                                                                                          <  2:2461494/150‑1 (MQ=255)
                gATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCgcgc                                                                                                                                       >  1:2107041/1‑150 (MQ=255)
                   ggTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCaa                                                                                                                                    >  2:5139018/1‑150 (MQ=255)
                     ttGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCCTAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATc                                                                                                                                  >  2:4417804/1‑150 (MQ=255)
                        tgtgGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATATTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCt                                                                                                                               >  1:4761285/1‑150 (MQ=255)
                        tgtgGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCt                                                                                                                               <  2:598902/150‑1 (MQ=255)
                          tgGCTAAATTTTTTTAAAACTAATACTCCACACCATGCAAATAATGATCACAATATAGACAAACAGCATTGTTGAAGTGACTTTTTTACTGGCATCCAGCGACTGACCAAACAGATAGTTACACGCCGTGAGTCCCGCGCCAATCCCTaa                                                                                                                             <  1:457405/150‑1 (MQ=255)
                           ggCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAg                                                                                                                            >  1:1670453/1‑150 (MQ=255)
                           ggCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAg                                                                                                                            >  1:3700078/1‑150 (MQ=255)
                            gCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAgg                                                                                                                           <  1:4329073/150‑1 (MQ=255)
                            gCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAgg                                                                                                                           <  1:3331301/150‑1 (MQ=255)
                                  ttGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAatca                                                                                                                     <  2:2558100/150‑1 (MQ=255)
                                    gTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAatcatc                                                                                                                   >  2:112660/1‑150 (MQ=255)
                                            aCGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAAtttttt                                                                                                           <  1:2132449/150‑1 (MQ=255)
                                             cGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTa                                                                                                          <  2:1545950/150‑1 (MQ=255)
                                                  aCTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCt                                                                                                     >  1:4716994/1‑150 (MQ=255)
                                                     ccACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGAt                                                                                                  <  1:2585705/150‑1 (MQ=255)
                                                     ccACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGAt                                                                                                  <  2:5052576/150‑1 (MQ=255)
                                                     ccACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGAt                                                                                                  <  1:5009181/150‑1 (MQ=255)
                                                       acacCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGc                                                                                                >  2:1781713/1‑150 (MQ=255)
                                                         acCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGt                                                                                              >  2:515159/1‑150 (MQ=255)
                                                         acCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGt                                                                                              >  2:3093169/1‑150 (MQ=255)
                                                           cATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTc                                                                                            >  2:1593282/1‑150 (MQ=255)
                                                               ccAATAATGATCACAATATAGACAAACAGCATAGTTGACGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTaaa                                                                                        <  2:853734/150‑1 (MQ=255)
                                                                 aataatGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTaaaaa                                                                                      >  2:5028570/1‑150 (MQ=255)
                                                                 aataatGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTaaaaa                                                                                      >  2:1508879/1‑150 (MQ=255)
                                                                  ataatGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAAc                                                                                     >  1:2223321/1‑150 (MQ=255)
                                                                   taatGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAAca                                                                                    <  1:946760/150‑1 (MQ=255)
                                                                    aatGATCACAAGATAGACACACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAAcag                                                                                   <  2:4337141/150‑1 (MQ=255)
                                                                     atGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAAcagc                                                                                  >  1:3812660/1‑150 (MQ=255)
                                                                     atGATCACAATAGAGACAACCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAAGCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAAcagc                                                                                  <  2:2225951/150‑1 (MQ=255)
                                                                       gATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGTCATCCAGCGACAGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCAAAAGGCAATCATCAATTTTTTAACGCAGATGCGTTCTAAAAAcagcag                                                                                >  1:3970949/1‑150 (MQ=255)
                                                                       gATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAAcagcag                                                                                <  2:1541660/150‑1 (MQ=255)
                                                                        aTCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGt                                                                               <  2:696794/150‑1 (MQ=255)
                                                                              atatAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAg                                                                         >  2:1950848/1‑150 (MQ=255)
                                                                               tataGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGa                                                                        >  1:3747488/1‑150 (MQ=255)
                                                                                 taGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACa                                                                      <  2:3639809/150‑1 (MQ=255)
                                                                                  aGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAg                                                                     <  2:1053486/150‑1 (MQ=255)
                                                                                      aaaCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTaa                                                                 <  2:2714739/150‑1 (MQ=255)
                                                                                      aaaCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTaa                                                                 >  2:4192151/1‑150 (MQ=255)
                                                                                      aaaCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTaa                                                                 <  1:450083/150‑1 (MQ=255)
                                                                                      aaaCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTaa                                                                 >  1:2060705/1‑150 (MQ=255)
                                                                                      aaaCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTaa                                                                 >  2:4521954/1‑150 (MQ=255)
                                                                                       aaCAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCCGTAACAAGACAGGTaaa                                                                >  2:1301055/1‑150 (MQ=255)
                                                                                          aGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTaaaaaa                                                             >  1:1053343/1‑150 (MQ=255)
                                                                                            cATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGt                                                           >  1:1443967/1‑150 (MQ=255)
                                                                                               aGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGAc                                                        >  1:1275578/1‑150 (MQ=255)
                                                                                               aGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGCTAAAAAAGAGAc                                                        >  2:3751830/1‑150 (MQ=255)
                                                                                                   gAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATg                                                    <  1:3057603/150‑1 (MQ=255)
                                                                                                    aaGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATgg                                                   >  2:4001992/1‑150 (MQ=255)
                                                                                                    aaGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATgg                                                   >  1:991600/1‑150 (MQ=255)
                                                                                                       tGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTCAAAAAGTGACGATGGTGa                                                >  1:812036/1‑150 (MQ=255)
                                                                                                         aCGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATg                                              <  2:3685833/150‑1 (MQ=255)
                                                                                                          cGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATgg                                             <  2:1678633/150‑1 (MQ=255)
                                                                                                            tATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATgggg                                           >  2:4759414/1‑150 (MQ=255)
                                                                                                              tttACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGaa                                         <  1:4601303/150‑1 (MQ=255)
                                                                                                               ttACTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGaaa                                        >  1:1047888/1‑150 (MQ=255)
                                                                                                                  cTGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACtt                                     <  2:3184134/150‑1 (MQ=255)
                                                                                                                   tGGCATCCAGCGACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGAc                                                        <  2:3372348/130‑1 (MQ=255)
                                                                                                                           aGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGACACTTTCGCCAAGc                            >  2:3210260/1‑150 (MQ=255)
                                                                                                                           aGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTGTCGCCAAGc                            >  2:1780713/1‑150 (MQ=255)
                                                                                                                            gCGACTGACCAAACAGATAGTTACAAGCCGTGAGTCCCGCGCCACTCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCt                           <  2:4605045/150‑1 (MQ=255)
                                                                                                                             cGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTGACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACCGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTg                          <  1:1512482/150‑1 (MQ=255)
                                                                                                                             cGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTg                          >  2:2103908/1‑150 (MQ=255)
                                                                                                                              gACTGACCAAACAGATAGTTAACCGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTgg                         <  1:1114859/150‑1 (MQ=255)
                                                                                                                              gACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTgg                         >  2:2879900/1‑150 (MQ=255)
                                                                                                                              gACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAAACCCAAAGGCAATCCTCAATAATTAAACGCTGAAGCGTGCTAAAAACAGCAGAAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGACACTAAAGCCAAGcacg                         >  1:4068728/1‑147 (MQ=255)
                                                                                                                              gACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTgg                         >  2:4700121/1‑150 (MQ=255)
                                                                                                                              gACTGACCAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTgg                         >  1:3050310/1‑150 (MQ=255)
                                                                                                                                  gACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGa                     >  2:3903148/1‑150 (MQ=255)
                                                                                                                                  gACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGa                     >  1:5150667/1‑150 (MQ=255)
                                                                                                                                    ccAAACAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAtt                   >  1:2941225/1‑150 (MQ=255)
                                                                                                                                     cAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAttt                  <  1:3519494/150‑1 (MQ=255)
                                                                                                                                      aaaCAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGAtttt                 <  1:3102418/150‑1 (MQ=255)
                                                                                                                                       aaCAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGGAAAAAAGTGACGATGGTGATGGGGACACTTTCGCCAAGCTGGGCGATTTTg                >  1:4945462/1‑150 (MQ=255)
                                                                                                                                       aaCAGATAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTg                <  2:774286/150‑1 (MQ=255)
                                                                                                                                         cAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCTTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGAc              <  1:17335/150‑1 (MQ=255)
                                                                                                                                         cAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGAc              <  1:2635695/150‑1 (MQ=255)
                                                                                                                                           gATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGt            <  2:3511377/150‑1 (MQ=255)
                                                                                                                                            aTAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGtt           <  2:1828762/150‑1 (MQ=255)
                                                                                                                                            aTAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGtt           <  1:1852529/150‑1 (MQ=255)
                                                                                                                                            aTAGTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGtt           <  1:4449278/150‑1 (MQ=255)
                                                                                                                                              aGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTaa         <  2:662311/150‑1 (MQ=255)
                                                                                                                                               gTTAAAAGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAAc        >  2:537668/1‑150 (MQ=255)
                                                                                                                                                  aaaCGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACcag     >  1:1965021/1‑150 (MQ=255)
                                                                                                                                                  aaaCGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACcag     >  2:2614816/1‑150 (MQ=255)
                                                                                                                                                    aCGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACcagca   <  1:4255939/150‑1 (MQ=255)
                                                                                                                                                     cGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACcagcag  >  2:3966669/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACGTAATATCTTGCGATGGTTGTGTGGCGAAATTGTTTTAAAACGAATACTCCACACCATGCCAATAATGATCACAATATAGACAAACAGCATAGTTGAAGTGACGTATTTACTGGCATCCAGCGACTGACCAAACAGATAGTTAAACGCCGTGAGTCCCGCGCCAATCCCTAAGGCAATCATCAATTTTTTAACGCTGATGCGTTCTAAAAACAGCAGTAACAAGACAGGTAAAAAAGTGACGATGGTGATGGGGAAACTTTCGCCAAGCTGGGCGATTTTGACGTTAACCAGCAG  >  CP009273/2124010‑2124308

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: