Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,175,193 A→T 9.4% E36V (GAA→GTA)  cobB → deacetylase of acs and cheY, chemotaxis regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,175,1930AT 9.4% 868.8 / 65.5 277E36V (GAA→GTA) cobBdeacetylase of acs and cheY, chemotaxis regulator
Reads supporting (aligned to +/- strand):  ref base A (143/108);  new base T (4/22);  total (147/130)
Fisher's exact test for biased strand distribution p-value = 4.81e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 284 total aligned reads displayed.

GCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGATC  >  CP009273/1175044‑1175342
                                                                                                                                                     |                                                                                                                                                     
 cGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGaa                                                                                                                                                      <  2:2107983/150‑1 (MQ=255)
     aGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGAtg                                                                                                                                                  >  1:1739632/1‑150 (MQ=255)
        gTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTgg                                                                                                                                               >  2:241969/1‑150 (MQ=255)
          aTTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCa                                                                                                                                             <  1:507130/150‑1 (MQ=255)
          aTTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCa                                                                                                                                             >  2:1905971/1‑150 (MQ=255)
           tttttttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCaa                                                                                                                                            >  2:4202223/1‑150 (MQ=255)
               tttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGcc                                                                                                                                        <  1:1421194/150‑1 (MQ=255)
                ttCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCg                                                                                                                                       <  1:4153415/150‑1 (MQ=255)
                 tCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCgg                                                                                                                                      <  1:2174779/150‑1 (MQ=255)
                 tCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCgg                                                                                                                                      <  1:183184/150‑1 (MQ=255)
                   agagATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGaa                                                                                                                                    >  1:895229/1‑150 (MQ=255)
                   agagATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGaa                                                                                                                                    <  1:1668013/150‑1 (MQ=255)
                   agagATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGaa                                                                                                                                    <  1:3807593/150‑1 (MQ=255)
                     agaTAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAgg                                                                                                                                  >  2:4175356/1‑150 (MQ=255)
                       aTAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGtt                                                                                                                                <  2:1046632/150‑1 (MQ=255)
                       aTAAAGAGGTGCCGGAAGCAATCGAAAAACCAAGAGTACTCGTACTGACCGGGGCAGGAATTTCTGCGGATTCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGtt                                                                                                                                <  2:4189209/150‑1 (MQ=255)
                         aaaGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCACCGCCGGAAGGTTTc                                                                                                                              >  2:2787528/1‑150 (MQ=255)
                         aaaGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTc                                                                                                                              <  1:2593823/150‑1 (MQ=255)
                          aaGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCg                                                                                                                             <  2:4392723/150‑1 (MQ=255)
                          aaGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCg                                                                                                                             >  2:4634992/1‑150 (MQ=255)
                            gtggtgCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGAt                                                                                                                           <  2:3975176/150‑1 (MQ=255)
                            gtggtgCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGAt                                                                                                                           <  2:4534629/150‑1 (MQ=255)
                                tgCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATcgcg                                                                                                                       <  1:881424/150‑1 (MQ=255)
                                  cctgAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGAt                                                                                                                     <  2:2739192/147‑1 (MQ=255)
                                  ccGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGAt                                                                                                                     >  1:1023586/1‑150 (MQ=255)
                                   cGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATc                                                                                                                    <  2:2370646/150‑1 (MQ=255)
                                    ggAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATcc                                                                                                                   >  1:3161430/1‑150 (MQ=255)
                                    ggAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATcc                                                                                                                   >  1:692699/1‑150 (MQ=255)
                                    ggAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATcc                                                                                                                   >  2:3093028/1‑150 (MQ=255)
                                        tcAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGACTCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGaa                                                                                                               >  2:746174/2‑150 (MQ=255)
                                        gCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGaa                                                                                                               <  2:3274488/150‑1 (MQ=255)
                                         cAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAAc                                                                                                              >  2:298601/1‑150 (MQ=255)
                                         cAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAAc                                                                                                              >  2:677961/1‑150 (MQ=255)
                                          aaTGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACt                                                                                                             >  2:836638/1‑150 (MQ=255)
                                          aaTGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACt                                                                                                             <  2:3754407/150‑1 (MQ=255)
                                           aTGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTg                                                                                                            <  1:4159197/150‑1 (MQ=255)
                                              gAAAAACCAAGAGTACTCTTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTTGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTg                                                                                                         >  2:1799675/1‑150 (MQ=255)
                                                aaaaCCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCa                                                                                                       >  1:3724286/1‑150 (MQ=255)
                                                aaaaCCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCa                                                                                                       >  2:4096837/1‑150 (MQ=255)
                                                 aaaCCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCaa                                                                                                      <  1:3592527/150‑1 (MQ=255)
                                                  aaCCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAg                                                                                                     >  1:4381716/1‑150 (MQ=255)
                                                  aaCCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAg                                                                                                     >  2:697011/1‑150 (MQ=255)
                                                     cAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGt                                                                                                  <  2:2252335/150‑1 (MQ=255)
                                                     cAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGt                                                                                                  >  1:1515766/1‑150 (MQ=255)
                                                      aaGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGtt                                                                                                 >  2:444731/1‑150 (MQ=255)
                                                      aaGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGtt                                                                                                 <  2:1707637/150‑1 (MQ=255)
                                                       cgagTTCTCGTCCTGACAGGGGCAGGACTTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGttt                                                                                                <  1:3221202/149‑1 (MQ=255)
                                                       agagTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGttt                                                                                                >  2:2912560/1‑150 (MQ=255)
                                                        gagTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGtttt                                                                                               >  1:960825/1‑150 (MQ=255)
                                                        gagTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGtttt                                                                                               >  1:3379952/1‑150 (MQ=255)
                                                           tACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTtata                                                                                            <  2:3391819/150‑1 (MQ=255)
                                                              tCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACg                                                                                         <  2:4329253/150‑1 (MQ=255)
                                                                gTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGcc                                                                                       >  2:2958490/1‑150 (MQ=255)
                                                                gTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGcc                                                                                       <  1:2041174/150‑1 (MQ=255)
                                                                  aCTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcg                                                                                     <  1:1201193/150‑1 (MQ=255)
                                                                   cTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgt                                                                                    <  2:53683/150‑1 (MQ=255)
                                                                   cTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgt                                                                                    >  2:2538124/1‑150 (MQ=255)
                                                                    tGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgtc                                                                                   >  1:2149741/1‑150 (MQ=255)
                                                                     gACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgtcg                                                                                  <  2:2840874/150‑1 (MQ=255)
                                                                     gACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCcgtcg                                                                                  >  2:2623617/1‑150 (MQ=255)
                                                                       cAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCgtcgtc                                                                                >  2:1399567/1‑150 (MQ=255)
                                                                             cAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGc                                                                          >  2:3652616/1‑150 (MQ=255)
                                                                              aGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGCAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCt                                                                         >  2:2517582/1‑150 (MQ=255)
                                                                                          ggAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGa                                                             >  1:2557349/1‑150 (MQ=255)
                                                                                               cAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTc                                                        <  1:4361492/150‑1 (MQ=255)
                                                                                                    aTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCt                                                   >  1:1146800/1‑150 (MQ=255)
                                                                                                    aTTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAAAAATTCAGCCt                                                   >  1:2469096/1‑150 (MQ=255)
                                                                                                       cGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGGTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAc                                                >  2:4514354/1‑150 (MQ=255)
                                                                                                        gTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAAcg                                               >  1:3465924/1‑150 (MQ=255)
                                                                                                             tttCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAATCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCATCCATAAATTCAGCCTAACGCcgcg                                          <  2:297893/150‑1 (MQ=255)
                                                                                                               tCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCa                                        >  2:1278107/1‑150 (MQ=255)
                                                                                                                  cgccgcAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCt                                     >  2:435256/1‑150 (MQ=255)
                                                                                                                   gccgcAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCtt                                    <  2:3835044/150‑1 (MQ=255)
                                                                                                                   gccgcAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCtt                                    >  2:3900477/1‑150 (MQ=255)
                                                                                                                     cgcAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTgc                                  >  1:4625489/1‑150 (MQ=255)
                                                                                                                       cAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTgcgc                                >  1:26888/1‑150 (MQ=255)
                                                                                                                        aGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCt                               <  1:2623617/150‑1 (MQ=255)
                                                                                                                              ccTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTaa                         >  2:3565283/1‑150 (MQ=255)
                                                                                                                              ccTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTaa                         >  2:3673046/1‑150 (MQ=255)
                                                                                                                                 gtgGGAAGAACATCGGGTTGTAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACt                      <  1:3900004/150‑1 (MQ=255)
                                                                                                                                   gggAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGc                    <  1:1745207/150‑1 (MQ=255)
                                                                                                                                   gggAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGc                    >  1:2487973/1‑150 (MQ=255)
                                                                                                                                    ggAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCa                   >  1:4068748/1‑150 (MQ=255)
                                                                                                                                    ggAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCa                   >  1:2780034/1‑150 (MQ=255)
                                                                                                                                        gaaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGAc               >  2:492130/1‑150 (MQ=255)
                                                                                                                                        gaaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGAc               >  2:4465151/1‑150 (MQ=255)
                                                                                                                                        gaaCATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGAGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGAc               >  1:2313379/1‑150 (MQ=255)
                                                                                                                                           cATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGcc            >  2:2444028/1‑150 (MQ=255)
                                                                                                                                           cATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCAGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGcc            >  1:2205485/1‑150 (MQ=255)
                                                                                                                                            aTCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGccc           >  1:886082/1‑150 (MQ=255)
                                                                                                                                             tCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCt          <  1:1867239/150‑1 (MQ=255)
                                                                                                                                             tCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCt          <  2:2730991/150‑1 (MQ=255)
                                                                                                                                              cGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTc         >  2:1562126/1‑150 (MQ=255)
                                                                                                                                              cGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTc         >  1:2955537/1‑150 (MQ=255)
                                                                                                                                               gggTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCg        <  2:895233/150‑1 (MQ=255)
                                                                                                                                                 gTTGAAGATGTGGCATCGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGc      <  2:3573914/150‑1 (MQ=255)
                                                                                                                                                 gTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGc      <  1:2420079/150‑1 (MQ=255)
                                                                                                                                                     aaGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGATc  <  1:601621/150‑1 (MQ=255)
                                                                                                                                                     aaGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGATc  <  1:3944789/150‑1 (MQ=255)
                                                                                                                                                     aaGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGATc  <  1:3118734/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTCAGCGTATTTTTTTCAGAGATAAAGTGGTGCCGGAAGCAATGGAAAAACCAAGAGTACTCGTACTGACAGGGGCAGGAATTTCTGCGGAATCAGGTATTCGTACCTTTCGCGCCGCAGATGGCCTGTGGGAAGAACATCGGGTTGAAGATGTGGCAACGCCGGAAGGTTTCGATCGCGATCCTGAACTGGTGCAAGCGTTTTATAACGCCCGTCGTCGACAGCTGCAGCAGCCAGAAATTCAGCCTAACGCCGCGCATCTTGCGCTGGCTAAACTGCAAGACGCCCTCGGCGATC  >  CP009273/1175044‑1175342

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: