Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 646,796 C→A 6.9% D173Y (GAC→TAC)  citC ← [citrate [pro‑3S]‑lyase] ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273646,7960CA 6.9% 1221.9 / 48.9 317D173Y (GAC→TAC) citC[citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (154/141);  new base A (12/10);  total (166/151)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 320 total aligned reads displayed.

GGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTC  >  CP009273/646648‑646943
                                                                                                                                                    |                                                                                                                                                   
  gAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGc                                                                                                                                                   >  2:2555964/1‑149 (MQ=255)
    aCGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCAc                                                                                                                                                 <  1:4432641/149‑1 (MQ=255)
       tAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACtgt                                                                                                                                              >  1:4820500/1‑149 (MQ=255)
       tAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACtgt                                                                                                                                              >  1:3207312/1‑149 (MQ=255)
       tAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACtgt                                                                                                                                              <  2:4710055/149‑1 (MQ=255)
        aGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTAGCACtgtg                                                                                                                                             >  2:2048318/1‑149 (MQ=255)
         gcgcGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGc                                                                                                                                            >  2:3836451/1‑149 (MQ=255)
         gcgcGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTAGCACTGTGc                                                                                                                                            >  2:1326609/1‑149 (MQ=255)
             gggAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATAGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCa                                                                                                                                        <  1:3131309/149‑1 (MQ=255)
               gagaTGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGc                                                                                                                                      <  2:596611/149‑1 (MQ=255)
                   tgatgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGt                                                                                                                                  <  1:3327261/149‑1 (MQ=255)
                    gatgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGtt                                                                                                                                 >  1:4313690/1‑149 (MQ=255)
                      tgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTAGCACTGTGCCGCAGCCTGTTGa                                                                                                                               >  2:189/1‑149 (MQ=255)
                       gTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGaa                                                                                                                              >  1:3710182/1‑149 (MQ=255)
                          ttCGGAGCCACGATGCCCAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCa                                                                                                                           <  2:3467580/149‑1 (MQ=255)
                          ttCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCa                                                                                                                           >  2:2963171/1‑149 (MQ=255)
                            cGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGa                                                                                                                         <  2:2388327/149‑1 (MQ=255)
                                gCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAAc                                                                                                                     <  1:4777671/149‑1 (MQ=255)
                                gCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAAc                                                                                                                     <  2:2245296/149‑1 (MQ=255)
                                  cACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACgg                                                                                                                   >  1:3016377/1‑149 (MQ=255)
                                  cACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACgg                                                                                                                   >  2:2332545/1‑149 (MQ=255)
                                   aCGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGt                                                                                                                  <  1:1389322/149‑1 (MQ=255)
                                      aTGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGAc                                                                                                               <  2:285980/149‑1 (MQ=255)
                                        gCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCa                                                                                                             <  2:4010786/149‑1 (MQ=255)
                                        gCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCa                                                                                                             >  1:3518084/1‑149 (MQ=255)
                                         cacaGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCAt                                                                                                            <  2:395792/149‑1 (MQ=255)
                                              tCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTa                                                                                                       <  1:659255/149‑1 (MQ=255)
                                              tCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTa                                                                                                       <  2:139962/149‑1 (MQ=255)
                                                    gTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGa                                                                                                 >  2:4221288/1‑149 (MQ=255)
                                                    gTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGa                                                                                                 >  2:197130/1‑149 (MQ=255)
                                                    gTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGa                                                                                                 <  2:3926741/149‑1 (MQ=255)
                                                    gTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGa                                                                                                 >  1:1082694/1‑149 (MQ=255)
                                                     tGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGAt                                                                                                >  1:4914386/1‑149 (MQ=255)
                                                     tGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGAt                                                                                                >  2:3771039/1‑149 (MQ=255)
                                                     tGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTAGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGAt                                                                                                <  1:823436/149‑1 (MQ=255)
                                                         atatCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAAGCAGATAACGGTGACCATTCGTAAAGGGATTGGc                                                                                            >  1:125956/1‑149 (MQ=255)
                                                              ggcggTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTtca                                                                                       <  1:5201056/149‑1 (MQ=255)
                                                                  gTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATcac                                                                                   <  1:1126480/149‑1 (MQ=255)
                                                                  gTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATcac                                                                                   >  2:5155154/1‑149 (MQ=255)
                                                                    gCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACaa                                                                                 >  1:1550063/1‑149 (MQ=255)
                                                                      cTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATg                                                                               >  1:2034804/1‑149 (MQ=255)
                                                                       ttttAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGc                                                                              >  1:3619490/1‑149 (MQ=255)
                                                                       ttttAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGc                                                                              <  1:2481837/149‑1 (MQ=255)
                                                                           aaCACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGcc                                                                          >  1:184745/1‑149 (MQ=255)
                                                                            acacCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCa                                                                         >  2:2346928/1‑149 (MQ=255)
                                                                                 aaaTCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCtt                                                                    >  1:3929705/1‑149 (MQ=255)
                                                                                 aaaTCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCtt                                                                    <  1:64112/149‑1 (MQ=255)
                                                                                    tCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGtt                                                                 >  2:3317783/1‑149 (MQ=255)
                                                                                    tCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGtt                                                                 >  2:718030/1‑149 (MQ=255)
                                                                                    tCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGtt                                                                 >  1:3797185/1‑149 (MQ=255)
                                                                                     cccgCCGGACTTTATAGGTGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCTCACTGTGCAGCAGCATGTTGAATCAGATAACGGTGACCATTCGTAAATGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTc                                                                <  2:4514032/146‑1 (MQ=255)
                                                                                       aGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTccc                                                              <  2:3613819/149‑1 (MQ=255)
                                                                                       aGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTccc                                                              >  2:1665218/1‑149 (MQ=255)
                                                                                        gCCGGTCTTCATAGGGGTAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCt                                                             >  1:4688583/1‑149 (MQ=255)
                                                                                        gCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCt                                                             <  2:3279848/149‑1 (MQ=255)
                                                                                           ggTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGa                                                          >  2:4324289/1‑149 (MQ=255)
                                                                                           ggTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGa                                                          <  1:1683941/149‑1 (MQ=255)
                                                                                           ggTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTAGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGa                                                          >  1:4960545/1‑149 (MQ=255)
                                                                                             tCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATg                                                        >  2:2085066/1‑149 (MQ=255)
                                                                                              cTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGa                                                       >  2:538326/1‑149 (MQ=255)
                                                                                               ttCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGAc                                                      <  2:3679567/149‑1 (MQ=255)
                                                                                                   tAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGaaa                                                  <  2:2852977/149‑1 (MQ=255)
                                                                                                     ggggAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAtt                                                >  1:3868857/1‑149 (MQ=255)
                                                                                                      gggAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAttt                                               >  2:2839943/1‑149 (MQ=255)
                                                                                                        gAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATACCGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAttttt                                             <  2:4467607/149‑1 (MQ=255)
                                                                                                        gAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAttttt                                             >  1:2656722/1‑149 (MQ=255)
                                                                                                           gCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCa                                          <  2:4487990/149‑1 (MQ=255)
                                                                                                             gTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGc                                        <  2:105529/149‑1 (MQ=255)
                                                                                                              tGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAACTTTTTTCAGCg                                       >  2:4779190/1‑149 (MQ=255)
                                                                                                              tGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCg                                       >  2:1091752/1‑149 (MQ=255)
                                                                                                               gaagaaTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGa                                      <  1:3629632/149‑1 (MQ=255)
                                                                                                               gaagaaTCTTCTTTGACTAAAAACAGATGCAACCAGTAGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGa                                      <  1:189/149‑1 (MQ=255)
                                                                                                                 agaaTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGAtt                                    >  1:3673431/1‑149 (MQ=255)
                                                                                                                  gaaTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCCATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTc                                   >  1:539313/1‑149 (MQ=255)
                                                                                                                    aTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCgg                                 >  1:2874919/1‑149 (MQ=255)
                                                                                                                    aTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCgg                                 >  1:3419413/1‑149 (MQ=255)
                                                                                                                     tcttctTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGc                                >  2:559679/1‑149 (MQ=255)
                                                                                                                      cttcttTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCa                               <  2:3198417/149‑1 (MQ=255)
                                                                                                                       ttcttTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCAt                              >  2:1588199/1‑149 (MQ=255)
                                                                                                                         cttTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAg                            >  1:4060361/1‑149 (MQ=255)
                                                                                                                          tttGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGc                           >  2:179196/1‑149 (MQ=255)
                                                                                                                          tttGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGc                           >  1:788926/1‑149 (MQ=255)
                                                                                                                            tGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCATAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGt                         >  2:844664/1‑149 (MQ=255)
                                                                                                                               cTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTc                      <  1:2539135/149‑1 (MQ=255)
                                                                                                                                tAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCa                     >  2:240662/1‑149 (MQ=255)
                                                                                                                                  aaaaCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGt                   <  2:2441828/149‑1 (MQ=255)
                                                                                                                                  aaaaCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGt                   <  2:4946103/149‑1 (MQ=255)
                                                                                                                                   aaaCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTc                  >  1:2146481/1‑149 (MQ=255)
                                                                                                                                       aGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGt              >  1:1494554/1‑149 (MQ=255)
                                                                                                                                       aGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGt              <  1:837730/149‑1 (MQ=255)
                                                                                                                                       aGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGt              <  1:2364471/149‑1 (MQ=255)
                                                                                                                                       aGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGt              <  2:5137732/149‑1 (MQ=255)
                                                                                                                                             aaCCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCt        >  1:553513/1‑149 (MQ=255)
                                                                                                                                             aaCCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCt        <  1:1602975/149‑1 (MQ=255)
                                                                                                                                              aCCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTg       >  1:2063722/1‑149 (MQ=255)
                                                                                                                                              aCCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTg       >  2:1269517/1‑149 (MQ=255)
                                                                                                                                               ccAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGt      >  2:1134809/1‑149 (MQ=255)
                                                                                                                                                 aGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGttt    <  1:3020455/149‑1 (MQ=255)
                                                                                                                                                 aGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGttt    >  1:3263068/1‑149 (MQ=255)
                                                                                                                                                   tCGCACTGTGGCGCAGCATGTTGTATAAGATAACGGTGACCATTCGTAAAGGTATTGTAGTTCATCACAATGCAGCCAATCTTGTTACCTGGATGACGAATTTTTTTCAGCGATTCGGCATTGCGTTTCAGTCGCGTGGCGCTGTTTTc  <  1:3903941/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACACCAAATCCAGCCGGTCTTCATAGGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTC  >  CP009273/646648‑646943

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: