Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,336,681 A→T 24.1% Y155N (TAC→AAC)  fumB ← anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,336,6810AT24.1% 694.5 / 246.6 327Y155N (TAC→AAC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base A (132/116);  new base T (33/46);  total (165/162)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 336 total aligned reads displayed.

AGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGC  >  CP009273/4336532‑4336830
                                                                                                                                                     |                                                                                                                                                     
aGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTa                                                                                                                                                       >  2:3627969/1‑150 (MQ=255)
aGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTa                                                                                                                                                       >  2:3628043/1‑150 (MQ=255)
aGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTAAGTGCCGGTGTTGACCTCTTTGTa                                                                                                                                                       >  2:1737517/1‑150 (MQ=255)
    agcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGACTTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATg                                                                                                                                                   >  2:4493215/1‑150 (MQ=255)
    agcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATg                                                                                                                                                   <  1:5047085/150‑1 (MQ=255)
    agcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATg                                                                                                                                                   >  2:3961368/1‑150 (MQ=255)
     gcagGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGt                                                                                                                                                  >  2:1060521/1‑150 (MQ=255)
       agGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTcc                                                                                                                                                <  2:1776121/150‑1 (MQ=255)
        gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCa                                                                                                                                               >  2:1491516/1‑150 (MQ=255)
        gggCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCa                                                                                                                                               <  1:4480534/150‑1 (MQ=255)
                  ttCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCAtt                                                                                                                                     >  1:1013848/1‑150 (MQ=255)
                     cTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCtg                                                                                                                                  <  1:2086767/150‑1 (MQ=255)
                       ggTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCtgtg                                                                                                                                <  2:3523053/150‑1 (MQ=255)
                         tAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGaa                                                                                                                              <  2:3136956/150‑1 (MQ=255)
                              aTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgcg                                                                                                                         <  1:473612/150‑1 (MQ=255)
                              aTACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAgcg                                                                                                                         <  1:1474339/150‑1 (MQ=255)
                                 cGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAg                                                                                                                      <  1:970958/150‑1 (MQ=255)
                                    tttGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGtt                                                                                                                   <  1:1912993/150‑1 (MQ=255)
                                    tttGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGtt                                                                                                                   >  1:2202660/1‑150 (MQ=255)
                                    tttGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGtt                                                                                                                   >  1:1273747/1‑150 (MQ=255)
                                      tgttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTAt                                                                                                                 <  1:1091032/150‑1 (MQ=255)
                                       gttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATc                                                                                                                <  2:2350075/150‑1 (MQ=255)
                                        ttgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCt                                                                                                               >  1:3212362/1‑150 (MQ=255)
                                        ttgGCAGAGCCGCCGCCTTTCGAAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCt                                                                                                               >  2:3389413/1‑150 (MQ=255)
                                         tgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCtt                                                                                                              >  2:3273621/1‑150 (MQ=255)
                                         tgGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCtt                                                                                                              <  2:2906687/150‑1 (MQ=255)
                                         aggCAGAGCCGCCGACATTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCTCGTACAGGTAGATATGCGCAGGCAGGATAGTGACGGTGTTGACATCTTTGTACATGACCAGCTCCGCAATCTGTGAATAGCGCAGGTAATCtt                                                                                                              >  2:4737215/2‑150 (MQ=255)
                                          ggCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTc                                                                                                             <  1:738234/150‑1 (MQ=255)
                                           gCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCg                                                                                                            <  2:1226484/150‑1 (MQ=255)
                                            cAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGa                                                                                                           <  2:592301/150‑1 (MQ=255)
                                             agagCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGat                                                                                                          <  1:4644875/150‑1 (MQ=255)
                                             agagCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGat                                                                                                          >  2:401870/1‑150 (MQ=255)
                                                 ccgccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAg                                                                                                      <  1:2704384/150‑1 (MQ=255)
                                                 ccgccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAg                                                                                                      >  2:1378168/1‑150 (MQ=255)
                                                 ccgccgccTTTCGCAACGAAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTAAATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAg                                                                                                      >  2:1936999/1‑150 (MQ=255)
                                                   gccgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTCAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgt                                                                                                    <  2:1162100/150‑1 (MQ=255)
                                                     cgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgtgt                                                                                                  <  1:721537/150‑1 (MQ=255)
                                                     cgccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgtgt                                                                                                  <  1:1722957/150‑1 (MQ=255)
                                                      gccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGtt                                                                                                 >  2:2654213/1‑150 (MQ=255)
                                                      gccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGtt                                                                                                 >  1:521351/1‑150 (MQ=255)
                                                       ccTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTta                                                                                                <  2:4778763/150‑1 (MQ=255)
                                                        cTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtat                                                                                               >  1:3687878/1‑150 (MQ=255)
                                                        cTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtat                                                                                               >  2:2116188/1‑150 (MQ=255)
                                                         tttCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtata                                                                                              <  2:2530674/150‑1 (MQ=255)
                                                         tttCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTtata                                                                                              <  2:4949754/150‑1 (MQ=255)
                                                           tCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGa                                                                                            >  1:3801206/1‑150 (MQ=255)
                                                             gcaacgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACg                                                                                          >  1:2543872/1‑150 (MQ=255)
                                                              caacgcaaAGGAATTTGTACTCATCGCCATCAACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGATATCTTCGATATAGGTGTTATAGACGc                                                                                         >  1:1330570/1‑150 (MQ=255)
                                                                acgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCt                                                                                       >  2:71238/1‑150 (MQ=255)
                                                                acgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCt                                                                                       <  2:3785696/150‑1 (MQ=255)
                                                                acgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCt                                                                                       >  2:3208597/1‑150 (MQ=255)
                                                                 cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                      <  1:3139725/150‑1 (MQ=255)
                                                                 cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                      >  1:287053/1‑150 (MQ=255)
                                                                 cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                      <  2:2845041/150‑1 (MQ=255)
                                                                 cgcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCtt                                                                                      <  1:1566245/150‑1 (MQ=255)
                                                                  gcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCttt                                                                                     <  2:5071053/150‑1 (MQ=255)
                                                                  gcaaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCttt                                                                                     >  2:2495431/1‑150 (MQ=255)
                                                                   caaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTc                                                                                    <  2:4606550/150‑1 (MQ=255)
                                                                   caaAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTc                                                                                    <  1:823267/150‑1 (MQ=255)
                                                                       ggAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACa                                                                                <  1:4195339/150‑1 (MQ=255)
                                                                       ggAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACa                                                                                <  1:4578794/150‑1 (MQ=255)
                                                                       ggAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACa                                                                                <  1:4863424/150‑1 (MQ=255)
                                                                         aaTTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGc                                                                              <  2:1925733/150‑1 (MQ=255)
                                                                             tGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGttt                                                                          <  1:499148/150‑1 (MQ=255)
                                                                                  tCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtca                                                                     >  1:236250/1‑150 (MQ=255)
                                                                                  tCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtca                                                                     <  1:4840732/150‑1 (MQ=255)
                                                                                    aTCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTtcatc                                                                   >  2:4722142/1‑150 (MQ=255)
                                                                                      cGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAc                                                                 >  1:4045169/1‑150 (MQ=255)
                                                                                          aTCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgcc                                                             <  2:3301601/150‑1 (MQ=255)
                                                                                           tCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccg                                                            >  1:1237313/1‑150 (MQ=255)
                                                                                             tACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccgcc                                                          >  1:3828773/1‑150 (MQ=255)
                                                                                             tACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccgcc                                                          <  2:4296761/150‑1 (MQ=255)
                                                                                              nccGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCGTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCAccgccgccg                                                         >  2:805292/2‑150 (MQ=255)
                                                                                               ccGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCgg                                                        <  2:4114393/150‑1 (MQ=255)
                                                                                                  cgTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTcc                                                     >  1:304321/1‑150 (MQ=255)
                                                                                                    tACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCcac                                                   <  2:1631831/150‑1 (MQ=255)
                                                                                                    tACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCcac                                                   >  1:4789015/1‑150 (MQ=255)
                                                                                                     aCAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCcaca                                                  <  1:3397105/150‑1 (MQ=255)
                                                                                                       aGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcg                                                >  2:1859590/1‑150 (MQ=255)
                                                                                                       aGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACAcg                                                <  2:2961593/150‑1 (MQ=255)
                                                                                                           cGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCt                                            <  1:1021827/150‑1 (MQ=255)
                                                                                                           cGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCt                                            >  2:2144416/1‑150 (MQ=255)
                                                                                                            gATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGACTAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTACGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTg                                           >  1:2294949/1‑150 (MQ=255)
                                                                                                             aTTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTgg                                          <  2:5060031/150‑1 (MQ=255)
                                                                                                               ttGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGcc                                        >  1:2909455/1‑150 (MQ=255)
                                                                                                                tGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCt                                       >  1:1178746/1‑150 (MQ=255)
                                                                                                                tGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCt                                       <  1:904933/150‑1 (MQ=255)
                                                                                                                tGCGCAGGCAGGGTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCt                                       >  2:418840/1‑150 (MQ=255)
                                                                                                                  cgcAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt                                     <  1:2750141/150‑1 (MQ=255)
                                                                                                                      ggcaggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTa                                 <  2:535380/150‑1 (MQ=255)
                                                                                                                          ggTTAGTGCCGGTGTTGACCTCTTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcga                             >  2:2015413/1‑150 (MQ=255)
                                                                                                                          ggTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcga                             >  2:4264095/1‑150 (MQ=255)
                                                                                                                            ttAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACcgacg                           <  2:4079865/150‑1 (MQ=255)
                                                                                                                                tGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACgatga                       <  1:560942/150‑1 (MQ=255)
                                                                                                                                         tGACCTCTTTGTTCAGGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGc              <  1:3513568/150‑1 (MQ=255)
                                                                                                                                         tGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGc              >  1:4158813/1‑150 (MQ=255)
                                                                                                                                             ctctTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGt          >  2:2883472/1‑150 (MQ=255)
                                                                                                                                               ctTTGTTCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTAt        >  1:1138091/1‑150 (MQ=255)
                                                                                                                                                tttGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATc       >  1:3995450/1‑150 (MQ=255)
                                                                                                                                                tttGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATc       >  1:1433435/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGATTTGCGCAGGCAGGTTAGTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGC  >  CP009273/4336532‑4336830

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: