Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,184,128 C→T 6.4% S46S (AGC→AGT yqiI → fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,184,1280CT 6.4% 1175.3 / 44.9 332S46S (AGC→AGTyqiIfimbrial protein
Reads supporting (aligned to +/- strand):  ref base C (145/165);  new base T (4/17);  total (149/183)
Fisher's exact test for biased strand distribution p-value = 1.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 336 total aligned reads displayed.

AGATCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCAAGAC  >  CP009273/3183980‑3184275
                                                                                                                                                    |                                                                                                                                                   
  aTCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGa                                                                                                                                                   <  2:1138057/149‑1 (MQ=255)
    cAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATg                                                                                                                                                 <  2:409993/149‑1 (MQ=255)
    cAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATg                                                                                                                                                 <  1:2098564/149‑1 (MQ=255)
     aataatATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATgg                                                                                                                                                >  2:2474043/1‑149 (MQ=255)
     aataatATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATgg                                                                                                                                                >  2:4955511/1‑149 (MQ=255)
      ataataTGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGt                                                                                                                                               <  1:646273/149‑1 (MQ=255)
      ataataTGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGt                                                                                                                                               >  2:1831146/1‑149 (MQ=255)
      ataataTGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGt                                                                                                                                               <  2:4787096/149‑1 (MQ=255)
         atatGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaaa                                                                                                                                            >  1:4923176/1‑149 (MQ=255)
          tatGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaaaa                                                                                                                                           >  2:4416865/1‑149 (MQ=255)
            tGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTaaaaaa                                                                                                                                         <  2:5042270/149‑1 (MQ=255)
              cgcTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAgg                                                                                                                                       <  2:375547/149‑1 (MQ=255)
                  aCTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGAt                                                                                                                                   <  2:4020867/149‑1 (MQ=255)
                      gTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATaaaa                                                                                                                               >  2:1813335/1‑149 (MQ=255)
                      gTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATaaaa                                                                                                                               >  1:765285/1‑149 (MQ=255)
                      gTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGGTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATaaaa                                                                                                                               >  2:4534445/1‑149 (MQ=255)
                        tAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATc                                                                                                                             <  2:2439688/149‑1 (MQ=255)
                         aaTTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCt                                                                                                                            <  1:356600/149‑1 (MQ=255)
                           ttGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTgg                                                                                                                          >  1:4135425/1‑149 (MQ=255)
                                attaCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAg                                                                                                                     >  1:1018774/1‑149 (MQ=255)
                                attaCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAg                                                                                                                     <  1:2746160/149‑1 (MQ=255)
                                 ttaCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGc                                                                                                                    >  1:1286556/1‑149 (MQ=255)
                                 ttaCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGc                                                                                                                    >  1:2487502/1‑149 (MQ=255)
                                   aCTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAg                                                                                                                  >  2:4851511/1‑149 (MQ=255)
                                    cTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGt                                                                                                                 >  1:326622/1‑149 (MQ=255)
                                    cTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGt                                                                                                                 <  2:866831/149‑1 (MQ=255)
                                        ttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAat                                                                                                             >  2:3590129/1‑149 (MQ=255)
                                        ttttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAat                                                                                                             <  2:1408104/149‑1 (MQ=255)
                                          ttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAatat                                                                                                           <  1:2564968/149‑1 (MQ=255)
                                          ttATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAatat                                                                                                           <  1:2382633/149‑1 (MQ=255)
                                           tATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAtt                                                                                                          >  1:414216/1‑149 (MQ=255)
                                            aTGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAttt                                                                                                         >  1:5014444/1‑149 (MQ=255)
                                             tGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAtttt                                                                                                        <  1:2238374/149‑1 (MQ=255)
                                             tGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAtttt                                                                                                        <  1:5081022/149‑1 (MQ=255)
                                             tGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATAtttt                                                                                                        >  1:4423414/1‑149 (MQ=255)
                                              gtgTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTa                                                                                                       <  2:3483936/147‑1 (MQ=255)
                                               ggTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTAt                                                                                                      <  2:2835587/149‑1 (MQ=255)
                                                gTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATg                                                                                                     >  1:3979500/1‑149 (MQ=255)
                                                gTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATg                                                                                                     >  1:4622082/1‑149 (MQ=255)
                                                      gTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGtt                                                                                               <  1:1910332/149‑1 (MQ=255)
                                                      gTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGtt                                                                                               <  2:2988733/149‑1 (MQ=255)
                                                      gTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGAGGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGtt                                                                                               >  2:4672530/1‑149 (MQ=255)
                                                       ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGttt                                                                                              <  2:2080661/149‑1 (MQ=255)
                                                       ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGttt                                                                                              >  2:462965/1‑149 (MQ=255)
                                                       ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGttt                                                                                              >  1:1183915/1‑149 (MQ=255)
                                                       ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGttt                                                                                              >  1:2616153/1‑149 (MQ=255)
                                                       ttCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGttt                                                                                              <  1:1956219/149‑1 (MQ=255)
                                                        tCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGtttt                                                                                             <  1:4070078/149‑1 (MQ=255)
                                                         cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTc                                                                                            <  1:939619/149‑1 (MQ=255)
                                                         cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAATGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTc                                                                                            >  2:1501753/1‑149 (MQ=255)
                                                         cATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTc                                                                                            >  1:5127600/1‑149 (MQ=255)
                                                          aTTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCt                                                                                           <  2:2232195/149‑1 (MQ=255)
                                                          aTTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCt                                                                                           >  2:5138527/1‑149 (MQ=255)
                                                           ttACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTg                                                                                          >  1:2673254/1‑149 (MQ=255)
                                                            tACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGt                                                                                         <  2:1715960/149‑1 (MQ=255)
                                                              ccTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGtaa                                                                                       >  2:2005807/1‑149 (MQ=255)
                                                                   aTGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATg                                                                                  >  1:1809863/1‑149 (MQ=255)
                                                                      ggCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCAc                                                                               <  1:4551417/149‑1 (MQ=255)
                                                                       gCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACt                                                                              >  2:3263008/1‑149 (MQ=255)
                                                                          aTGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTAcc                                                                           <  2:3857038/149‑1 (MQ=255)
                                                                          aTGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTAcc                                                                           >  1:2473672/1‑149 (MQ=255)
                                                                          aTGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTAcc                                                                           >  1:4651309/1‑149 (MQ=255)
                                                                           tGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCa                                                                          <  1:3888901/149‑1 (MQ=255)
                                                                            ggACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAg                                                                         <  2:4878174/149‑1 (MQ=255)
                                                                            ggACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAg                                                                         <  2:124508/149‑1 (MQ=255)
                                                                               cTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATTCACTACCAGGcc                                                                      >  1:3534597/1‑149 (MQ=255)
                                                                                  cTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCAt                                                                   <  1:4470812/149‑1 (MQ=255)
                                                                                    aTGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCAtct                                                                 >  2:155649/1‑149 (MQ=255)
                                                                                           aCATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGgag                                                          >  1:737247/1‑149 (MQ=255)
                                                                                                aaGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCACAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAgaagaa                                                     >  2:528090/1‑149 (MQ=255)
                                                                                                    tGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATAcc                                                 <  2:599996/149‑1 (MQ=255)
                                                                                                     tgAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATAccc                                                <  1:4180874/148‑1 (MQ=255)
                                                                                                     ggAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATAccc                                                >  1:323885/1‑149 (MQ=255)
                                                                                                        aaTACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATc                                             <  1:2550941/149‑1 (MQ=255)
                                                                                                         aTACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCt                                            <  2:1880692/149‑1 (MQ=255)
                                                                                                          tACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTg                                           >  2:426847/1‑149 (MQ=255)
                                                                                                           aCCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGc                                          <  1:3185356/149‑1 (MQ=255)
                                                                                                           aCCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGAAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGc                                          >  2:3327757/1‑149 (MQ=255)
                                                                                                               ttGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGTGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATgg                                      <  2:759748/149‑1 (MQ=255)
                                                                                                               ttGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATgg                                      >  2:4174535/1‑149 (MQ=255)
                                                                                                               ttGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATgg                                      >  1:4741316/1‑149 (MQ=255)
                                                                                                                  tCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCa                                   >  1:1377689/1‑149 (MQ=255)
                                                                                                                  tCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCa                                   >  1:529964/1‑149 (MQ=255)
                                                                                                                  tCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCa                                   >  1:1595496/1‑149 (MQ=255)
                                                                                                                    aTTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGAAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCaaa                                 >  2:1182706/1‑149 (MQ=255)
                                                                                                                      tGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAata                               <  1:2006614/149‑1 (MQ=255)
                                                                                                                         ttACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATaa                            <  2:4344216/149‑1 (MQ=255)
                                                                                                                          tACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAat                           <  2:4388631/149‑1 (MQ=255)
                                                                                                                           aCGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAatg                          <  2:2595036/149‑1 (MQ=255)
                                                                                                                              tACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAatgatg                       >  1:1869325/1‑149 (MQ=255)
                                                                                                                              tACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAatgatg                       <  2:2321228/149‑1 (MQ=255)
                                                                                                                                  aaGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGAtgttgt                   >  2:3179897/1‑149 (MQ=255)
                                                                                                                                   aGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATgttgtt                  >  2:1977451/1‑149 (MQ=255)
                                                                                                                                     aaTCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAg                <  2:307619/149‑1 (MQ=255)
                                                                                                                                        ctctATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCAt             >  1:2839230/1‑149 (MQ=255)
                                                                                                                                         tctATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATc            >  1:3660552/1‑149 (MQ=255)
                                                                                                                                              ccccAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTc       <  1:3323488/149‑1 (MQ=255)
                                                                                                                                                 cAGCGATGGTAAAAAAGGTGATAACATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACGACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCAAg    <  2:4041993/149‑1 (MQ=255)
                                                                                                                                                  aGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCAAGa   >  2:5110291/1‑149 (MQ=255)
                                                                                                                                                   gCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCAAGAc  >  2:3698878/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGATCAATAATATGCGCTACTTGTTAATTGTTATTACTTTTTTTATGGGTTTTAGTTCATTACCTGCATGGGCAATGGACTGCTATGCTGAACATGAAGGTGGAAATACCGTTGTCATTGGTTACGTACCAAGAATCTCTATCCCCAGCGATGGTAAAAAAGGTGATAAAATCTGGCAAAGCAGTGAATATTTTATGAATGTTTTCTGTAATAATGCACTACCAGGCCCATCTCCAGGAGAAGAATACCCATCTGCATGGGCAAATATAATGATGTTGTTAGCATCAGGTCAAGAC  >  CP009273/3183980‑3184275

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: