Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,995,152 T→A 12.4% E40V (GAA→GTA)  fliA ← RNA polymerase, sigma 28 (sigma F) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,995,1520TA12.4% 1142.7 / 130.5 348E40V (GAA→GTA) fliARNA polymerase, sigma 28 (sigma F) factor
Reads supporting (aligned to +/- strand):  ref base T (120/185);  new base A (19/24);  total (139/209)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 354 total aligned reads displayed.

CCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATC  >  CP009273/1995006‑1995300
                                                                                                                                                  |                                                                                                                                                    
   aGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACg                                                                                                                                                 >  1:1773050/1‑149 (MQ=255)
        cGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTcgc                                                                                                                                            >  1:4185015/1‑149 (MQ=255)
         ggCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTcgcg                                                                                                                                           <  1:1556231/149‑1 (MQ=255)
          gCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCgg                                                                                                                                          >  1:1637196/1‑149 (MQ=255)
          gCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGACGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCgg                                                                                                                                          >  1:3215981/1‑149 (MQ=255)
           cTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCggg                                                                                                                                         <  1:1583716/149‑1 (MQ=255)
            tGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGc                                                                                                                                        <  2:3573874/149‑1 (MQ=255)
                    tCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAAc                                                                                                                                <  2:4271910/149‑1 (MQ=255)
                    tCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAAc                                                                                                                                <  1:3028093/149‑1 (MQ=255)
                     cATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAAcc                                                                                                                               >  2:4283233/1‑149 (MQ=255)
                      aTCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCt                                                                                                                              <  1:3684120/149‑1 (MQ=255)
                      aTCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCt                                                                                                                              <  1:4455440/149‑1 (MQ=255)
                        ccAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGc                                                                                                                            >  1:4002699/1‑149 (MQ=255)
                        ccAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGc                                                                                                                            >  1:5168127/1‑149 (MQ=255)
                          aGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAg                                                                                                                          >  2:2261762/1‑149 (MQ=255)
                           gCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAgg                                                                                                                         <  1:5208335/149‑1 (MQ=255)
                               aGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCa                                                                                                                     >  1:2660461/1‑149 (MQ=255)
                               aGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCa                                                                                                                     >  2:2745663/1‑149 (MQ=255)
                                gctccACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCaa                                                                                                                    <  2:5175031/146‑1 (MQ=255)
                                gcgcCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCaa                                                                                                                    >  1:2845686/1‑149 (MQ=255)
                                gcgcCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCaa                                                                                                                    <  2:1181813/149‑1 (MQ=255)
                                    cACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCt                                                                                                                <  1:4380336/149‑1 (MQ=255)
                                       ggATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCg                                                                                                             <  1:991649/149‑1 (MQ=255)
                                         aTACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTg                                                                                                           >  2:260563/1‑149 (MQ=255)
                                          tACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGa                                                                                                          >  1:2690251/1‑149 (MQ=255)
                                          tACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGa                                                                                                          <  1:1893575/149‑1 (MQ=255)
                                               tGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCAc                                                                                                     >  1:4626265/1‑149 (MQ=255)
                                               tGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCAc                                                                                                     <  1:4209349/149‑1 (MQ=255)
                                                  aCTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAg                                                                                                  >  2:2197218/1‑149 (MQ=255)
                                                  aCTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAg                                                                                                  >  2:699574/1‑149 (MQ=255)
                                                    tGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCg                                                                                                <  2:775409/149‑1 (MQ=255)
                                                        tAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGAc                                                                                            >  2:4058145/1‑149 (MQ=255)
                                                         aaGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACa                                                                                           <  2:3601798/149‑1 (MQ=255)
                                                         aaGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACa                                                                                           <  2:2420046/149‑1 (MQ=255)
                                                          aGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACAt                                                                                          <  1:563006/149‑1 (MQ=255)
                                                           gTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATa                                                                                         <  1:4280785/149‑1 (MQ=255)
                                                           gTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATa                                                                                         <  1:585562/149‑1 (MQ=255)
                                                            ttGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATaa                                                                                        <  2:1851041/149‑1 (MQ=255)
                                                            ttGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATaa                                                                                        <  2:1851088/149‑1 (MQ=255)
                                                            ttGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATaa                                                                                        <  2:4467246/149‑1 (MQ=255)
                                                               tAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGc                                                                                     >  2:3461680/1‑149 (MQ=255)
                                                                aaaTGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCt                                                                                    >  2:4846908/1‑149 (MQ=255)
                                                                aaaTGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCt                                                                                    <  2:3705706/149‑1 (MQ=255)
                                                                  aTGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGc                                                                                  <  2:4798771/149‑1 (MQ=255)
                                                                   tGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGAGCGTCACGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCACTGCTTCGTGACGCACCAGCGGGACATAACGCTGcc                                                                                 >  2:4942358/1‑149 (MQ=255)
                                                                     ccGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcac                                                                               <  1:1262074/149‑1 (MQ=255)
                                                                     ccGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcac                                                                               >  2:1734278/1‑149 (MQ=255)
                                                                      cGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcaca                                                                              <  1:4836958/149‑1 (MQ=255)
                                                                        ttCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGc                                                                            >  1:2543627/1‑149 (MQ=255)
                                                                         tCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCg                                                                           <  1:605594/149‑1 (MQ=255)
                                                                          ccTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGa                                                                          >  2:3064443/1‑149 (MQ=255)
                                                                          ccTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGa                                                                          <  2:654608/149‑1 (MQ=255)
                                                                            ttGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgt                                                                        <  1:2735659/149‑1 (MQ=255)
                                                                            ttGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgt                                                                        <  1:3573636/149‑1 (MQ=255)
                                                                               tAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGtt                                                                     >  1:2601117/1‑149 (MQ=255)
                                                                                 gggCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTa                                                                   <  2:3648454/149‑1 (MQ=255)
                                                                                 gggCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTa                                                                   <  1:2078801/149‑1 (MQ=255)
                                                                                 gggCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTa                                                                   <  1:3865328/149‑1 (MQ=255)
                                                                                    cGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATcc                                                                <  2:5231042/149‑1 (MQ=255)
                                                                                     gTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCa                                                               <  1:3139884/149‑1 (MQ=255)
                                                                                      tCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCAt                                                              >  1:2877427/1‑149 (MQ=255)
                                                                                              ttCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCtt                                                      <  1:4306303/149‑1 (MQ=255)
                                                                                               tCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTc                                                     <  2:3485372/149‑1 (MQ=255)
                                                                                                   cGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCg                                                 <  2:4766266/149‑1 (MQ=255)
                                                                                                   cGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCg                                                 <  1:1738672/149‑1 (MQ=255)
                                                                                                    ggCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCgg                                                <  1:727631/149‑1 (MQ=255)
                                                                                                    ggCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCgg                                                >  2:2282288/1‑149 (MQ=255)
                                                                                                    ggCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCgg                                                <  2:4420856/149‑1 (MQ=255)
                                                                                                     gCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGt                                               >  2:2992962/1‑149 (MQ=255)
                                                                                                      cATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGta                                              <  2:3907893/149‑1 (MQ=255)
                                                                                                      cATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGta                                              >  1:3400410/1‑149 (MQ=255)
                                                                                                      cATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGta                                              <  1:1689945/149‑1 (MQ=255)
                                                                                                       aTTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGtat                                             >  2:1703941/1‑149 (MQ=255)
                                                                                                        ttAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGtata                                            >  2:438183/1‑149 (MQ=255)
                                                                                                        ttAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGtata                                            <  1:3716658/149‑1 (MQ=255)
                                                                                                           aGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATagag                                         <  1:2181703/149‑1 (MQ=255)
                                                                                                              aaCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGa                                      >  1:1374839/1‑149 (MQ=255)
                                                                                                               aCCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                     <  2:2883889/149‑1 (MQ=255)
                                                                                                                cccAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAt                                    <  1:4775848/149‑1 (MQ=255)
                                                                                                                 ccAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAtt                                   <  1:4507219/149‑1 (MQ=255)
                                                                                                                    aTGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                >  1:1039346/1‑149 (MQ=255)
                                                                                                                    aTGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                >  1:834537/1‑149 (MQ=255)
                                                                                                                      gccgccCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                              <  1:3492184/149‑1 (MQ=255)
                                                                                                                      gccgccCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                              <  1:5004528/149‑1 (MQ=255)
                                                                                                                      gccgccCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                              <  1:2430448/149‑1 (MQ=255)
                                                                                                                         gcccgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATaa                           >  1:4736179/1‑149 (MQ=255)
                                                                                                                          cccgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATaaa                          <  1:3351541/149‑1 (MQ=255)
                                                                                                                           ccgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAc                         >  1:2217266/1‑149 (MQ=255)
                                                                                                                           ccgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAc                         <  1:4717385/149‑1 (MQ=255)
                                                                                                                            cgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                        <  2:5094908/149‑1 (MQ=255)
                                                                                                                             gccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAg                       <  1:1350751/149‑1 (MQ=255)
                                                                                                                             gccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAg                       >  2:1547886/1‑149 (MQ=255)
                                                                                                                               cTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGGATAGAGTGAATTCACGATAAACAGcc                     >  2:1388180/1‑149 (MQ=255)
                                                                                                                                 gTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCt                   >  2:3946139/1‑149 (MQ=255)
                                                                                                                                  tAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTg                  <  2:2709915/149‑1 (MQ=255)
                                                                                                                                   aGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGc                 <  2:1776530/149‑1 (MQ=255)
                                                                                                                                               aGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGtt     <  1:4181589/149‑1 (MQ=255)
                                                                                                                                                gTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCTGCGGGACATAACGCTGCCCCTGCGAGTGTTTATCCATTACAGCTTCAGCGGGAGAGAGTGAATTCACGATAATCAGCCCTGCGTTATATGAGTTa    >  1:2315243/1‑149 (MQ=255)
                                                                                                                                                 ttCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt   <  2:3186739/149‑1 (MQ=255)
                                                                                                                                                 ttCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt   >  2:1312225/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATC  >  CP009273/1995006‑1995300

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: