Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,867,062 C→A 40.5% R258L (CGC→CTC)  dgoK ← 2‑oxo‑3‑deoxygalactonate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,867,0620CA40.5% 350.4 / inf 291R258L (CGC→CTC) dgoK2‑oxo‑3‑deoxygalactonate kinase
Reads supporting (aligned to +/- strand):  ref base C (84/89);  new base A (47/71);  total (131/160)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 295 total aligned reads displayed.

CGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGG  >  CP009273/3866915‑3867209
                                                                                                                                                   |                                                                                                                                                   
cgtagGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAgagc                                                                                                                                                   >  2:4580290/1‑150 (MQ=255)
   cgGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGcgcgcg                                                                                                                                                >  1:1417654/1‑150 (MQ=255)
     ggCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGt                                                                                                                                              <  2:1020362/150‑1 (MQ=255)
       cGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCa                                                                                                                                            <  1:3108885/150‑1 (MQ=255)
        gTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAg                                                                                                                                           <  1:2852443/150‑1 (MQ=255)
         tAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGc                                                                                                                                          >  1:1233845/1‑150 (MQ=255)
         tAATACAGCGCATAATTGCGATCTGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCCGCCTGAAATTCAGTGTCGCCCGCCACAGCAGTCACGTAGCAACCCATCGCCTGAAAGGCTTTCTGGTAGCGCGCGGTCatc                                                                                                                                          <  1:1966846/150‑3 (MQ=255)
          aaTACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCg                                                                                                                                         >  1:3418145/1‑150 (MQ=255)
          aaTACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCg                                                                                                                                         <  2:3749594/150‑1 (MQ=255)
          aaTACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCg                                                                                                                                         >  2:3286976/1‑150 (MQ=255)
              ccGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGt                                                                                                                                     >  2:3469496/1‑150 (MQ=255)
               cgcgCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGtt                                                                                                                                    <  2:3577761/150‑1 (MQ=255)
                  gcAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCg                                                                                                                                 <  1:4586424/150‑1 (MQ=255)
                    aaaaTGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGc                                                                                                                               >  1:2927198/1‑150 (MQ=255)
                     aaaTGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTTTACCTGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCg                                                                                                                              <  1:1143362/150‑1 (MQ=255)
                     aaaTGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCg                                                                                                                              >  1:2531693/1‑150 (MQ=255)
                     aaaTGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCg                                                                                                                              <  2:4199234/150‑1 (MQ=255)
                      aaTGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGa                                                                                                                             >  1:1258259/1‑150 (MQ=255)
                        tGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACa                                                                                                                           <  2:1389659/150‑1 (MQ=255)
                        tGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACa                                                                                                                           <  2:1389583/150‑1 (MQ=255)
                          gCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAg                                                                                                                         >  2:2297903/1‑150 (MQ=255)
                           cGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAgg                                                                                                                        >  2:2928228/1‑150 (MQ=255)
                           cGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAgg                                                                                                                        >  2:3674918/1‑150 (MQ=255)
                           cGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAgg                                                                                                                        <  1:844582/150‑1 (MQ=255)
                           cGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGCTCCGGCGACAAgg                                                                                                                        >  2:3674913/1‑150 (MQ=255)
                               cAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCAGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGa                                                                                                                    >  2:1668064/1‑150 (MQ=255)
                                     gAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCgt                                                                                                              >  2:3425784/1‑150 (MQ=255)
                                     gAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCgt                                                                                                              <  1:2387043/150‑1 (MQ=255)
                                        tttagtttGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGtt                                                                                                           <  2:3022462/150‑1 (MQ=255)
                                            gtttGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATg                                                                                                       <  1:2078501/150‑1 (MQ=255)
                                            gtttGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATg                                                                                                       <  1:3842388/150‑1 (MQ=255)
                                              ttGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACCAGGGTGATGGCGTGTTGATggg                                                                                                     >  2:3372129/1‑150 (MQ=255)
                                               tGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGc                                                                                                    >  1:1692807/1‑150 (MQ=255)
                                                 ccACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCca                                                                                                  <  2:614186/150‑1 (MQ=255)
                                                 ccACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCca                                                                                                  >  1:2416753/1‑150 (MQ=255)
                                                   aCTGCATGAGCGATGCTCCTTATCCCAGCCTGAAATGCCGTGTCGCACGACACCGCCGTCACGTCGCAACCCCTCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCcaca                                                                                                <  2:1582203/150‑1 (MQ=255)
                                                   aCTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCcaca                                                                                                >  2:3830801/1‑150 (MQ=255)
                                                    cTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACAt                                                                                               >  1:322315/1‑150 (MQ=255)
                                                     tGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATa                                                                                              <  1:4340099/150‑1 (MQ=255)
                                                      gCATGAGCGATGCTCCTTATACCATCCTGAAATGCAGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGATGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAg                                                                                             >  2:4402209/1‑150 (MQ=255)
                                                        aTGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTc                                                                                           <  2:485125/150‑1 (MQ=255)
                                                        aTGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTc                                                                                           >  2:173895/1‑150 (MQ=255)
                                                          gAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTcgc                                                                                         >  2:4494815/1‑150 (MQ=255)
                                                           aGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTcgcg                                                                                        <  1:3389718/150‑1 (MQ=255)
                                                           aGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTcgcg                                                                                        >  2:74912/1‑150 (MQ=255)
                                                               aTGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATa                                                                                    <  2:2007167/150‑1 (MQ=255)
                                                               aTGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATa                                                                                    <  2:1177337/150‑1 (MQ=255)
                                                               aTGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATa                                                                                    <  1:1891899/150‑1 (MQ=255)
                                                                 gCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACt                                                                                  <  1:615180/150‑1 (MQ=255)
                                                                 gCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACt                                                                                  <  1:4423409/150‑1 (MQ=255)
                                                                 gCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACt                                                                                  >  2:2136656/1‑150 (MQ=255)
                                                                  cTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTg                                                                                 >  2:621687/1‑150 (MQ=255)
                                                                  cTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTg                                                                                 >  2:3128120/1‑150 (MQ=255)
                                                                     cTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCAGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCg                                                                              <  2:2598158/150‑1 (MQ=255)
                                                                      ttATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGa                                                                             >  1:2120397/1‑150 (MQ=255)
                                                                      ttATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGa                                                                             >  1:1599384/1‑150 (MQ=255)
                                                                       tataCCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGAc                                                                            >  1:2417165/1‑150 (MQ=255)
                                                                          aCCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACctc                                                                         >  2:89425/1‑150 (MQ=255)
                                                                           ccAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACctct                                                                        <  2:1211494/150‑1 (MQ=255)
                                                                           ccAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACctct                                                                        <  1:239246/150‑1 (MQ=255)
                                                                            cAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTg                                                                       >  1:2667341/1‑150 (MQ=255)
                                                                             aGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgc                                                                      <  1:3271203/150‑1 (MQ=255)
                                                                             aGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgc                                                                      <  1:420233/150‑1 (MQ=255)
                                                                              gCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgcg                                                                     >  1:1011273/1‑150 (MQ=255)
                                                                              gCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgcg                                                                     <  1:3206392/150‑1 (MQ=255)
                                                                               ccTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTgcgc                                                                    >  1:4188892/1‑150 (MQ=255)
                                                                                  gAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCaa                                                                 <  1:188339/150‑1 (MQ=255)
                                                                                  gAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCaa                                                                 >  1:1566055/1‑150 (MQ=255)
                                                                                  gAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCaa                                                                 >  1:3953974/1‑150 (MQ=255)
                                                                                    aaTGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATc                                                               >  1:3516308/1‑150 (MQ=255)
                                                                                    aaTGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATc                                                               <  2:1110185/150‑1 (MQ=255)
                                                                                    aaTGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATc                                                               <  1:2190500/150‑1 (MQ=255)
                                                                                        ccGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATcaaca                                                           >  1:1783849/1‑150 (MQ=255)
                                                                                         cGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATcaacaa                                                          <  2:4506497/150‑1 (MQ=255)
                                                                                              cgcccgccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCag                                                     >  2:848831/1‑150 (MQ=255)
                                                                                                cccgccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCagag                                                   <  2:94794/150‑1 (MQ=255)
                                                                                                cccgccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCagag                                                   >  2:3942496/1‑150 (MQ=255)
                                                                                                  cgccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCagagag                                                 >  2:362528/1‑150 (MQ=255)
                                                                                                   gccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAgagaga                                                <  1:771193/150‑1 (MQ=255)
                                                                                                    ccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGaa                                               >  1:3693042/1‑150 (MQ=255)
                                                                                                      aCCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAAt                                             <  2:1312591/150‑1 (MQ=255)
                                                                                                         gccgTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCg                                          <  1:274712/150‑1 (MQ=255)
                                                                                                           cgTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCt                                        >  2:4057868/1‑150 (MQ=255)
                                                                                                            gTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTg                                       <  1:964306/150‑1 (MQ=255)
                                                                                                                  tCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGt                                 <  2:2773452/150‑1 (MQ=255)
                                                                                                                   cGCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGtt                                >  2:3911880/1‑150 (MQ=255)
                                                                                                                    gCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTc                               >  1:415743/1‑150 (MQ=255)
                                                                                                                    gCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTc                               <  1:205984/150‑1 (MQ=255)
                                                                                                                    gCAACCCATCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTc                               <  2:537723/150‑1 (MQ=255)
                                                                                                                      aaCCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTcgc                             >  2:1911098/1‑150 (MQ=255)
                                                                                                                           aTCGCCTGAAAGGCTTGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGa                        <  1:2876831/150‑1 (MQ=255)
                                                                                                                            tCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGaa                       <  1:92606/150‑1 (MQ=255)
                                                                                                                            tCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGaa                       <  1:3663494/150‑1 (MQ=255)
                                                                                                                               ccTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAgtg                    >  1:1468383/1‑150 (MQ=255)
                                                                                                                               ccTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAgtg                    <  2:1937185/150‑1 (MQ=255)
                                                                                                                                      ggCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAg             <  2:2525799/150‑1 (MQ=255)
                                                                                                                                         ttGCTGGTAGAGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGcac          <  2:2579140/150‑1 (MQ=255)
                                                                                                                                               gTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGa    >  1:3808839/1‑150 (MQ=255)
                                                                                                                                                 aGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCTGGGAAGTGTTCCCAGCACATGCGAgg  <  2:1481107/150‑1 (MQ=255)
                                                                                                                                                 aGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAgg  <  1:1007081/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGTCGGGCGTAATACCGCGCAAAATGGCGATCAGCGGGAGTTTAGTTTGCCACTGCATGAGCGATGCTCCTTATACCAGCCTGAAATGCCGTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGG  >  CP009273/3866915‑3867209

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: ‑