Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,311,385 T→A 5.9% Q32L (CAG→CTG)  trpA ← tryptophan synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,311,3850TA 5.9% 1037.6 / 41.9 290Q32L (CAG→CTG) trpAtryptophan synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (117/155);  new base A (5/13);  total (122/168)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

ACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCG  >  CP009273/1311237‑1311533
                                                                                                                                                    |                                                                                                                                                    
 cATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTg                                                                                                                                                     >  2:2379259/1‑149 (MQ=255)
      tgCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCaa                                                                                                                                                <  2:4045763/149‑1 (MQ=255)
       gCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAAt                                                                                                                                               >  2:2796409/1‑149 (MQ=255)
         cGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGc                                                                                                                                             <  2:600196/149‑1 (MQ=255)
          ggAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGcc                                                                                                                                            >  1:682602/1‑149 (MQ=255)
          ggAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGcc                                                                                                                                            <  1:2318021/149‑1 (MQ=255)
             gTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCgg                                                                                                                                         >  1:333923/1‑149 (MQ=255)
             gTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATACGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCgg                                                                                                                                         >  1:359495/1‑149 (MQ=255)
              tCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGa                                                                                                                                        >  2:2943242/1‑149 (MQ=255)
                acacCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATc                                                                                                                                      >  1:427954/1‑149 (MQ=255)
                acacCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATc                                                                                                                                      >  2:101485/1‑149 (MQ=255)
                acacCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATc                                                                                                                                      >  1:428376/1‑149 (MQ=255)
                 cacCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCa                                                                                                                                     >  1:3465239/1‑149 (MQ=255)
                 cacCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATAGATAATTTTCAATGACTGCTCAATGACCGGATCa                                                                                                                                     >  1:3399272/1‑149 (MQ=255)
                  acCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCAc                                                                                                                                    <  2:1467207/149‑1 (MQ=255)
                   ccTGCCGCAAAGGCGCGCAGAGTGGCGTTTTTAATCGTCGGGCCATCCGCCCGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCAcc                                                                                                                                   >  1:1439476/1‑149 (MQ=255)
                     tGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGa                                                                                                                                 <  1:4905079/149‑1 (MQ=255)
                         gCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGGATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGt                                                                                                                             >  1:2735199/1‑149 (MQ=255)
                         gCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATACGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGt                                                                                                                             >  1:137266/1‑149 (MQ=255)
                            aaGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGAc                                                                                                                          >  2:2092412/1‑149 (MQ=255)
                            aaGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGAc                                                                                                                          <  1:1457238/149‑1 (MQ=255)
                                 gcgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAg                                                                                                                     <  1:4483558/149‑1 (MQ=255)
                                 gcgcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAg                                                                                                                     <  1:2420384/149‑1 (MQ=255)
                                   gcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGa                                                                                                                   <  1:1707751/149‑1 (MQ=255)
                                   gcAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGa                                                                                                                   <  2:1676139/149‑1 (MQ=255)
                                    cAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGaa                                                                                                                  <  2:1946007/149‑1 (MQ=255)
                                        gTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGaa                                                                                                              <  2:4666833/149‑1 (MQ=255)
                                         tGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAAt                                                                                                             <  2:3036990/149‑1 (MQ=255)
                                         tGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAAt                                                                                                             >  1:3686181/1‑149 (MQ=255)
                                         tGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAAt                                                                                                             <  2:4348294/149‑1 (MQ=255)
                                          ggCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATg                                                                                                            >  2:660771/1‑149 (MQ=255)
                                               tttGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCct                                                                                                       <  2:276836/149‑1 (MQ=255)
                                               tttGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGATTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGATCGAATGCGAct                                                                                                       <  2:276895/149‑1 (MQ=255)
                                                ttGAATCGTCGGGCCATCCGCCAGTGGGTCTGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCctt                                                                                                      <  2:4145268/149‑1 (MQ=255)
                                                  gAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCcttct                                                                                                    >  1:31024/1‑149 (MQ=255)
                                                   aaTCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCcttctt                                                                                                   <  1:3919192/149‑1 (MQ=255)
                                                    aTCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCttt                                                                                                  >  1:822813/1‑149 (MQ=255)
                                                    aTCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCttt                                                                                                  <  2:2640308/149‑1 (MQ=255)
                                                           ggCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTcc                                                                                           >  1:2815501/1‑149 (MQ=255)
                                                            gCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCt                                                                                          <  2:3917312/149‑1 (MQ=255)
                                                             ccATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCtt                                                                                         <  1:1900021/149‑1 (MQ=255)
                                                                   gCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTg                                                                                   <  1:3250383/149‑1 (MQ=255)
                                                                    ccAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTgg                                                                                  <  2:1424372/149‑1 (MQ=255)
                                                                    ccAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTgg                                                                                  <  2:3353060/149‑1 (MQ=255)
                                                                     cAGTGGGTCGGAGTAGGGGATTCCTAACTCCAGCGCGTCAGCACCGGCTTCATTTAGCGTATCGATAATTTTCTATGACTGCTCAATGCCCGGATCACAGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTggg                                                                                 <  1:4334560/149‑1 (MQ=255)
                                                                      aGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGc                                                                                <  1:4161349/149‑1 (MQ=255)
                                                                          ggTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAAc                                                                            <  2:681243/149‑1 (MQ=255)
                                                                           gTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACa                                                                           <  2:3826739/149‑1 (MQ=255)
                                                                            tCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACag                                                                          <  1:640711/149‑1 (MQ=255)
                                                                             cGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACaga                                                                         <  1:474877/149‑1 (MQ=255)
                                                                              ggAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACagag                                                                        >  2:958376/1‑149 (MQ=255)
                                                                              ggAGAAGGGGATAACTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACagag                                                                        >  2:4699857/1‑149 (MQ=255)
                                                                               gagaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAgaga                                                                       <  2:793977/149‑1 (MQ=255)
                                                                                agaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGAt                                                                      <  2:3303622/149‑1 (MQ=255)
                                                                                 gaAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGAtt                                                                     >  2:3372075/1‑149 (MQ=255)
                                                                                  aaGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTc                                                                    <  1:2490448/149‑1 (MQ=255)
                                                                                   aGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCg                                                                   <  1:2682318/149‑1 (MQ=255)
                                                                                     gggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTa                                                                 <  2:2463652/149‑1 (MQ=255)
                                                                                     gggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTa                                                                 <  2:4637987/149‑1 (MQ=255)
                                                                                     gggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTa                                                                 <  2:3927148/149‑1 (MQ=255)
                                                                                      ggATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAg                                                                <  2:3008178/149‑1 (MQ=255)
                                                                                          aCCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGtt                                                            <  2:1089398/149‑1 (MQ=255)
                                                                                            cTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTcc                                                          <  1:2544087/149‑1 (MQ=255)
                                                                                             tAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCa                                                         <  2:4244552/149‑1 (MQ=255)
                                                                                             tAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCa                                                         >  2:1327249/1‑149 (MQ=255)
                                                                                                cTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCa                                                      <  2:3690921/149‑1 (MQ=255)
                                                                                                 tCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAg                                                     >  1:4707923/1‑149 (MQ=255)
                                                                                                  ccAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGa                                                    >  1:197499/1‑149 (MQ=255)
                                                                                                   cAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACAGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGAt                                                   <  2:1573805/149‑1 (MQ=255)
                                                                                                        cgTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTccc                                              <  1:1047397/149‑1 (MQ=255)
                                                                                                          tCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCt                                            >  1:1100198/1‑149 (MQ=255)
                                                                                                          tCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCt                                            >  1:1766322/1‑149 (MQ=255)
                                                                                                             gCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGt                                         >  2:3642177/1‑149 (MQ=255)
                                                                                                                 cGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCtt                                     <  2:1760202/149‑1 (MQ=255)
                                                                                                                  ggCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCttt                                    <  2:3864697/149‑1 (MQ=255)
                                                                                                                   gCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTc                                   <  1:2091336/149‑1 (MQ=255)
                                                                                                                    cTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCa                                  >  2:2612158/1‑149 (MQ=255)
                                                                                                                      tCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTACTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCaaa                                >  1:1600807/1‑149 (MQ=255)
                                                                                                                       cAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCaaaa                               <  2:131868/149‑1 (MQ=255)
                                                                                                                       cAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCaaaa                               <  1:3166800/149‑1 (MQ=255)
                                                                                                                       cAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCaaaa                               <  2:358048/149‑1 (MQ=255)
                                                                                                                        aaTTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAat                              <  2:1227882/149‑1 (MQ=255)
                                                                                                                         aTTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAata                             >  2:508456/1‑149 (MQ=255)
                                                                                                                         aTTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAata                             <  1:3379310/149‑1 (MQ=255)
                                                                                                                          ttAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAatat                            <  2:2907880/149‑1 (MQ=255)
                                                                                                                           tAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATc                           >  2:3151982/1‑149 (MQ=255)
                                                                                                                            aGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCg                          <  2:899923/149‑1 (MQ=255)
                                                                                                                                tATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGaa                      >  1:4022809/1‑149 (MQ=255)
                                                                                                                                tATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGaa                      <  1:4685642/149‑1 (MQ=255)
                                                                                                                                 aTCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAAc                     >  1:5168684/1‑149 (MQ=255)
                                                                                                                                  tCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACg                    >  1:2728213/1‑149 (MQ=255)
                                                                                                                                   cGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACgg                   >  2:2182355/1‑149 (MQ=255)
                                                                                                                                     aTAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTg                 >  2:1063558/1‑149 (MQ=255)
                                                                                                                                     aTAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTg                 >  2:498176/1‑149 (MQ=255)
                                                                                                                                         ttttCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGa             <  2:4806760/149‑1 (MQ=255)
                                                                                                                                         ttttCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGa             <  1:2857179/149‑1 (MQ=255)
                                                                                                                                          tttCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGAt            <  1:4943331/149‑1 (MQ=255)
                                                                                                                                           ttCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATg           <  2:2695074/149‑1 (MQ=255)
                                                                                                                                              aaTGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAATGGTGAAGATGTCt        <  1:958376/149‑1 (MQ=255)
                                                                                                                                              aaTGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCt        <  1:3740289/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCG  >  CP009273/1311237‑1311533

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: