Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,080,618 T→A 94.2% 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,6180TA94.2% 841.5 / 32.2 241V185E (GTG→GAG) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (9/5);  new base A (117/110);  total (126/115)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 241 total aligned reads displayed.

AGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATC  >  CP009273/3080472‑3080763
                                                                                                                                                  |                                                                                                                                                 
aGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCt                                                                                                                                                >  1:4279101/1‑150 (MQ=255)
  cacctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCttt                                                                                                                                              <  1:2023082/150‑1 (MQ=255)
   acctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGACAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGATATCGATGCTGTCGAGCTTTc                                                                                                                                             >  1:879957/1‑150 (MQ=255)
    cctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTcc                                                                                                                                            <  1:1935641/150‑1 (MQ=255)
     ctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCa                                                                                                                                           >  1:3032877/1‑150 (MQ=255)
     ctatcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCa                                                                                                                                           >  1:229513/1‑150 (MQ=255)
        tcacctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTc                                                                                                                                        >  1:4646208/1‑150 (MQ=255)
            ctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCa                                                                                                                                    >  1:2993426/1‑150 (MQ=255)
            ctatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCa                                                                                                                                    <  2:728354/150‑1 (MQ=255)
             tatGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCAc                                                                                                                                   <  1:2731950/150‑1 (MQ=255)
              atGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCAct                                                                                                                                  >  1:4166793/1‑150 (MQ=255)
              atGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCAct                                                                                                                                  >  1:3444296/1‑150 (MQ=255)
                gCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCActct                                                                                                                                <  2:2458161/150‑1 (MQ=255)
                 cacaCCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTg                                                                                                                               >  1:2323561/1‑150 (MQ=255)
                 cacaCCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTg                                                                                                                               <  1:4761505/150‑1 (MQ=255)
                  acacCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGa                                                                                                                              <  2:1642065/150‑1 (MQ=255)
                   cacCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGaa                                                                                                                             <  1:130594/150‑1 (MQ=255)
                         cTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATc                                                                                                                       >  1:3656715/1‑150 (MQ=255)
                         cTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATc                                                                                                                       <  2:2603622/150‑1 (MQ=255)
                         cTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATc                                                                                                                       >  1:1227545/1‑150 (MQ=255)
                         cTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATc                                                                                                                       >  1:1849759/1‑150 (MQ=255)
                          tGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGATGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGTTCg                                                                                                                      >  1:1975668/1‑150 (MQ=255)
                          tGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCg                                                                                                                      >  1:2193732/1‑150 (MQ=255)
                             tNCAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGAcc                                                                                                                   >  1:1368885/1‑150 (MQ=255)
                               cAGCGTCAGGCTGAGGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAg                                                                                                                 >  2:1046665/1‑150 (MQ=255)
                                  cGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGaaa                                                                                                              >  1:1681251/1‑150 (MQ=255)
                                   gTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAAt                                                                                                             <  2:439321/150‑1 (MQ=255)
                                     cAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCCGAAATCg                                                                                                           >  1:938077/1‑150 (MQ=255)
                                      aGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGc                                                                                                          >  1:603005/1‑150 (MQ=255)
                                        gCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTg                                                                                                        >  1:2397870/1‑150 (MQ=255)
                                           gAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaa                                                                                                     >  2:285335/1‑150 (MQ=255)
                                           gAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaa                                                                                                     <  1:3097706/150‑1 (MQ=255)
                                              gTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCaagaa                                                                                                  >  2:2683925/1‑150 (MQ=255)
                                                 cGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAcgcg                                                                                               >  2:4222200/1‑150 (MQ=255)
                                                 cGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCg                                                                                               >  1:3346054/1‑150 (MQ=255)
                                                   tAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGt                                                                                             >  1:4858035/1‑150 (MQ=255)
                                                    aaaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTa                                                                                            <  1:4043513/150‑1 (MQ=255)
                                                    aaaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTa                                                                                            >  1:4814654/1‑150 (MQ=255)
                                                     aaaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTaa                                                                                           <  2:3501222/150‑1 (MQ=255)
                                                      aaaCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAAt                                                                                          <  2:4330745/150‑1 (MQ=255)
                                                        aCGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATgg                                                                                        >  2:4160889/1‑150 (MQ=255)
                                                          ggCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGaa                                                                                      >  2:627102/1‑150 (MQ=255)
                                                            cACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTACTGGAaga                                                                                    >  1:3766879/1‑150 (MQ=255)
                                                              ctctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTCATGGAAgaga                                                                                  >  1:4798380/1‑150 (MQ=255)
                                                               tctGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGat                                                                                 >  2:3340877/1‑150 (MQ=255)
                                                                     gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                           <  1:650920/150‑1 (MQ=255)
                                                                     gTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCaa                                                                           <  2:4844230/150‑1 (MQ=255)
                                                                      tGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAg                                                                          >  1:606272/1‑150 (MQ=255)
                                                                       ggCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGc                                                                         <  2:1456107/150‑1 (MQ=255)
                                                                       ggCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGc                                                                         <  2:3586490/150‑1 (MQ=255)
                                                                          tGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCaa                                                                      >  2:4702570/1‑150 (MQ=255)
                                                                          tGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCaa                                                                      >  2:2875102/1‑150 (MQ=255)
                                                                          tGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCaa                                                                      >  1:3689788/1‑150 (MQ=255)
                                                                           gcgcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAAt                                                                     <  2:2633444/150‑1 (MQ=255)
                                                                            cgcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAAtt                                                                    >  2:3377392/1‑150 (MQ=255)
                                                                             gcCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTc                                                                   <  1:323375/150‑1 (MQ=255)
                                                                              cccGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCt                                                                  >  1:4372036/1‑150 (MQ=255)
                                                                                cGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGc                                                                >  2:2326266/1‑150 (MQ=255)
                                                                                cGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGc                                                                >  2:1790142/1‑150 (MQ=255)
                                                                                  gACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGccc                                                              >  1:4202923/1‑150 (MQ=255)
                                                                                  gACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGccc                                                              <  2:3207504/150‑1 (MQ=255)
                                                                                       gAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGa                                                         >  2:1306849/1‑150 (MQ=255)
                                                                                       gAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGa                                                         <  1:1031758/150‑1 (MQ=255)
                                                                                        aaaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGaa                                                        <  1:1525057/150‑1 (MQ=255)
                                                                                        aaaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGaa                                                        >  1:3541520/1‑150 (MQ=255)
                                                                                         aaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAAt                                                       <  2:5121822/150‑1 (MQ=255)
                                                                                         aaaGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAAt                                                       <  2:4629506/150‑1 (MQ=255)
                                                                                           aGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATgg                                                     <  2:1878225/150‑1 (MQ=255)
                                                                                             ccAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCt                                                   >  1:639544/1‑150 (MQ=255)
                                                                                             ccAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCt                                                   <  2:2753814/150‑1 (MQ=255)
                                                                                               aGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGa                                                 >  2:4545708/1‑150 (MQ=255)
                                                                                                ggTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGAc                                                >  1:422289/1‑150 (MQ=255)
                                                                                                  tGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACtt                                              >  1:2367437/1‑150 (MQ=255)
                                                                                                   gACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTc                                             <  1:5109766/150‑1 (MQ=255)
                                                                                                      ttttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGc                                          <  2:925654/150‑1 (MQ=255)
                                                                                                      ttttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGc                                          <  1:281931/150‑1 (MQ=255)
                                                                                                       tttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGcc                                         >  1:4963849/1‑150 (MQ=255)
                                                                                                       tttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGcc                                         >  1:2675863/1‑150 (MQ=255)
                                                                                                        ttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGcca                                        <  1:1397029/150‑1 (MQ=255)
                                                                                                        ttCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGcca                                        <  1:5140330/150‑1 (MQ=255)
                                                                                                              aTGACGACGGCAATATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAt                                  <  2:2830832/150‑1 (MQ=255)
                                                                                                              aTGACGACGGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAt                                  >  1:2758000/1‑150 (MQ=255)
                                                                                                                  cgacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttctt                              >  1:1458867/1‑150 (MQ=255)
                                                                                                                   gacgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAAttcttc                             <  2:3052728/150‑1 (MQ=255)
                                                                                                                    acgGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtca                            <  1:2943126/150‑1 (MQ=255)
                                                                                                                     cgGCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcat                           <  2:383545/150‑1 (MQ=255)
                                                                                                                       gCAAAATCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTtcatca                         >  1:25333/1‑150 (MQ=255)
                                                                                                                         aaaaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAc                       >  1:623543/1‑150 (MQ=255)
                                                                                                                         aaaaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAc                       >  1:3155559/1‑150 (MQ=255)
                                                                                                                          aaaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAcc                      >  2:3425692/1‑150 (MQ=255)
                                                                                                                           aaTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAAccc                     >  1:351018/1‑150 (MQ=255)
                                                                                                                            aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                    <  2:3032877/150‑1 (MQ=255)
                                                                                                                            aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                    <  1:3943247/150‑1 (MQ=255)
                                                                                                                            aTCGTTGGTATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCg                    <  2:3611777/150‑1 (MQ=255)
                                                                                                                                    tATCGATGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcg            <  2:1629811/150‑1 (MQ=255)
                                                                                                                                         aTGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg       <  2:4290497/150‑1 (MQ=255)
                                                                                                                                         aTGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtgtcg       >  1:4601365/1‑150 (MQ=255)
                                                                                                                                          tGCTGTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcgt      <  1:5192043/150‑1 (MQ=255)
                                                                                                                                              gTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc  <  2:892105/150‑1 (MQ=255)
                                                                                                                                              gTCGAGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATc  <  2:2059083/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCACCTATCACCTATGCACACCGTCTGGTACAGCGTCAGGCTGAAGTGCGTAAAAACGGCACTCTGCCGTGGCTGCGCCCGGACGCGAAAAGCCAGGTGACTTTTCAGTATGACGACGGCAAAATCGTTGGTATCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATC  >  CP009273/3080472‑3080763

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: