Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,618 T→A 79.3% V185E (GTG→GAG)  metK → S‑adenosylmethionine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,080,6180TA79.3% 868.3 / 219.6 352V185E (GTG→GAG) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (28/45);  new base A (132/147);  total (160/192)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 352 total aligned reads displayed.

CCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGT  >  CP009273/3080470‑3080766
                                                                                                                                                    |                                                                                                                                                    
    cacctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCttt                                                                                                                                                 >  1:2840050/1‑150 (MQ=255)
    cacctatcacctagGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCttt                                                                                                                                                 <  1:632925/150‑1 (MQ=255)
        tatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCAc                                                                                                                                             >  1:2404796/1‑150 (MQ=255)
            acctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAg                                                                                                                                         >  2:4668251/1‑150 (MQ=255)
                  gCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCActct                                                                                                                                   >  2:3539551/1‑150 (MQ=255)
                  gCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCActct                                                                                                                                   <  1:550074/150‑1 (MQ=255)
                      acCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAag                                                                                                                               >  1:2169250/1‑150 (MQ=255)
                        cGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTGTCAGTATGACGACGGCCAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTCAagag                                                                                                                             >  2:5605079/1‑150 (MQ=255)
                           cTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATc                                                                                                                          <  2:1347448/150‑1 (MQ=255)
                              gTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGAc                                                                                                                       >  1:3899417/1‑150 (MQ=255)
                              gTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGAc                                                                                                                       >  2:471316/1‑150 (MQ=255)
                               tACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGAcc                                                                                                                      <  1:4979682/150‑1 (MQ=255)
                               tACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGAcc                                                                                                                      <  1:4156521/150‑1 (MQ=255)
                                   gCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGaa                                                                                                                  <  2:1843863/150‑1 (MQ=255)
                                     gTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAAt                                                                                                                <  1:4914682/150‑1 (MQ=255)
                                         ggCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCt                                                                                                            >  1:2626365/1‑150 (MQ=255)
                                         ggCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCt                                                                                                            >  2:348183/1‑150 (MQ=255)
                                          gCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTg                                                                                                           <  2:1231227/150‑1 (MQ=255)
                                          gCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTg                                                                                                           <  2:5088833/150‑1 (MQ=255)
                                            tGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCa                                                                                                         >  2:1893793/1‑150 (MQ=255)
                                            tGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCa                                                                                                         <  1:4280431/150‑1 (MQ=255)
                                            tGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCa                                                                                                         <  1:5389904/150‑1 (MQ=255)
                                              aaGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaag                                                                                                       <  1:1464162/150‑1 (MQ=255)
                                              aaGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaag                                                                                                       <  2:825821/150‑1 (MQ=255)
                                              aaGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaag                                                                                                       >  2:3632226/1‑150 (MQ=255)
                                              aaGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaag                                                                                                       >  2:4570306/1‑150 (MQ=255)
                                              aaGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaag                                                                                                       <  2:5477736/150‑1 (MQ=255)
                                               aGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaaga                                                                                                      <  2:4957535/150‑1 (MQ=255)
                                               aGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCcaga                                                                                                      >  2:5089657/1‑150 (MQ=255)
                                               aGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCcaga                                                                                                      >  2:5410543/1‑150 (MQ=255)
                                                gTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaagaa                                                                                                     <  1:5022012/150‑1 (MQ=255)
                                                 tGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaagaag                                                                                                    <  1:1388252/150‑1 (MQ=255)
                                                   cGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCg                                                                                                  <  1:1887812/150‑1 (MQ=255)
                                                    gTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCgg                                                                                                 >  1:1837379/1‑150 (MQ=255)
                                                      aaaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGACGAGATCGACCAGAAATCGCTGCAAGAAGCGGTa                                                                                               >  2:5168415/1‑150 (MQ=255)
                                                       aaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTaa                                                                                              <  1:3123464/150‑1 (MQ=255)
                                                       aaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTaa                                                                                              <  1:5626184/150‑1 (MQ=255)
                                                          aCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTAGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGACGCGGTAATgg                                                                                           >  2:5350375/1‑150 (MQ=255)
                                                           cGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGa                                                                                          >  2:2077199/1‑150 (MQ=255)
                                                            ggCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGaa                                                                                         <  2:2970075/150‑1 (MQ=255)
                                                             gCACTCTGCCGTGGCTGCGCCCGGACGCGAACAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAag                                                                                        <  2:3205803/150‑1 (MQ=255)
                                                             gCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAag                                                                                        <  2:1564216/150‑1 (MQ=255)
                                                               aCTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAagag                                                                                      <  2:3957761/150‑1 (MQ=255)
                                                                ctctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAgaga                                                                                     >  2:4930475/1‑150 (MQ=255)
                                                                 tctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGat                                                                                    >  1:4033882/1‑150 (MQ=255)
                                                                  ctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGatc                                                                                   >  1:1642530/1‑150 (MQ=255)
                                                                    gCCGTGGCTGCGCCCGGACGCGAACAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGatcat                                                                                 <  2:3349010/150‑1 (MQ=255)
                                                                      cGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGAtcatca                                                                               >  2:2621150/1‑150 (MQ=255)
                                                                       gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                              <  1:1381929/150‑1 (MQ=255)
                                                                       gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                              <  1:2335190/150‑1 (MQ=255)
                                                                         ggCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGc                                                                            >  2:1815669/1‑150 (MQ=255)
                                                                            tGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCaa                                                                         <  1:1739576/150‑1 (MQ=255)
                                                                               gcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTc                                                                      >  2:1111695/1‑150 (MQ=255)
                                                                               gcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTc                                                                      <  1:1772422/150‑1 (MQ=255)
                                                                               gcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGAAGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTc                                                                      <  2:2498068/150‑1 (MQ=255)
                                                                                  cGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTTAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGc                                                                   <  2:632056/150‑1 (MQ=255)
                                                                                    gACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGccc                                                                 >  1:966001/1‑150 (MQ=255)
                                                                                    gACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGccc                                                                 >  2:2373725/1‑150 (MQ=255)
                                                                                       gcgAAAAGCCAGGTTACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCt                                                              <  2:1027102/150‑1 (MQ=255)
                                                                                       gcgAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCt                                                              <  2:309921/150‑1 (MQ=255)
                                                                                           aaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAAt                                                          <  1:524482/150‑1 (MQ=255)
                                                                                              gCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGc                                                       <  2:842721/150‑1 (MQ=255)
                                                                                               ccAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCt                                                      >  2:404803/1‑150 (MQ=255)
                                                                                               ccAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCt                                                      >  2:5217538/1‑150 (MQ=255)
                                                                                                    tGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGACTGGCTGACtt                                                 >  2:1641905/1‑150 (MQ=255)
                                                                                                            cAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGccacc                                         >  2:4182224/1‑150 (MQ=255)
                                                                                                             aGTATGACGACGGCAAACTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCACTTCTGCCCGCTGAATGGCTGACTTCTGccacca                                        <  2:915858/150‑1 (MQ=255)
                                                                                                             aGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGccacca                                        <  1:1293267/150‑1 (MQ=255)
                                                                                                             aGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGccacca                                        >  2:680038/1‑150 (MQ=255)
                                                                                                              gTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCaa                                       <  1:4701882/150‑1 (MQ=255)
                                                                                                              gTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCaa                                       <  1:1066731/150‑1 (MQ=255)
                                                                                                              gTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCaa                                       <  1:4588770/150‑1 (MQ=255)
                                                                                                               tATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCaaa                                      >  2:2689349/1‑150 (MQ=255)
                                                                                                                     gacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttcttc                                <  1:1916121/150‑1 (MQ=255)
                                                                                                                        ggCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcatc                             <  1:5130194/150‑1 (MQ=255)
                                                                                                                        ggCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcatc                             <  1:2077467/150‑1 (MQ=255)
                                                                                                                           aaaaTCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAc                          >  1:3453100/1‑150 (MQ=255)
                                                                                                                           aaaaTCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTAAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAc                          >  2:1777458/1‑150 (MQ=255)
                                                                                                                            aaaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCCAcc                         >  1:1608790/1‑150 (MQ=255)
                                                                                                                             aaTCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAccc                        <  2:2031380/150‑1 (MQ=255)
                                                                                                                             aaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAccc                        <  1:909869/150‑1 (MQ=255)
                                                                                                                              aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                       <  1:2708893/150‑1 (MQ=255)
                                                                                                                                 gTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAcc                    <  2:1093222/150‑1 (MQ=255)
                                                                                                                                  ttGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCg                   <  1:737145/150‑1 (MQ=255)
                                                                                                                                  ttGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCg                   <  1:2849435/150‑1 (MQ=255)
                                                                                                                                  ttGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCg                   >  2:3809282/1‑150 (MQ=255)
                                                                                                                                   tGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCgg                  >  2:1482750/1‑150 (MQ=255)
                                                                                                                                   tGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCGTCATCAACCCGACCgg                  >  1:1963676/1‑150 (MQ=255)
                                                                                                                                     gTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtc                >  1:3062558/1‑150 (MQ=255)
                                                                                                                                         cGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgttt            <  2:4161470/150‑1 (MQ=255)
                                                                                                                                         cGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgttt            <  1:5473852/150‑1 (MQ=255)
                                                                                                                                          gATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc           <  1:3364113/150‑1 (MQ=255)
                                                                                                                                           aTGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg          <  2:1170905/150‑1 (MQ=255)
                                                                                                                                           aTGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg          <  1:4372463/150‑1 (MQ=255)
                                                                                                                                           aTGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg          <  2:4421876/150‑1 (MQ=255)
                                                                                                                                           aTGCTGTAGAGCTTTCCACTCAGCACTCTGAAGAGATCGACAAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg          <  2:2137389/150‑1 (MQ=255)
                                                                                                                                             gCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgtt        <  2:2850831/150‑1 (MQ=255)
                                                                                                                                               tGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGGCTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAt      >  1:375137/1‑150 (MQ=255)
                                                                                                                                                gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc     <  1:4179788/150‑1 (MQ=255)
                                                                                                                                                   gTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGt  >  1:1012192/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGT  >  CP009273/3080470‑3080766

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: