Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,741 G→C 73.6% A168A (GCG→GCC talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,7410GC73.6% 543.9 / 230.7 265A168A (GCG→GCCtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base C (111/84);  total (146/119)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-07

Only 100 of 265 total aligned reads displayed.

CGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGT  >  CP009273/8597‑8888
                                                                                                                                                |                                                                                                                                                   
 gATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCg                                                                                                                                               <  1:790979/150‑1 (MQ=255)
      tGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCgt                                                                                                                                          >  1:691827/1‑150 (MQ=255)
       ggcatcaGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCgtc                                                                                                                                         >  1:121822/7‑149 (MQ=255)
               cAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGAt                                                                                                                                 >  2:304496/1‑150 (MQ=255)
                    aTCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGc                                                                                                                            <  2:835260/150‑1 (MQ=255)
                     tCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGcc                                                                                                                           <  2:26426/150‑1 (MQ=255)
                      cgcaTTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCg                                                                                                                          >  2:843451/4‑150 (MQ=255)
                      cGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCg                                                                                                                          >  1:634055/1‑150 (MQ=255)
                          ttCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGttcg                                                                                                                      >  2:490149/1‑148 (MQ=255)
                              gATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTgg                                                                                                                  >  2:510435/1‑150 (MQ=255)
                               aTCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCcgttcgtt                                                                                                                    <  2:299513/147‑1 (MQ=255)
                                tCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGcc                                                                                                                <  1:495758/150‑1 (MQ=255)
                                tCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGcc                                                                                                                <  1:311894/150‑1 (MQ=255)
                                  aaaCTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCg                                                                                                               >  1:240259/1‑149 (MQ=255)
                                    aCTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTAt                                                                                                            >  1:220416/1‑150 (MQ=255)
                                      tGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTc                                                                                                          >  2:643372/1‑150 (MQ=255)
                                      tGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTAtcc                                                                                                          >  2:812848/1‑148 (MQ=255)
                                        gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                                                                                                        <  1:510435/144‑2 (MQ=255)
                                        gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCtc                                                                                                        <  1:441445/144‑2 (MQ=255)
                                        gcctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCt                                                                                                         <  1:668993/143‑1 (MQ=255)
                                        gCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCtt                                                                                                        <  2:756412/150‑1 (MQ=255)
                                          ctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGa                                                                                                      >  2:620653/5‑150 (MQ=255)
                                          ctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGa                                                                                                      >  1:483897/5‑150 (MQ=255)
                                          ctccaCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTAGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGa                                                                                                      >  2:255125/5‑150 (MQ=255)
                                             tACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTg                                                                                                   >  1:271913/1‑150 (MQ=255)
                                             caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAc                                                                                                     <  1:141530/147‑1 (MQ=255)
                                             caCCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGAc                                                                                                     <  1:725985/147‑1 (MQ=255)
                                               ccTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTgg                                                                                                  >  2:537955/1‑149 (MQ=255)
                                                cTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTAGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTa                                                                                                >  2:661418/1‑150 (MQ=255)
                                                cTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTa                                                                                                >  2:10781/1‑150 (MQ=255)
                                                  ggCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTAc                                                                                               >  1:386805/1‑149 (MQ=255)
                                                   gCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACa                                                                                              <  1:689551/149‑1 (MQ=255)
                                                    cAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACaaa                                                                                            <  1:109455/150‑1 (MQ=255)
                                                     aGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAg                                                                                           <  1:284679/150‑1 (MQ=255)
                                                     aGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAg                                                                                           >  1:771754/1‑150 (MQ=255)
                                                         tATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGaa                                                                                       <  2:189461/150‑1 (MQ=255)
                                                           tCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATa                                                                                     <  1:857424/150‑1 (MQ=255)
                                                             cGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATAcc                                                                                   <  2:437692/150‑1 (MQ=255)
                                                                 ctgcAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGAt                                                                                <  2:334964/149‑1 (MQ=255)
                                                                 ccgctgAACAGCTGGAAAAAGAAGGCATTAACTGTACCCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGaca                                                                               <  2:405/145‑3 (MQ=255)
                                                                 ccgctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGac                                                                                >  1:343662/6‑148 (MQ=255)
                                                                   gctgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGac                                                                                >  1:526804/4‑146 (MQ=255)
                                                                    cAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaa                                                                              <  1:293735/148‑1 (MQ=255)
                                                                     tgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACaaga                                                                            >  2:670078/2‑149 (MQ=255)
                                                                     tgAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACaaga                                                                            >  2:315424/2‑149 (MQ=255)
                                                                      gAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaagaaa                                                                          <  1:181502/150‑1 (MQ=255)
                                                                        aCAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATaagaaaga                                                                        <  2:725580/150‑1 (MQ=255)
                                                                         cAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAagag                                                                       <  2:681745/150‑1 (MQ=255)
                                                                          aGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGt                                                                      <  1:851569/150‑1 (MQ=255)
                                                                          aGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGt                                                                      <  2:601821/150‑1 (MQ=255)
                                                                          aGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGt                                                                      <  1:620653/150‑1 (MQ=255)
                                                                           gCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTa                                                                     <  1:573653/150‑1 (MQ=255)
                                                                            cTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTAc                                                                    <  2:526804/150‑1 (MQ=255)
                                                                            cTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTAc                                                                    <  2:390341/150‑1 (MQ=255)
                                                                            cTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTAc                                                                    <  2:909591/150‑1 (MQ=255)
                                                                            cTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTAc                                                                    >  2:702870/1‑150 (MQ=255)
                                                                              ggAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACg                                                                   >  1:894698/1‑149 (MQ=255)
                                                                               gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGca                                                                 >  2:961380/1‑149 (MQ=255)
                                                                               gAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGca                                                                 >  2:737096/1‑149 (MQ=255)
                                                                                aaaaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGcac                                                                >  2:460612/1‑148 (MQ=255)
                                                                                 aaaaGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCAcc                                                               <  2:449347/150‑1 (MQ=255)
                                                                                    agaagGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCgg                                                             >  1:450367/1‑149 (MQ=255)
                                                                                    agaagGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCgg                                                             >  1:924555/1‑149 (MQ=255)
                                                                                     gaagGCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCa                                                           <  2:13626/150‑1 (MQ=255)
                                                                                        ggCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCagaa                                                        <  1:537955/150‑1 (MQ=255)
                                                                                        ggCATTAACTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCagaa                                                        <  1:555921/150‑1 (MQ=255)
                                                                                         gCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCagaag                                                       <  2:240259/150‑1 (MQ=255)
                                                                                         gCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCagaag                                                       <  1:549671/150‑1 (MQ=255)
                                                                                            ttaaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATc                                                    >  1:256319/3‑150 (MQ=255)
                                                                                              aaCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATcc                                                   <  1:490149/149‑1 (MQ=255)
                                                                                               aCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCgg                                                 <  1:670078/150‑1 (MQ=255)
                                                                                               aCTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCgg                                                 <  1:407843/150‑1 (MQ=255)
                                                                                               aCTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCgg                                                 >  2:477572/1‑150 (MQ=255)
                                                                                                cTGTAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCggg                                                >  1:821478/1‑150 (MQ=255)
                                                                                                   tAACCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGt                                             <  1:315424/150‑1 (MQ=255)
                                                                                                     aCCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTg                                            <  2:551605/149‑1 (MQ=255)
                                                                                                     aCCTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGAGCGTgg                                           >  2:172701/1‑150 (MQ=255)
                                                                                                     aCCTGACCCTGCTGTTCCCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTg                                            <  1:167985/149‑1 (MQ=255)
                                                                                                      ccTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGt                                          >  2:270024/1‑150 (MQ=255)
                                                                                                       cTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGtt                                         <  2:587848/150‑1 (MQ=255)
                                                                                                       cTGACGCTGCTGTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGtt                                         >  1:146418/1‑150 (MQ=255)
                                                                                                          aCCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCt                                      >  2:87378/1‑150 (MQ=255)
                                                                                                           cgctgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTccg                                     <  1:371669/148‑3 (MQ=255)
                                                                                                            gctgctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGt                                    <  1:373842/149‑1 (MQ=255)
                                                                                                             ctgctgTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTa                                   >  1:674597/1‑150 (MQ=255)
                                                                                                                ctgTTCTCCTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGtcag                                 <  2:282456/149‑4 (MQ=255)
                                                                                                                 tgTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCt                                >  1:196215/1‑149 (MQ=255)
                                                                                                                   ttCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTg                               <  1:724285/148‑1 (MQ=255)
                                                                                                                        cTTCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTa                        <  1:459867/150‑1 (MQ=255)
                                                                                                                         ttCGCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTAc                       >  2:508650/1‑150 (MQ=255)
                                                                                                                            gCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAg                    <  1:758561/150‑1 (MQ=255)
                                                                                                                            gCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTAcga                     <  2:16227/149‑3 (MQ=255)
                                                                                                                            gCTCAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAg                    <  2:146418/150‑1 (MQ=255)
                                                                                                                               cAGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGtat                 >  1:261662/1‑149 (MQ=255)
                                                                                                                                aGGCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGtat                 <  2:899024/149‑2 (MQ=255)
                                                                                                                                 ggCTCGCGCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGtatta               >  1:593345/1‑150 (MQ=255)
                                                                                                                                       cgCTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaagca         <  1:39840/149‑4 (MQ=255)
                                                                                                                                         cTTGCGCCGAAGCGGGCGTCTTCCTGATTTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACaagcagc       <  2:632906/150‑6 (MQ=255)
                                                                                                                                            gtgCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACg    >  1:921337/1‑150 (MQ=255)
                                                                                                                                              gCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGt  >  2:774251/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGCTCAGGCTCGTGCTTGTGCGGAAGCGGGCGTGTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGT  >  CP009273/8597‑8888

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: