Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,969,710 G→T 91.2% G25C (GGT→TGT)  wzxE → O‑antigen translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,969,7100GT91.2% 1161.4 / 96.4 385G25C (GGT→TGT) wzxEO‑antigen translocase
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base T (202/149);  total (219/166)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 385 total aligned reads displayed.

TGCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCG  >  CP009273/3969563‑3969856
                                                                                                                                                   |                                                                                                                                                  
  ccGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTaa                                                                                                                                                <  2:2129189/150‑1 (MQ=255)
     cTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGtt                                                                                                                                             <  2:1909032/150‑1 (MQ=255)
          cTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTgg                                                                                                                                        >  1:3282850/1‑150 (MQ=255)
             cAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCgg                                                                                                                                     >  1:186673/1‑150 (MQ=255)
              aaCCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGgt                                                                                                                                    <  2:1607207/150‑1 (MQ=255)
               aCCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGgtg                                                                                                                                   >  1:437568/1‑150 (MQ=255)
                 cTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTc                                                                                                                                 >  2:4212751/1‑150 (MQ=255)
                     cGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAttt                                                                                                                             <  2:1679461/150‑1 (MQ=255)
                      gCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTg                                                                                                                            >  1:481091/1‑150 (MQ=255)
                          gTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTcc                                                                                                                        >  1:1147926/1‑150 (MQ=255)
                          gTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTcc                                                                                                                        >  1:4152951/1‑150 (MQ=255)
                           tCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCcg                                                                                                                       >  1:2964401/1‑150 (MQ=255)
                                cAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCggg                                                                                                                  >  1:581956/1‑150 (MQ=255)
                                cAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCggg                                                                                                                  >  1:4441582/1‑150 (MQ=255)
                                 aGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCgggg                                                                                                                 >  1:3595838/1‑150 (MQ=255)
                                 aGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCgggg                                                                                                                 >  1:4625532/1‑150 (MQ=255)
                                 aGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCgggg                                                                                                                 >  2:4147807/1‑150 (MQ=255)
                                    gTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCtt                                                                                                              >  1:1701101/1‑150 (MQ=255)
                                     tACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTg                                                                                                             <  2:4750036/150‑1 (MQ=255)
                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                           <  2:3248435/150‑1 (MQ=255)
                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                           >  1:3058531/1‑150 (MQ=255)
                                        ggTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGc                                                                                                          >  1:233162/1‑150 (MQ=255)
                                        ggTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGc                                                                                                          >  2:1764336/1‑150 (MQ=255)
                                         gTAATTGCGTCTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTACTTGTGGTCGGCGTCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCt                                                                                                         <  1:3383445/150‑1 (MQ=255)
                                            aTTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggc                                                                                                      <  1:868891/150‑1 (MQ=255)
                                             ttGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcg                                                                                                     <  2:1994004/150‑1 (MQ=255)
                                                        ttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCg                                                                                          >  1:1669926/1‑150 (MQ=255)
                                                        ttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCg                                                                                          <  1:2220101/150‑1 (MQ=255)
                                                         tgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGc                                                                                         >  2:4360398/1‑150 (MQ=255)
                                                           aaCTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCa                                                                                       >  1:3744272/1‑150 (MQ=255)
                                                             ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGt                                                                                     >  1:3984767/1‑150 (MQ=255)
                                                             ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGt                                                                                     >  1:1461638/1‑150 (MQ=255)
                                                              tactTTTCCAGAGATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGtt                                                                                    <  2:3380636/150‑1 (MQ=255)
                                                                 ttttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAt                                                                                 >  2:2808809/1‑150 (MQ=255)
                                                                 ttttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAt                                                                                 >  2:4300655/1‑150 (MQ=255)
                                                                    tCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAc                                                                              <  2:2726315/150‑1 (MQ=255)
                                                                     ccTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAcc                                                                             >  2:3130443/1‑150 (MQ=255)
                                                                        gATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTg                                                                          >  1:228966/1‑150 (MQ=255)
                                                                         atatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGc                                                                         >  1:92618/1‑150 (MQ=255)
                                                                         atatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGc                                                                         >  1:4575445/1‑150 (MQ=255)
                                                                          tatGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCt                                                                        >  1:990390/1‑150 (MQ=255)
                                                                           atGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTc                                                                       >  2:3504034/1‑150 (MQ=255)
                                                                            tGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCg                                                                      <  2:2567887/150‑1 (MQ=255)
                                                                               cGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCg                                                                   <  1:2792922/150‑1 (MQ=255)
                                                                                 ttGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTg                                                                 >  1:182072/1‑150 (MQ=255)
                                                                                    gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                              <  1:3099069/150‑1 (MQ=255)
                                                                                     cAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTg                                                             >  2:667018/1‑150 (MQ=255)
                                                                                      aaaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGc                                                            <  1:2114918/150‑1 (MQ=255)
                                                                                      aaaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGc                                                            >  1:3623188/1‑150 (MQ=255)
                                                                                        aaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCg                                                          <  2:3205975/150‑1 (MQ=255)
                                                                                        aaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCg                                                          >  2:869114/1‑150 (MQ=255)
                                                                                          gCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                        <  1:2477463/150‑1 (MQ=255)
                                                                                          gCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCggg                                                        <  1:3845835/150‑1 (MQ=255)
                                                                                           cGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                                       <  2:1347386/150‑1 (MQ=255)
                                                                                            gTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                      <  1:324606/150‑1 (MQ=255)
                                                                                             tCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                     <  1:788751/150‑1 (MQ=255)
                                                                                             tCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                     <  1:2683359/150‑1 (MQ=255)
                                                                                              ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTg                                                    <  2:3999396/150‑1 (MQ=255)
                                                                                              ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTg                                                    <  1:3158764/150‑1 (MQ=255)
                                                                                                  gtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAt                                                <  2:2908927/150‑1 (MQ=255)
                                                                                                    ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                              <  1:1069087/150‑1 (MQ=255)
                                                                                                    ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                              >  2:3026010/1‑150 (MQ=255)
                                                                                                    ggACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCt                                              <  1:2389996/150‑1 (MQ=255)
                                                                                                     gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                             >  1:4018083/1‑150 (MQ=255)
                                                                                                     gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGTAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                             <  2:4660756/150‑1 (MQ=255)
                                                                                                      aCGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                            <  2:2231596/150‑1 (MQ=255)
                                                                                                      aCGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                            >  1:1261736/1‑150 (MQ=255)
                                                                                                      aCGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                            <  1:1574789/150‑1 (MQ=255)
                                                                                                       cggcggCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTa                                           <  2:233162/150‑1 (MQ=255)
                                                                                                       cggcggCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTa                                           >  2:1729456/1‑150 (MQ=255)
                                                                                                         gcggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAAc                                         <  2:1847956/150‑1 (MQ=255)
                                                                                                          cggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACg                                        >  1:2875775/1‑150 (MQ=255)
                                                                                                           ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                       <  1:4280701/150‑1 (MQ=255)
                                                                                                            gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                                      <  2:1528265/150‑1 (MQ=255)
                                                                                                            gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                                      >  2:2747851/1‑150 (MQ=255)
                                                                                                               aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                   <  2:1241472/150‑1 (MQ=255)
                                                                                                               aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                   <  1:4697813/150‑1 (MQ=255)
                                                                                                                  acacTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAcc                                >  1:2548029/1‑150 (MQ=255)
                                                                                                                    acTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                              >  1:3198087/1‑150 (MQ=255)
                                                                                                                    acTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                              <  2:48113/150‑1 (MQ=255)
                                                                                                                        gTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATAc                          >  2:2284930/1‑150 (MQ=255)
                                                                                                                        gTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATAc                          <  2:3802434/150‑1 (MQ=255)
                                                                                                                           aaGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt                       <  1:2645501/150‑1 (MQ=255)
                                                                                                                           aaGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt                       >  2:3318976/1‑150 (MQ=255)
                                                                                                                            aGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTg                      >  1:2204150/1‑150 (MQ=255)
                                                                                                                             gATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                     >  2:512767/1‑150 (MQ=255)
                                                                                                                             gATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGc                     >  1:152032/1‑150 (MQ=255)
                                                                                                                                  gTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGt                >  2:4564138/1‑150 (MQ=255)
                                                                                                                                   tGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTa               >  1:2016759/1‑150 (MQ=255)
                                                                                                                                    gCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAc              <  2:3609986/150‑1 (MQ=255)
                                                                                                                                     ccGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAcc             <  2:450465/150‑1 (MQ=255)
                                                                                                                                      cGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCa            >  2:4227492/1‑150 (MQ=255)
                                                                                                                                       gggTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCAt           >  1:3439079/1‑150 (MQ=255)
                                                                                                                                          ttACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt        >  2:3849655/1‑150 (MQ=255)
                                                                                                                                          ttACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt        >  1:3820921/1‑150 (MQ=255)
                                                                                                                                          ttACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCAGGCGGCAAAATTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt        >  1:1408974/1‑150 (MQ=255)
                                                                                                                                             cTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAAt     <  2:3103911/150‑1 (MQ=255)
                                                                                                                                              tGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATc    >  2:2441534/1‑150 (MQ=255)
                                                                                                                                               ggTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc   >  1:3344086/1‑150 (MQ=255)
                                                                                                                                                gTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCg  <  1:2924336/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCG  >  CP009273/3969563‑3969856

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: