Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,325,982 G→T 86.6% T226T (ACG→ACT topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,325,9820GT86.6% 540.3 / 85.6 241T226T (ACG→ACTtopADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (15/18);  new base T (88/115);  total (108/133)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 241 total aligned reads displayed.

TGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAC  >  CP009273/1325837‑1326127
                                                                                                                                                 |                                                                                                                                                 
tGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAcc                                                                                                                                                <  1:3526027/149‑1 (MQ=255)
tGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTAcc                                                                                                                                                >  1:3529227/1‑149 (MQ=255)
 ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCa                                                                                                                                               >  2:2546175/1‑149 (MQ=255)
 ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCa                                                                                                                                               <  1:3455281/149‑1 (MQ=255)
 ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCa                                                                                                                                               >  2:2184814/1‑149 (MQ=255)
 ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCa                                                                                                                                               >  2:2184662/1‑149 (MQ=255)
 ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCa                                                                                                                                               <  2:4341651/149‑1 (MQ=255)
  gTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCAc                                                                                                                                              >  1:2873302/1‑149 (MQ=255)
     tCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGGGAGCGTGAAATTAAATCGTTCGTGCCGGAAGAGTTCGGGGAAGTCGATGCCAGCACGACCAcccc                                                                                                                                           >  1:1616705/1‑149 (MQ=255)
       gccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACAACCACGCCaa                                                                                                                                         >  1:2884626/1‑148 (MQ=255)
        ccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATc                                                                                                                                        >  2:1540394/1‑149 (MQ=255)
        ccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATc                                                                                                                                        <  2:1126798/149‑1 (MQ=255)
         cgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGGGAAATTAAAGCGTTAGTGCCGGAAGAGTTCTGGGAAGTCGATGCCCGCACAACCACGCCATCt                                                                                                                                       >  1:4423815/1‑149 (MQ=255)
           ctgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTgg                                                                                                                                     <  2:649635/149‑1 (MQ=255)
           ctgctATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTgg                                                                                                                                     <  2:3708542/149‑1 (MQ=255)
                 tGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGc                                                                                                                               <  2:4016743/149‑1 (MQ=255)
                   gAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGt                                                                                                                             <  2:4606425/149‑1 (MQ=255)
                    aaaaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGtt                                                                                                                            >  2:1145028/1‑149 (MQ=255)
                     aaaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTTGCGGTTCGCCTGGTGGTCGATCGTGATCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATTCCAGCACTACCACGCCATCATGTGAAGAGTTg                                                                                                                           >  1:4153668/1‑149 (MQ=255)
                     aaaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTg                                                                                                                           >  1:1763894/1‑149 (MQ=255)
                      aaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTgg                                                                                                                          >  2:3925223/1‑149 (MQ=255)
                      aaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTgg                                                                                                                          <  1:1935227/149‑1 (MQ=255)
                      aaaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTgg                                                                                                                          >  2:3086851/1‑149 (MQ=255)
                        aGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTTAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCg                                                                                                                        >  1:749656/1‑149 (MQ=255)
                           tCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTa                                                                                                                     <  1:4274808/149‑1 (MQ=255)
                             gCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGGTCGCCGGGGGGTCGAGCGTGAGCGGGAAAGTAAAGCGGTCGTGCCGGAAGAGTTCTGGGACGTCGATGCCCCCACGACCCCGCCCACTGGTGAAGCGTTGGCGGTACa                                                                                                                   >  1:1243725/1‑149 (MQ=255)
                              cTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGGGAAATTAAAGCGTTCGTGCCGGGAGAGTTCTGGGCAGTCGATGCCAGCACAACCACGCCCTCTGGTGAAGCGGTGGCGTTACAg                                                                                                                  >  1:1386075/1‑149 (MQ=255)
                               tCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAgg                                                                                                                 <  2:4174054/149‑1 (MQ=255)
                                cGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGt                                                                                                                <  1:3251045/149‑1 (MQ=255)
                                 gTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTg                                                                                                               <  2:3225060/149‑1 (MQ=255)
                                  tGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGa                                                                                                              <  1:2958163/149‑1 (MQ=255)
                                       tGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcat                                                                                                         <  1:4122318/149‑1 (MQ=255)
                                       tGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcat                                                                                                         >  2:2164236/1‑149 (MQ=255)
                                         tCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACtcatca                                                                                                       <  2:3779602/149‑1 (MQ=255)
                                           tGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGa                                                                                                     >  1:4691536/1‑149 (MQ=255)
                                            gCCGGTCGTGTGCAGTCGGTGGCGGATCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGAGTTCGTGCCGGAAGAGTTCAGGGAAGTCGATTCCAGCACAACCACGCCATCTGGAGAAGCGTAGGAGTTACAGGTTAATCAACAGaa                                                                                                    >  2:1576900/1‑149 (MQ=255)
                                                gTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGAc                                                                                                >  1:4332350/1‑149 (MQ=255)
                                                 tCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACa                                                                                               <  2:4096097/149‑1 (MQ=255)
                                                    tgtgCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAAc                                                                                            >  2:4546427/1‑149 (MQ=255)
                                                      tgCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCg                                                                                          <  1:2354857/149‑1 (MQ=255)
                                                         aGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtc                                                                                       <  2:1309122/149‑1 (MQ=255)
                                                          gTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtcc                                                                                      <  1:4520798/149‑1 (MQ=255)
                                                           tCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccg                                                                                     >  2:713012/1‑149 (MQ=255)
                                                              gTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTtccgtcc                                                                                  >  2:2016212/1‑149 (MQ=255)
                                                                ggCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCgg                                                                                <  1:4433056/149‑1 (MQ=255)
                                                                  cGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTc                                                                              <  2:1283681/149‑1 (MQ=255)
                                                                  cGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTc                                                                              <  2:1055570/149‑1 (MQ=255)
                                                                  cGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTc                                                                              <  2:4089783/149‑1 (MQ=255)
                                                                  cGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTc                                                                              <  1:3381661/149‑1 (MQ=255)
                                                                       cGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTcaacaa                                                                         >  2:280333/1‑149 (MQ=255)
                                                                         ccTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAg                                                                       >  1:1011385/1‑149 (MQ=255)
                                                                          cTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGa                                                                      <  1:2484876/149‑1 (MQ=255)
                                                                            ggtggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAAc                                                                    <  2:1531629/149‑1 (MQ=255)
                                                                               ggtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACaaa                                                                 >  2:55865/1‑149 (MQ=255)
                                                                                gtCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGACGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAAc                                                                >  2:1461201/1‑149 (MQ=255)
                                                                                    agcgtgagcgtgAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAg                                                            <  1:2027056/149‑1 (MQ=255)
                                                                                    agcgtgagcgtgAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAg                                                            >  1:4094405/1‑149 (MQ=255)
                                                                                     gcgtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAgg                                                           >  1:4334681/1‑149 (MQ=255)
                                                                                          agcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCg                                                      <  2:285736/149‑1 (MQ=255)
                                                                                           gcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCgg                                                     <  1:1583038/149‑1 (MQ=255)
                                                                                            cgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGt                                                    <  1:2463495/149‑1 (MQ=255)
                                                                                            cgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGt                                                    <  2:1335481/149‑1 (MQ=255)
                                                                                              tgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTaa                                                  <  1:3853655/149‑1 (MQ=255)
                                                                                              tgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTaa                                                  <  2:4228648/149‑1 (MQ=255)
                                                                                               taaaTTAAAGCGTTCGTGCCTGAAGAGTTCTGGGAAGTCGATGCCAGCACTAACACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAg                                                 <  1:429576/148‑1 (MQ=255)
                                                                                                aaaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGt                                                >  1:3030102/1‑149 (MQ=255)
                                                                                                 aaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTc                                               <  1:4144006/149‑1 (MQ=255)
                                                                                                 aaTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTc                                               >  2:2004100/1‑149 (MQ=255)
                                                                                                  aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTct                                              >  2:3022193/1‑149 (MQ=255)
                                                                                                   ttAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctg                                             >  1:2745853/1‑149 (MQ=255)
                                                                                                    tAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgc                                            <  1:2014625/149‑1 (MQ=255)
                                                                                                     aaaGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgct                                           >  1:29079/1‑149 (MQ=255)
                                                                                                      aaGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTctgctg                                          >  2:2320097/1‑149 (MQ=255)
                                                                                                          gTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaa                                      >  2:1691396/1‑149 (MQ=255)
                                                                                                           ttCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGaaaa                                     >  1:212000/1‑149 (MQ=255)
                                                                                                               tGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcg                                 <  2:2761867/149‑1 (MQ=255)
                                                                                                                gCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcgc                                <  2:4087575/149‑1 (MQ=255)
                                                                                                                  cGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcgcgc                              <  2:4661266/149‑1 (MQ=255)
                                                                                                                  cGGAAGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcgcgc                              <  1:421301/149‑1 (MQ=255)
                                                                                                                     aaGAGTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTAc                           <  1:4500202/149‑1 (MQ=255)
                                                                                                                     aaGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTAc                           <  1:863435/149‑1 (MQ=255)
                                                                                                                      agagTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACa                          <  2:3297930/149‑1 (MQ=255)
                                                                                                                       gagTTCTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAg                         >  2:208382/1‑149 (MQ=255)
                                                                                                                            cTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGc                    <  2:1939589/149‑1 (MQ=255)
                                                                                                                            cTGGGAAGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGc                    <  1:751961/149‑1 (MQ=255)
                                                                                                                              gggAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTg                  <  2:3479430/149‑1 (MQ=255)
                                                                                                                                  aGTCGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAAc              >  1:1808252/1‑149 (MQ=255)
                                                                                                                                  aGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAAc              >  1:1163839/1‑149 (MQ=255)
                                                                                                                                  aGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAAc              >  1:1164052/1‑149 (MQ=255)
                                                                                                                                   gTCGATGCCAGCACTACCACGCCATCTGGCGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACg             >  1:64392/1‑149 (MQ=255)
                                                                                                                                     cGATGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTg           <  1:2250096/149‑1 (MQ=255)
                                                                                                                                       aTGCCAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGaa         >  1:1294393/1‑149 (MQ=255)
                                                                                                                                          ccAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGAc      >  1:2578853/1‑149 (MQ=255)
                                                                                                                                           cAGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACa     <  2:1829264/149‑1 (MQ=255)
                                                                                                                                            aGCACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACaa    >  1:4672025/1‑149 (MQ=255)
                                                                                                                                              cACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  2:2114427/149‑1 (MQ=255)
                                                                                                                                              cACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  2:1813445/149‑1 (MQ=255)
                                                                                                                                              cACTACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  1:972122/149‑1 (MQ=255)
                                                                                                                                              cACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  2:1601245/149‑1 (MQ=255)
                                                                                                                                              cACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAc  <  2:1485609/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAC  >  CP009273/1325837‑1326127

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: