Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,573 A→C 62.3% A112A (GCA→GCC)  talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,5730AC62.3% 185.3 / 252.3 183A112A (GCA→GCC) talBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base A (31/38);  new base C (48/66);  total (79/104)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 185 total aligned reads displayed.

AACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGC  >  CP009273/8427‑8722
                                                                                                                                                  |                                                                                                                                                     
agcGATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCcaag                                                                                                                                                    <  2:191909/148‑2 (MQ=255)
 aCGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAAc                                                                                                                                                   <  2:239312/150‑1 (MQ=255)
  cGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACg                                                                                                                                                  >  2:293735/1‑150 (MQ=255)
   gATCGCGCGCAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCg                                                                                                                                                  >  1:294147/1‑149 (MQ=255)
    aTCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGcc                                                                                                                                                >  2:455152/1‑150 (MQ=255)
     tCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCt                                                                                                                                               <  1:627530/150‑1 (MQ=255)
       gcgcgcAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgcatca                                                                                                                                             >  2:10675/1‑146 (MQ=255)
       gcgcgcAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAgcgcatca                                                                                                                                             >  1:319030/1‑146 (MQ=255)
        cgcgcAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGAt                                                                                                                                            <  1:701776/150‑1 (MQ=255)
         gcgcAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATc                                                                                                                                           <  1:829808/150‑1 (MQ=255)
          cgcAGCAGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCatcatc                                                                                                                                           >  1:874030/1‑149 (MQ=255)
           gcagcaGATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCaa                                                                                                                                         <  2:511138/150‑1 (MQ=255)
             agcagATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAAc                                                                                                                                       <  1:662169/150‑1 (MQ=255)
             agcagATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAAc                                                                                                                                       <  1:80526/150‑1 (MQ=255)
               cagATCCTCGATGCTTCCGATAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGctg                                                                                                                                     <  2:987855/150‑2 (MQ=255)
                    cGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACaa                                                                                                                                >  1:647144/1‑150 (MQ=255)
                          cgcgACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGc                                                                                                                          <  1:324860/150‑1 (MQ=255)
                              tccgatAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCggg                                                                                                                       >  1:687404/2‑149 (MQ=255)
                               ccgatAAGCTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCgggca                                                                                                                     >  1:835260/1‑148 (MQ=255)
                                cgacAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTAt                                                                                                                    <  2:107370/150‑1 (MQ=255)
                                 gacAAACTGGCAGTAAATATTGGTCTGGACATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTAt                                                                                                                    >  2:790979/1‑149 (MQ=255)
                                  acAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCAAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTAt                                                                                                                    >  2:382544/1‑148 (MQ=255)
                                   taagcTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGGGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAg                                                                                                                 >  1:13713/5‑150 (MQ=255)
                                       cTGGCGGTGAACATCGGTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAAc                                                                                                             >  1:15351/1‑150 (MQ=255)
                                          gCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACg                                                                                                            <  2:152594/148‑1 (MQ=255)
                                              tAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTa                                                                                                      <  2:614303/150‑1 (MQ=255)
                                               gaacatcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCGCGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCAt                                                                                                     >  2:112324/8‑150 (MQ=255)
                                                aacatcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCAt                                                                                                     <  2:319030/143‑1 (MQ=255)
                                                 atatTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTc                                                                                                   <  1:519728/150‑1 (MQ=255)
                                                  catcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCt                                                                                                  >  1:101974/5‑150 (MQ=255)
                                                   atcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCt                                                                                                  >  2:984350/4‑149 (MQ=255)
                                                   atcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCt                                                                                                  >  2:572725/4‑149 (MQ=255)
                                                    tcggTCTCGAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGa                                                                                                <  2:709733/148‑1 (MQ=255)
                                                          tcgAAATTCTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAAc                                                                                          >  1:415496/3‑150 (MQ=255)
                                                            gAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACt                                                                                         <  2:870582/149‑1 (MQ=255)
                                                               attcTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGcc                                                                                     >  2:716414/4‑149 (MQ=255)
                                                                ttcTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGcct                                                                                    <  1:938442/148‑3 (MQ=255)
                                                                ttcTGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGcct                                                                                    <  1:429394/148‑3 (MQ=255)
                                                                   tGAAACTGATCCCGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTcca                                                                                 >  1:945407/1‑148 (MQ=255)
                                                                         tGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTgg                                                                            >  1:151833/1‑149 (MQ=255)
                                                                           gTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGc                                                                           <  1:234725/148‑1 (MQ=255)
                                                                            tcccGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAgg                                                                        <  2:645224/148‑1 (MQ=255)
                                                                               cGGGCCGTATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAgggca                                                                     >  2:130761/1‑148 (MQ=255)
                                                                                 ggCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATc                                                                   >  1:599510/1‑150 (MQ=255)
                                                                                    cGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGt                                                                >  2:313471/1‑150 (MQ=255)
                                                                                     gTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtg                                                               >  1:12935/1‑150 (MQ=255)
                                                                                      tATCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCgcgc                                                              >  1:181740/1‑150 (MQ=255)
                                                                                      tATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtgc                                                              <  1:796210/150‑1 (MQ=255)
                                                                                      tATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtgc                                                              >  1:553146/1‑150 (MQ=255)
                                                                                       aTCTCCACCGAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCgcgcc                                                             >  2:353990/1‑149 (MQ=255)
                                                                                          tcAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGtgctg                                                            >  1:560642/1‑148 (MQ=255)
                                                                                            caccgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGcgccgctg                                                         <  1:209908/145‑3 (MQ=255)
                                                                                              cTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAAc                                                      <  1:978900/150‑1 (MQ=255)
                                                                                               cgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGaa                                                       >  1:681745/2‑148 (MQ=255)
                                                                                               cgAAGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACa                                                     <  1:376037/149‑1 (MQ=255)
                                                                                                 aaGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACa                                                     <  1:767025/148‑1 (MQ=255)
                                                                                                  aGTGGACGCACGTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCt                                                  >  1:552727/1‑150 (MQ=255)
                                                                                                           gcgTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAaaga                                         <  1:707099/150‑1 (MQ=255)
                                                                                                            cgTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAaagaa                                        <  1:191245/150‑1 (MQ=255)
                                                                                                             gTCTCTCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAaagaag                                       <  2:549645/150‑1 (MQ=255)
                                                                                                                 ctCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAt                                   >  2:93557/2‑150 (MQ=255)
                                                                                                                  tCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGTCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGc                                     <  2:956759/147‑1 (MQ=255)
                                                                                                                  tCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCa                                    <  1:639691/148‑1 (MQ=255)
                                                                                                                  tCCTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCa                                    <  2:712743/148‑1 (MQ=255)
                                                                                                                    cTATGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAt                                   <  2:479924/147‑1 (MQ=255)
                                                                                                                    cTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCaa                                >  1:401619/1‑150 (MQ=255)
                                                                                                                    cTAAGACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCAtta                                 <  2:395826/149‑3 (MQ=255)
                                                                                                                        gACACCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTg                             >  2:726019/1‑149 (MQ=255)
                                                                                                                        gACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTg                             >  1:756412/1‑149 (MQ=255)
                                                                                                                        gACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTg                             >  1:737509/1‑149 (MQ=255)
                                                                                                                         acacCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTa                           <  1:460794/150‑1 (MQ=255)
                                                                                                                         acacCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTa                           <  2:553146/150‑1 (MQ=255)
                                                                                                                           acCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAAc                         <  2:923437/150‑1 (MQ=255)
                                                                                                                           acCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAAc                         <  2:43231/150‑1 (MQ=255)
                                                                                                                           acCGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAAc                         >  1:408321/1‑150 (MQ=255)
                                                                                                                             cGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                       <  1:902229/150‑1 (MQ=255)
                                                                                                                             cGAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCt                       <  2:573966/150‑1 (MQ=255)
                                                                                                                              gAAGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTg                      <  2:46284/150‑1 (MQ=255)
                                                                                                                                aGCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAc                    <  2:223235/150‑1 (MQ=255)
                                                                                                                                 gCGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAcg                   <  2:76606/150‑2 (MQ=255)
                                                                                                                                  cGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAcgc                  <  1:101803/150‑3 (MQ=255)
                                                                                                                                  cGTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGAcgc                  >  1:300901/1‑148 (MQ=255)
                                                                                                                                   gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                 >  1:405/1‑150 (MQ=255)
                                                                                                                                   gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                 <  1:349242/150‑1 (MQ=255)
                                                                                                                                   gTCTATCGCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGct                 <  1:112324/150‑1 (MQ=255)
                                                                                                                                    tctatcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctg                <  2:821156/144‑1 (MQ=255)
                                                                                                                                      aaTTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCctgct              >  1:157216/1‑150 (MQ=255)
                                                                                                                                        tcgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGctgctg             <  2:666305/147‑1 (MQ=255)
                                                                                                                                         cgCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt           >  2:231087/2‑150 (MQ=255)
                                                                                                                                          gCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGtt           <  1:565898/149‑1 (MQ=255)
                                                                                                                                          gCGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtc          <  1:572725/150‑1 (MQ=255)
                                                                                                                                           cGAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtct         <  1:39255/150‑1 (MQ=255)
                                                                                                                                           cGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtct         <  2:151833/150‑1 (MQ=255)
                                                                                                                                           cGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTtct         <  1:929977/150‑1 (MQ=255)
                                                                                                                                            gAAAGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTtctc        <  2:945407/150‑1 (MQ=255)
                                                                                                                                              aaGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCt      >  2:556175/1‑150 (MQ=255)
                                                                                                                                               aGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCtt     >  2:680939/1‑150 (MQ=255)
                                                                                                                                               aGCCAAGCGCATCATCAAGCTGTACAACGACGCGGGCATCAGCAACGATCGCATTCTGATCAAACTGGCCTCCACCTGGCAGGGCATTCGCGCCGCTGAACAGCTGGAAAAAGAAGGCATTAACTGTAACCTGACGCTGCTGTTCTCCtt     <  2:415496/150‑1 (MQ=255)
                                                                                                                                                 cAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCg   >  2:264589/1‑150 (MQ=255)
                                                                                                                                                  aaaaCGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCCCCTTCGc  >  2:167985/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AACGATCGCGCGCAGCAGATCGTGGACGCGACCGACAAACTGGCAGTAAATATTGGTCTGGAAATCCTGAAACTGGTTCCGGGCCGTATCTCAACTGAAGTTGATGCGCGTCTTTCCTATGACACCGAAGCGTCAATTGCGAAAGCAAAACGCCTGATCAAACTCTACAACGATGCTGGTATTAGCAACGATCGTATTCTGATCAAACTGGCTTCTACCTGGCAGGGTATCCGTGCTGCAGAACAGCTGGAAAAAGAAGGCATCAACTGTAACCTGACCCTGCTGTTCTCCTTCGC  >  CP009273/8427‑8722

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