Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,618 T→A 90.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,6180TA90.3% 986.4 / 77.0 298V185E (GTG→GAG) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base A (120/149);  total (135/163)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 298 total aligned reads displayed.

CAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCG  >  CP009273/3080471‑3080764
                                                                                                                                                   |                                                                                                                                                  
cAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGc                                                                                                                                                  >  2:825557/1‑150 (MQ=255)
 aGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCt                                                                                                                                                 >  1:902523/1‑150 (MQ=255)
  gCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCtt                                                                                                                                                <  2:2218713/150‑1 (MQ=255)
   cacctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCttt                                                                                                                                               >  2:2485947/1‑150 (MQ=255)
    acctctcacctatGCACACCGTCTTGTACAGCGTCAGGCTTAATTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTc                                                                                                                                              <  1:863633/150‑1 (MQ=255)
    acctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTc                                                                                                                                              >  2:2600332/1‑150 (MQ=255)
      ctatcacctatGCACTCCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCa                                                                                                                                            <  1:3606768/150‑1 (MQ=255)
      ctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCa                                                                                                                                            <  1:2863615/150‑1 (MQ=255)
        atcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACt                                                                                                                                          <  2:321914/150‑1 (MQ=255)
          cacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCa                                                                                                                                        <  2:681130/150‑1 (MQ=255)
             ctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCa                                                                                                                                     <  1:2819682/150‑1 (MQ=255)
              tatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCAc                                                                                                                                    >  2:3400533/1‑150 (MQ=255)
              tatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCAc                                                                                                                                    <  2:1499188/150‑1 (MQ=255)
                tGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCActc                                                                                                                                  <  1:2994749/150‑1 (MQ=255)
                    cacCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGaa                                                                                                                              >  1:3153436/1‑150 (MQ=255)
                     acCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAag                                                                                                                             <  2:3658047/150‑1 (MQ=255)
                       cGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAagag                                                                                                                           >  2:4172731/1‑150 (MQ=255)
                       cGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAagag                                                                                                                           >  2:1676773/1‑150 (MQ=255)
                           tGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCg                                                                                                                       >  2:953787/1‑150 (MQ=255)
                            ggTACAGCGTCAGGCTGACGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGACGAGATCGa                                                                                                                      >  2:1922563/1‑150 (MQ=255)
                             gTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGAc                                                                                                                     <  2:1472272/150‑1 (MQ=255)
                              tACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGAcc                                                                                                                    <  2:565586/150‑1 (MQ=255)
                               aCAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCa                                                                                                                   <  1:1487357/150‑1 (MQ=255)
                                cAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAg                                                                                                                  >  2:2862178/1‑150 (MQ=255)
                                cAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAg                                                                                                                  >  2:2496940/1‑150 (MQ=255)
                                 aGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGa                                                                                                                 >  1:1675431/1‑150 (MQ=255)
                                   cGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGaaa                                                                                                               <  1:3399741/150‑1 (MQ=255)
                                    gTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAAt                                                                                                              <  2:729945/150‑1 (MQ=255)
                                      cAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCg                                                                                                            <  2:1786030/150‑1 (MQ=255)
                                        ggCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCt                                                                                                          >  2:3031755/1‑150 (MQ=255)
                                         tcTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCTACCAGAAATCGCTg                                                                                                         <  1:321736/149‑1 (MQ=255)
                                          cTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGc                                                                                                        >  1:3723544/1‑150 (MQ=255)
                                           tGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCa                                                                                                       <  2:3674066/150‑1 (MQ=255)
                                           tGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCa                                                                                                       <  1:3129019/150‑1 (MQ=255)
                                             aaGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaag                                                                                                     >  1:1316586/1‑150 (MQ=255)
                                             aaGTGCGTAAAAACGGCACTCTGCCGGGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTGGTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaag                                                                                                     <  2:3203255/150‑1 (MQ=255)
                                                tGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaagaag                                                                                                  <  1:472004/150‑1 (MQ=255)
                                                 gCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGc                                                                                                 <  1:1223392/150‑1 (MQ=255)
                                                    tAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGt                                                                                              >  1:2106962/1‑150 (MQ=255)
                                                    tAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGt                                                                                              <  1:3775930/150‑1 (MQ=255)
                                                      aaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTaa                                                                                            >  2:2777680/1‑150 (MQ=255)
                                                       aaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAt                                                                                           <  1:2052274/150‑1 (MQ=255)
                                                       aaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAt                                                                                           <  1:3694859/150‑1 (MQ=255)
                                                       aaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAt                                                                                           <  1:2626154/150‑1 (MQ=255)
                                                       aaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAt                                                                                           <  1:1900325/150‑1 (MQ=255)
                                                        aaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATg                                                                                          <  2:1596963/150‑1 (MQ=255)
                                                         aCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATgg                                                                                         <  1:4305979/150‑1 (MQ=255)
                                                              aCTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAagag                                                                                    >  2:3519743/1‑150 (MQ=255)
                                                               ctctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAgaga                                                                                   >  2:52639/1‑150 (MQ=255)
                                                                 ctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAACTCGCTGCAAGAAGCGGTAATGGAAGAGatc                                                                                 >  2:1880413/1‑150 (MQ=255)
                                                                  tGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGatca                                                                                >  2:1054887/1‑150 (MQ=255)
                                                                       tGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAg                                                                           >  2:2146857/1‑150 (MQ=255)
                                                                        ggCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGc                                                                          >  1:1780340/1‑150 (MQ=255)
                                                                         gCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGcc                                                                         <  1:4267907/150‑1 (MQ=255)
                                                                              gcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTc                                                                    <  2:1825381/150‑1 (MQ=255)
                                                                               cccGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCt                                                                   <  1:2496940/150‑1 (MQ=255)
                                                                                ccGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTg                                                                  <  2:3083341/150‑1 (MQ=255)
                                                                                  ggACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGcc                                                                >  2:2949731/1‑150 (MQ=255)
                                                                                   gACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGccc                                                               >  2:4231651/1‑150 (MQ=255)
                                                                                   gACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGccc                                                               <  1:3631125/150‑1 (MQ=255)
                                                                                     cgcgAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGACATCGCTGCAAGAAGCGGGAATGGAAGAGATCATCAAGCCAATTCTGCCCGc                                                             >  2:2220217/1‑150 (MQ=255)
                                                                                      tcgAAATGCCAGGTGTCTTGTCAGTATGACGTCGGCAAAATCGTTGGTATCGAAGCTGTCGAGCTTTCCACTCTTCACTCTGAAGAGATCGTCCAGAAATGGCTGGTAGAAGCGGTAATGGAAGAGATCTTATAGCCAATTCTTCCCGCt                                                            <  2:2204089/149‑1 (MQ=255)
                                                                                      gcgAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCt                                                            <  1:4348563/150‑1 (MQ=255)
                                                                                         aaaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGaa                                                         >  1:3864594/1‑150 (MQ=255)
                                                                                          aaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAAt                                                        >  1:1641332/1‑150 (MQ=255)
                                                                                          aaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAAt                                                        >  1:2703786/1‑150 (MQ=255)
                                                                                          aaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAAt                                                        >  1:2329937/1‑150 (MQ=255)
                                                                                            aGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATgg                                                      >  1:775247/1‑150 (MQ=255)
                                                                                             gCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGc                                                     >  1:3985245/1‑150 (MQ=255)
                                                                                              ccAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCt                                                    <  1:1750611/150‑1 (MQ=255)
                                                                                                 ggTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGAc                                                 >  2:1568659/1‑150 (MQ=255)
                                                                                                 ggTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGAc                                                 >  1:3275941/1‑150 (MQ=255)
                                                                                                    gACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTc                                              >  1:3709062/1‑150 (MQ=255)
                                                                                                     aCTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCt                                             <  1:3953382/150‑1 (MQ=255)
                                                                                                     aCTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCt                                             <  1:3500581/150‑1 (MQ=255)
                                                                                                     aCTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCt                                             <  2:902523/150‑1 (MQ=255)
                                                                                                         ttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGcca                                         <  2:4081680/150‑1 (MQ=255)
                                                                                                            aGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGccacca                                      <  2:614195/150‑1 (MQ=255)
                                                                                                              tATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCaaa                                    >  1:792425/1‑150 (MQ=255)
                                                                                                               aTGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAt                                   <  1:4186326/150‑1 (MQ=255)
                                                                                                                  acgacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttct                                <  1:3534851/150‑1 (MQ=255)
                                                                                                                   cgacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttctt                               <  1:2654624/150‑1 (MQ=255)
                                                                                                                     acgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtca                             >  1:965873/1‑150 (MQ=255)
                                                                                                                      cgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcat                            <  2:3153436/150‑1 (MQ=255)
                                                                                                                       ggCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcatc                           <  2:4256281/150‑1 (MQ=255)
                                                                                                                        gCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACCCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcatca                          <  1:1065853/150‑1 (MQ=255)
                                                                                                                         cAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCaa                         <  2:2378737/150‑1 (MQ=255)
                                                                                                                             aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                     <  2:1274699/150‑1 (MQ=255)
                                                                                                                             aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                     >  1:2229215/1‑150 (MQ=255)
                                                                                                                             aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGAGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                     >  1:2770621/1‑150 (MQ=255)
                                                                                                                               cGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAc                   >  2:4352766/1‑150 (MQ=255)
                                                                                                                               cGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAc                   <  1:1449174/150‑1 (MQ=255)
                                                                                                                               cGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAc                   <  1:697381/150‑1 (MQ=255)
                                                                                                                                gTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGACGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCACTTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAcc                  >  2:3057517/1‑150 (MQ=255)
                                                                                                                                gTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAcc                  >  2:3799553/1‑150 (MQ=255)
                                                                                                                                   ggTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGt               <  1:1844195/150‑1 (MQ=255)
                                                                                                                                   ggTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGt               <  1:4033531/150‑1 (MQ=255)
                                                                                                                                        cGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgttt          <  1:1548954/150‑1 (MQ=255)
                                                                                                                                         gATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc         <  1:557024/150‑1 (MQ=255)
                                                                                                                                         gATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAGTCTTCATCAACCCGACCGGgcgttgc         >  2:3358761/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCG  >  CP009273/3080471‑3080764

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: