Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,601,793 C→A 6.4% S64Y (TCT→TAT)  tam → trans‑aconitate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,601,7930CA 6.4% 1411.9 / 44.0 361S64Y (TCT→TAT) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (151/187);  new base A (9/14);  total (160/201)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 378 total aligned reads displayed.

TGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTC  >  CP009273/1601644‑1601942
                                                                                                                                                     |                                                                                                                                                     
tGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtc                                                                                                                                                       >  2:1074699/1‑150 (MQ=255)
tGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtc                                                                                                                                                       <  2:133982/150‑1 (MQ=255)
 gAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtct                                                                                                                                                      <  2:4696056/150‑1 (MQ=255)
 gAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtct                                                                                                                                                      <  1:1645445/150‑1 (MQ=255)
  aaCGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtctc                                                                                                                                                     >  2:1612321/1‑150 (MQ=255)
  aaCGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGtctc                                                                                                                                                     <  2:4723521/150‑1 (MQ=255)
     gATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACACGTCTCCgg                                                                                                                                                  >  2:525317/1‑150 (MQ=255)
      aTCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGc                                                                                                                                                 <  1:420852/150‑1 (MQ=255)
      aTCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGc                                                                                                                                                 >  2:1554243/1‑150 (MQ=255)
        cgcgTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCga                                                                                                                                               >  1:4719101/1‑150 (MQ=255)
           gTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCgatga                                                                                                                                            >  2:1179328/1‑150 (MQ=255)
             ccGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGAtt                                                                                                                                          <  2:2527538/150‑1 (MQ=255)
             ccGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGAtt                                                                                                                                          <  1:227130/150‑1 (MQ=255)
                   gTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGaa                                                                                                                                    >  1:111495/1‑150 (MQ=255)
                    tGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAg                                                                                                                                   <  2:1239588/150‑1 (MQ=255)
                         ctgctTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAgcgcgc                                                                                                                              >  2:4335691/1‑150 (MQ=255)
                             ttGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTg                                                                                                                          <  1:2404825/150‑1 (MQ=255)
                             ttGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTg                                                                                                                          >  2:70657/1‑150 (MQ=255)
                              tGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGACGCGCGCAGTGc                                                                                                                         >  1:4854034/1‑150 (MQ=255)
                              tGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGc                                                                                                                         <  2:1657383/150‑1 (MQ=255)
                               gCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCt                                                                                                                        >  1:2502695/1‑150 (MQ=255)
                                  agagTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTg                                                                                                                     <  1:3506359/150‑1 (MQ=255)
                                  agagTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTg                                                                                                                     >  2:5039072/1‑150 (MQ=255)
                                      tGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAg                                                                                                                 >  1:2189992/1‑150 (MQ=255)
                                          gCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTg                                                                                                             <  2:3491910/150‑1 (MQ=255)
                                            tGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGcc                                                                                                           >  2:2490245/1‑150 (MQ=255)
                                             ggAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCa                                                                                                          <  1:2911713/150‑1 (MQ=255)
                                              gAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAg                                                                                                         <  2:2455572/150‑1 (MQ=255)
                                               aaaaTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGt                                                                                                        <  1:312709/150‑1 (MQ=255)
                                                aaaTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGtt                                                                                                       >  2:2125139/1‑150 (MQ=255)
                                                 aaTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTc                                                                                                      <  2:2775602/150‑1 (MQ=255)
                                                 aaTGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTc                                                                                                      >  2:3377602/1‑150 (MQ=255)
                                                    gTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTg                                                                                                   <  1:1707876/150‑1 (MQ=255)
                                                     tCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTgg                                                                                                  >  2:3635511/1‑150 (MQ=255)
                                                        aaTATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAg                                                                                               >  1:3669303/1‑150 (MQ=255)
                                                           atGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAg                                                                                            <  2:843443/150‑1 (MQ=255)
                                                             gTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGAc                                                                                          >  1:2335975/1‑150 (MQ=255)
                                                               tGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACAt                                                                                        <  1:161130/150‑1 (MQ=255)
                                                                gCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATc                                                                                       <  2:2296596/150‑1 (MQ=255)
                                                                 ccGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATcc                                                                                      >  1:4557707/1‑150 (MQ=255)
                                                                 ccGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATcc                                                                                      >  1:4000718/1‑150 (MQ=255)
                                                                   gcactgggCTGTGGCCCAGGTACCAGCACCGCCCTTCTACACCACCGCTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGc                                                                                    <  2:4365944/144‑1 (MQ=255)
                                                                   gATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGc                                                                                    >  1:2442065/1‑150 (MQ=255)
                                                                   gATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGc                                                                                    >  1:1436638/1‑150 (MQ=255)
                                                                    aTCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCa                                                                                   <  2:3153667/150‑1 (MQ=255)
                                                                       ttGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACt                                                                                <  1:874561/150‑1 (MQ=255)
                                                                        tGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTg                                                                               <  2:4149439/150‑1 (MQ=255)
                                                                        tGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTg                                                                               >  2:2097260/1‑150 (MQ=255)
                                                                         ggCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTgg                                                                              <  1:596831/150‑1 (MQ=255)
                                                                             gtgGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAAc                                                                          <  1:2009254/150‑1 (MQ=255)
                                                                              tgGCCCAGGTAACAGCACCGCCCTTCTACAACTACGTTGGCCTGCGGCCAGGATAACAGGCATGTACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTTCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAAcc                                                                         <  1:4240491/150‑1 (MQ=255)
                                                                               ggCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCg                                                                        <  2:4190896/150‑1 (MQ=255)
                                                                               ggCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCg                                                                        <  1:4572239/150‑1 (MQ=255)
                                                                                 cccAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGt                                                                      <  1:1858098/150‑1 (MQ=255)
                                                                                  ccAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCTGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTa                                                                     <  2:1461823/150‑1 (MQ=255)
                                                                                    aGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                   >  2:2660393/1‑150 (MQ=255)
                                                                                      gTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAgg                                                                 <  2:4558119/150‑1 (MQ=255)
                                                                                      gTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAgg                                                                 <  2:2198643/150‑1 (MQ=255)
                                                                                        aaCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCa                                                               <  1:282493/150‑1 (MQ=255)
                                                                                         aCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCAc                                                              <  1:1513923/150‑1 (MQ=255)
                                                                                          cAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACt                                                             >  2:4822366/1‑150 (MQ=255)
                                                                                             cACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGa                                                          <  1:1063536/150‑1 (MQ=255)
                                                                                               ccGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATc                                                        <  1:5108604/150‑1 (MQ=255)
                                                                                                cGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCTGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCt                                                       <  1:3836296/150‑1 (MQ=255)
                                                                                                cGCCCGTCTACAACAACCTTTGCCTGCGGCCAGGATAACAGGCCTCGACACCTCACCCGCCATGATTTCAGAAGCGCCCAGTTCGTTGCCAGACAGCCAGTTCGTGGGAGCAGACCTCCGCCACAGGCAACCGGTACAGGCACACCATca                                                       >  2:5096933/1‑149 (MQ=255)
                                                                                                  cccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGa                                                     >  1:983140/1‑150 (MQ=255)
                                                                                                  cccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGa                                                     <  1:1017049/150‑1 (MQ=255)
                                                                                                   ccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAt                                                    >  1:215226/1‑150 (MQ=255)
                                                                                                   ccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAt                                                    >  1:5020943/1‑150 (MQ=255)
                                                                                                     ttCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttt                                                  >  2:4455048/1‑150 (MQ=255)
                                                                                                       cTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttttt                                                <  2:3871087/150‑1 (MQ=255)
                                                                                                          caacaaCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCt                                             <  1:461511/150‑1 (MQ=255)
                                                                                                                 gTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCt                                      <  2:2010325/150‑1 (MQ=255)
                                                                                                                   tGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCa                                    >  1:2655827/1‑150 (MQ=255)
                                                                                                                      ccTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCTACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTg                                 <  1:2100341/150‑1 (MQ=255)
                                                                                                                        tGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCa                               <  2:4611233/150‑1 (MQ=255)
                                                                                                                        tGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCa                               <  2:1927331/150‑1 (MQ=255)
                                                                                                                           ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                            <  2:238818/150‑1 (MQ=255)
                                                                                                                            gCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATgg                           <  2:1897340/150‑1 (MQ=255)
                                                                                                                             ccAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGc                          >  1:3221188/1‑150 (MQ=255)
                                                                                                                              cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                         <  2:1115464/150‑1 (MQ=255)
                                                                                                                              cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                         <  1:1954297/150‑1 (MQ=255)
                                                                                                                              cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                         <  1:287897/150‑1 (MQ=255)
                                                                                                                                 gATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGcc                      <  1:3409667/150‑1 (MQ=255)
                                                                                                                                  aTAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGccc                     <  1:74338/150‑1 (MQ=255)
                                                                                                                                     aCAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAc                  >  2:3022403/1‑150 (MQ=255)
                                                                                                                                     aCAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAc                  >  1:4134509/1‑150 (MQ=255)
                                                                                                                                      cAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAcc                 <  2:289364/150‑1 (MQ=255)
                                                                                                                                      cAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGAcc                 >  2:696120/1‑150 (MQ=255)
                                                                                                                                       aGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCa                >  1:673309/1‑150 (MQ=255)
                                                                                                                                          cATCGACTCGTATCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTa             <  1:4829525/150‑1 (MQ=255)
                                                                                                                                            tCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACg           <  2:3762364/150‑1 (MQ=255)
                                                                                                                                              gACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa         <  1:1874694/150‑1 (MQ=255)
                                                                                                                                              gACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa         >  1:714446/1‑150 (MQ=255)
                                                                                                                                               aCTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAt        <  2:1804180/150‑1 (MQ=255)
                                                                                                                                                cTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAtt       <  2:3667231/150‑1 (MQ=255)
                                                                                                                                                 tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTg      >  2:3232792/1‑150 (MQ=255)
                                                                                                                                                 tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTg      <  1:1341860/150‑1 (MQ=255)
                                                                                                                                                   gTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGtt    >  2:4123719/1‑150 (MQ=255)
                                                                                                                                                    tctcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGttt   >  1:429462/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAACGATCGCGTCCGGCGGTGGAGCTGCTTGCCAGAGTGCCGCTGGAAAATGTCGAATATGTTGCCGATCTTGGCTGTGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTC  >  CP009273/1601644‑1601942

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: