Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,995,152 T→A 17.3% 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,1520TA17.3% 949.3 / 174.3 334E40V (GAA→GTA) fliARNA polymerase, sigma 28 (sigma F) factor
Reads supporting (aligned to +/- strand):  ref base T (126/150);  new base A (32/26);  total (158/176)
Fisher's exact test for biased strand distribution p-value = 1.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

CCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAT  >  CP009273/1995006‑1995299
                                                                                                                                                  |                                                                                                                                                   
cccAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTcc                                                                                                                                                   >  2:3674827/1‑149 (MQ=255)
cccAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTcc                                                                                                                                                   >  2:2583885/1‑149 (MQ=255)
cccAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTcc                                                                                                                                                   >  2:1338336/1‑149 (MQ=255)
       aCGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTcg                                                                                                                                            >  1:223443/1‑149 (MQ=255)
             gCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCa                                                                                                                                      <  1:7709/149‑1 (MQ=255)
               gAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGt                                                                                                                                    <  2:3258494/149‑1 (MQ=255)
                aaGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTc                                                                                                                                   >  1:2458924/1‑149 (MQ=255)
                  gTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGa                                                                                                                                 >  2:808160/1‑149 (MQ=255)
                      aTCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCt                                                                                                                             <  2:4682679/149‑1 (MQ=255)
                          aGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAg                                                                                                                         >  2:1476580/1‑149 (MQ=255)
                              tAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGgcgc                                                                                                                     >  2:2846550/1‑149 (MQ=255)
                               aGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCa                                                                                                                    >  2:778308/1‑149 (MQ=255)
                                gcgcCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCaa                                                                                                                   >  2:379727/1‑149 (MQ=255)
                                gcgcCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCaa                                                                                                                   >  2:1640785/1‑149 (MQ=255)
                                   ccACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGc                                                                                                                <  2:4250708/149‑1 (MQ=255)
                                   ccACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGc                                                                                                                <  1:919777/149‑1 (MQ=255)
                                      cGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTc                                                                                                             <  1:2672237/149‑1 (MQ=255)
                                       ggATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCg                                                                                                            >  2:4490682/1‑149 (MQ=255)
                                          tACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGa                                                                                                         >  1:1854849/1‑149 (MQ=255)
                                           aCGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGAc                                                                                                        >  1:2596951/1‑149 (MQ=255)
                                           aCGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGAc                                                                                                        <  2:527070/149‑1 (MQ=255)
                                              cTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCa                                                                                                     >  2:4560170/1‑149 (MQ=255)
                                               tGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCAc                                                                                                    <  2:2053176/149‑1 (MQ=255)
                                                   cTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGc                                                                                                >  1:554944/1‑149 (MQ=255)
                                                    tGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCg                                                                                               >  1:3554505/1‑149 (MQ=255)
                                                    tGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCg                                                                                               >  1:63694/1‑149 (MQ=255)
                                                        tAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGAc                                                                                           <  1:2210427/149‑1 (MQ=255)
                                                        tAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGAc                                                                                           >  1:4381369/1‑149 (MQ=255)
                                                          aGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACAt                                                                                         >  2:2675841/1‑149 (MQ=255)
                                                          aGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACAt                                                                                         <  1:3954648/149‑1 (MQ=255)
                                                            ttGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATaa                                                                                       <  2:4457460/149‑1 (MQ=255)
                                                              gTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACg                                                                                     >  2:3481910/1‑149 (MQ=255)
                                                                     ccGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCcac                                                                              <  1:1665312/149‑1 (MQ=255)
                                                                        ttCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGc                                                                           <  2:3011594/149‑1 (MQ=255)
                                                                           cTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAg                                                                        >  2:1978901/1‑149 (MQ=255)
                                                                            ttGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgt                                                                       >  1:1935134/1‑149 (MQ=255)
                                                                             tGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgtg                                                                      <  1:3674827/149‑1 (MQ=255)
                                                                             tGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgtg                                                                      <  1:1338336/149‑1 (MQ=255)
                                                                               tAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGtt                                                                    <  2:2924094/149‑1 (MQ=255)
                                                                                 gggCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTa                                                                  >  2:3595264/1‑149 (MQ=255)
                                                                                 gggCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTa                                                                  <  2:1082094/149‑1 (MQ=255)
                                                                                 gggCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTa                                                                  <  2:330591/149‑1 (MQ=255)
                                                                                      tCATAGCGTTCGACGGCATTAAGTAACCCAATGTCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCAt                                                             <  1:4617745/149‑1 (MQ=255)
                                                                                      tCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCAt                                                             >  2:2398339/1‑149 (MQ=255)
                                                                                         tAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGCCATAACGCTGCCACAGCGAGTGTTTATCCATTac                                                          >  1:2349720/1‑149 (MQ=255)
                                                                                         tAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTac                                                          <  1:67118/149‑1 (MQ=255)
                                                                                         tAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTac                                                          >  1:3091384/1‑149 (MQ=255)
                                                                                          aGCGTTCGACGGCAGTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTaca                                                         >  2:219127/1‑149 (MQ=255)
                                                                                           gCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTacac                                                        >  1:3703752/1‑149 (MQ=255)
                                                                                           gCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTacac                                                        <  1:2132381/149‑1 (MQ=255)
                                                                                              ttCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCtt                                                     <  1:2063474/149‑1 (MQ=255)
                                                                                               tCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTc                                                    <  2:2474121/149‑1 (MQ=255)
                                                                                               tCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTc                                                    <  2:2130039/149‑1 (MQ=255)
                                                                                                cGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCa                                                   <  1:3142163/149‑1 (MQ=255)
                                                                                                cGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACATTCa                                                   <  2:127373/149‑1 (MQ=255)
                                                                                                 gACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAg                                                  >  1:3849968/1‑149 (MQ=255)
                                                                                                 gACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAg                                                  >  2:18639/1‑149 (MQ=255)
                                                                                                   cGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCg                                                <  2:4261553/149‑1 (MQ=255)
                                                                                                    ggCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCgg                                               <  2:223958/149‑1 (MQ=255)
                                                                                                     gCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGt                                              >  2:3235297/1‑149 (MQ=255)
                                                                                                     gCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGt                                              <  2:377853/149‑1 (MQ=255)
                                                                                                      cATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGta                                             <  2:1707421/149‑1 (MQ=255)
                                                                                                      cATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGta                                             <  2:1392245/149‑1 (MQ=255)
                                                                                                       aTTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGtat                                            <  2:4176013/149‑1 (MQ=255)
                                                                                                        ttAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGtata                                           >  2:2751468/1‑149 (MQ=255)
                                                                                                        ttAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGtata                                           >  1:4558516/1‑149 (MQ=255)
                                                                                                          aaGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATaga                                         >  1:3565211/1‑149 (MQ=255)
                                                                                                           aGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATagag                                        >  2:3141546/1‑149 (MQ=255)
                                                                                                           aGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATagag                                        <  2:3975557/149‑1 (MQ=255)
                                                                                                           aGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATagag                                        <  2:3042878/149‑1 (MQ=255)
                                                                                                                cccAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAt                                   <  2:1204500/149‑1 (MQ=255)
                                                                                                                cccAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAt                                   >  1:1165331/1‑149 (MQ=255)
                                                                                                                 ccAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGACtt                                  >  2:1666551/1‑149 (MQ=255)
                                                                                                                 ccAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAtt                                  <  2:3920900/149‑1 (MQ=255)
                                                                                                                 ccAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAtt                                  <  2:1075286/149‑1 (MQ=255)
                                                                                                                  cAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                 <  1:198045/149‑1 (MQ=255)
                                                                                                                  cAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                 >  1:1115927/1‑149 (MQ=255)
                                                                                                                  cAATGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                 >  1:2912971/1‑149 (MQ=255)
                                                                                                                   aaTGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCa                                <  1:772033/149‑1 (MQ=255)
                                                                                                                   aaTGCCGCCCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCa                                >  1:1028592/1‑149 (MQ=255)
                                                                                                                      gccgccCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                             >  2:2302746/1‑149 (MQ=255)
                                                                                                                      gccgccCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                             <  2:4333827/149‑1 (MQ=255)
                                                                                                                      gccgccCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                             <  1:1250341/149‑1 (MQ=255)
                                                                                                                      gccgccCGCCTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                             <  1:2503538/149‑1 (MQ=255)
                                                                                                                         gcccgccTGTAGCAGATCGTCAAGTACCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATaa                          <  1:3476024/149‑1 (MQ=255)
                                                                                                                          cccgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATaaa                         <  1:4154932/149‑1 (MQ=255)
                                                                                                                           ccgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAc                        <  1:1069612/149‑1 (MQ=255)
                                                                                                                            cgccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                       <  1:1464242/149‑1 (MQ=255)
                                                                                                                             gccTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGGATAGAGTGAATTCACGATAAACAg                      >  2:3392901/1‑149 (MQ=255)
                                                                                                                                tGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                   <  2:4080735/149‑1 (MQ=255)
                                                                                                                                tGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                   <  1:688991/149‑1 (MQ=255)
                                                                                                                                tGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCCGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                   <  2:3850615/149‑1 (MQ=255)
                                                                                                                                 gTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCt                  <  1:3228727/149‑1 (MQ=255)
                                                                                                                                    gCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCg               <  2:635892/149‑1 (MQ=255)
                                                                                                                                     cAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt              <  2:2887469/149‑1 (MQ=255)
                                                                                                                                      aGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGtt             >  2:1896515/1‑149 (MQ=255)
                                                                                                                                       gATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTta            <  1:4244995/149‑1 (MQ=255)
                                                                                                                                       gATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTta            <  1:4126384/149‑1 (MQ=255)
                                                                                                                                       gATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTta            <  1:1759842/149‑1 (MQ=255)
                                                                                                                                         tCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtata          <  2:1030099/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCCAGTCACGGCTGCGAAGTTCATCCAGCATAGCGCCACGGATACGCTGCACTGCGTAAGTTGTAAATGCCGTTCCTTGTAGGGCGTCATAGCGTTCGACGGCATTAAGTAACCCAATGCCGCCCGCCTGTAGCAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAT  >  CP009273/1995006‑1995299

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: