Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,532,274 A→T 5.2% A223A (GCT→GCA cysM ← cysteine synthase B (O‑acetylserine sulfhydrolase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,532,2740AT 5.2% 949.1 / 31.6 295A223A (GCT→GCAcysMcysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base A (137/141);  new base T (10/6);  total (148/147)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 309 total aligned reads displayed.

GCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCC  >  CP009273/2532125‑2532423
                                                                                                                                                     |                                                                                                                                                     
gcgTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGaa                                                                                                                                                       <  2:2743444/150‑1 (MQ=255)
 cgTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAg                                                                                                                                                      <  2:1081630/150‑1 (MQ=255)
  gTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGc                                                                                                                                                     <  2:3054751/150‑1 (MQ=255)
    cAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGt                                                                                                                                                   >  2:3141287/1‑150 (MQ=255)
       ggTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGa                                                                                                                                                >  2:1011069/1‑150 (MQ=255)
        gTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGaa                                                                                                                                               >  2:2951471/1‑150 (MQ=255)
         ttAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGaaa                                                                                                                                              <  2:1182959/150‑1 (MQ=255)
            gCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATc                                                                                                                                           <  1:3647135/150‑1 (MQ=255)
              ttttGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCGGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATccc                                                                                                                                         >  1:4103431/1‑150 (MQ=255)
               tttGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATcccc                                                                                                                                        >  1:61841/1‑150 (MQ=255)
                ttGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCg                                                                                                                                       >  1:3800341/1‑150 (MQ=255)
                ttGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCg                                                                                                                                       >  1:3290973/1‑150 (MQ=255)
                 tGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCgg                                                                                                                                      >  1:419089/1‑150 (MQ=255)
                      cccGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCCGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACAAGAGAAGCGTTGAAAATCCCCGGCAGat                                                                                                                                 <  1:1486662/150‑1 (MQ=255)
                      cccGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGat                                                                                                                                 >  1:3278425/1‑150 (MQ=255)
                         gCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATAtt                                                                                                                              <  1:3097029/150‑1 (MQ=255)
                         gCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATAtt                                                                                                                              <  1:252007/150‑1 (MQ=255)
                           aGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTcc                                                                                                                            >  1:3511512/1‑150 (MQ=255)
                            gTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATAGTCCg                                                                                                                           <  1:3145989/150‑1 (MQ=255)
                              gCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTa                                                                                                                         >  1:3885816/1‑150 (MQ=255)
                               cTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGGGGAAAATCCCCGGCAGATATTCCGTAg                                                                                                                        >  2:1701471/1‑150 (MQ=255)
                                   ggCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCa                                                                                                                    >  2:2533537/1‑150 (MQ=255)
                                   ggCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCa                                                                                                                    >  1:722091/1‑150 (MQ=255)
                                   ggCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCa                                                                                                                    >  1:4733170/1‑150 (MQ=255)
                                     cAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAgc                                                                                                                  >  2:3085134/1‑150 (MQ=255)
                                      aaCCGCGCCGCCGGAGCTGACGCCACAGACTATTCCGTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAgcg                                                                                                                 <  2:2318319/150‑1 (MQ=255)
                                      aaCCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAgcg                                                                                                                 <  1:1401281/150‑1 (MQ=255)
                                       aCCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAgcgg                                                                                                                >  1:5027747/1‑150 (MQ=255)
                                         cgcgCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAgcggcg                                                                                                              <  1:4722360/150‑1 (MQ=255)
                                         cgcgCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAgcggcg                                                                                                              >  1:3272131/1‑150 (MQ=255)
                                            gccgccGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGATGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAAt                                                                                                           >  1:2530252/1‑150 (MQ=255)
                                             ccgccgGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATg                                                                                                          <  1:1892077/150‑1 (MQ=255)
                                               gccgGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGcc                                                                                                        >  1:921172/1‑150 (MQ=255)
                                                ccgGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCg                                                                                                       <  1:2570879/150‑1 (MQ=255)
                                                 cgGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCgg                                                                                                      >  2:3518666/1‑150 (MQ=255)
                                                  ggAGCTGACGCCACAGAATATTCCGTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCggg                                                                                                     <  1:5142525/150‑1 (MQ=255)
                                                    aGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGaa                                                                                                   >  2:4395331/1‑150 (MQ=255)
                                                      cTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtg                                                                                                 >  2:4895576/1‑150 (MQ=255)
                                                      cTGACGCCACAGAATATTCCGTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtg                                                                                                 <  2:3243892/150‑1 (MQ=255)
                                                       tGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtgc                                                                                                >  2:1478997/1‑150 (MQ=255)
                                                       tGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtgc                                                                                                <  1:4151765/150‑1 (MQ=255)
                                                        gACGCCCCAGCATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtgct                                                                                               <  2:5051323/150‑1 (MQ=255)
                                                           gCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtgctgct                                                                                            <  2:939862/150‑1 (MQ=255)
                                                            ccACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtgctgctg                                                                                           <  2:4719590/150‑1 (MQ=255)
                                                             cacaGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAAtgctgctgc                                                                                          >  1:892237/1‑150 (MQ=255)
                                                              acaGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGcc                                                                                         <  1:4507861/150‑1 (MQ=255)
                                                                aGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCct                                                                                       <  1:1619646/150‑1 (MQ=255)
                                                                   atatTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCtt                                                                                    <  1:1687367/150‑1 (MQ=255)
                                                                     atTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTcc                                                                                  <  2:1835759/150‑1 (MQ=255)
                                                                         cttccCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGtt                                                                              <  1:3144329/150‑1 (MQ=255)
                                                                           tccCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGc                                                                            <  2:4823799/150‑1 (MQ=255)
                                                                            cccGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCa                                                                           >  1:1816857/1‑150 (MQ=255)
                                                                            cccGCACCGCCAGTTCGCGCATGGTGATTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGATGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGTGAATGCTGCTGCCCTCTTCCGGTTGCa                                                                           <  1:2879962/150‑1 (MQ=255)
                                                                               gCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGc                                                                        <  2:264762/150‑1 (MQ=255)
                                                                                   cGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGAc                                                                    >  2:4288543/1‑150 (MQ=255)
                                                                                     ccAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACaa                                                                  >  1:2817506/1‑150 (MQ=255)
                                                                                      cAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAAt                                                                 <  1:4640472/150‑1 (MQ=255)
                                                                                       aGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATg                                                                >  2:383099/1‑150 (MQ=255)
                                                                                              gcATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCg                                                         <  1:3564602/150‑1 (MQ=255)
                                                                                              gcATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCg                                                         >  2:52924/1‑150 (MQ=255)
                                                                                              gcATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATAGTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCg                                                         <  1:2322579/150‑1 (MQ=255)
                                                                                                  ggTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGttt                                                     <  2:866305/150‑1 (MQ=255)
                                                                                                   gtgtTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTg                                                    <  1:3741842/150‑1 (MQ=255)
                                                                                                   gtgtTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTg                                                    <  2:2901105/150‑1 (MQ=255)
                                                                                                   gtgtTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTg                                                    <  1:1169791/150‑1 (MQ=255)
                                                                                                     gtTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGa                                                  <  1:4395331/150‑1 (MQ=255)
                                                                                                      ttttCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATAGTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGAt                                                 <  1:1536873/150‑1 (MQ=255)
                                                                                                      ttttCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTAGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGAt                                                 >  1:2090503/1‑150 (MQ=255)
                                                                                                       tttCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGAtt                                                <  2:1326033/150‑1 (MQ=255)
                                                                                                         tCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGt                                              <  1:2234328/150‑1 (MQ=255)
                                                                                                           cGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTc                                            <  1:4327867/150‑1 (MQ=255)
                                                                                                             cATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTcgc                                          >  1:260551/1‑150 (MQ=255)
                                                                                                                cgcgCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCa                                       <  1:1011069/150‑1 (MQ=255)
                                                                                                                   gcTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATaa                                    <  2:4073701/150‑1 (MQ=255)
                                                                                                                    cTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATaaa                                   <  2:4432880/150‑1 (MQ=255)
                                                                                                                    cTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATaaa                                   >  1:4895782/1‑150 (MQ=255)
                                                                                                                       atgaATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCg                                >  1:2498651/1‑150 (MQ=255)
                                                                                                                       atgaAGATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCg                                <  2:2018548/150‑1 (MQ=255)
                                                                                                                          aaTATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTga                             <  2:61841/150‑1 (MQ=255)
                                                                                                                            tatCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTgaga                           >  1:1969877/1‑150 (MQ=255)
                                                                                                                            tatCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATAGTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAAGGGTCACCGGTGTGGAGTGTTGGCGCATAAAGCGTgaga                           >  2:195041/1‑150 (MQ=255)
                                                                                                                             atCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGAc                          <  2:1455968/150‑1 (MQ=255)
                                                                                                                             atCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGAc                          <  2:4016401/150‑1 (MQ=255)
                                                                                                                                cAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAACGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGcc                       <  1:3023710/150‑1 (MQ=255)
                                                                                                                                   cACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGt                    <  2:3680781/150‑1 (MQ=255)
                                                                                                                                      cTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGAt                 >  1:3963094/1‑150 (MQ=255)
                                                                                                                                        cATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAg               <  2:2615217/150‑1 (MQ=255)
                                                                                                                                           ccaccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGc            >  2:3314843/1‑150 (MQ=255)
                                                                                                                                           ccaccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGc            >  1:4819130/1‑150 (MQ=255)
                                                                                                                                              ccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCgg         >  1:2446117/1‑150 (MQ=255)
                                                                                                                                              ccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCgg         >  1:3359568/1‑150 (MQ=255)
                                                                                                                                              ccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCgg         >  2:4365740/1‑150 (MQ=255)
                                                                                                                                              ccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCgg         <  2:1989878/150‑1 (MQ=255)
                                                                                                                                              ccaGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGGGATAGTGCCgg         >  2:781075/1‑150 (MQ=255)
                                                                                                                                                  agaAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGgtcgt     >  1:28713/1‑150 (MQ=255)
                                                                                                                                                   gaAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGgtcgtc    >  1:4336136/1‑150 (MQ=255)
                                                                                                                                                    aaGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTcc   >  1:1415626/1‑150 (MQ=255)
                                                                                                                                                    aaGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTcc   <  2:235853/150‑1 (MQ=255)
                                                                                                                                                    aaGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGGTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTcc   <  2:2108584/150‑1 (MQ=255)
                                                                                                                                                     aGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTccc  >  2:3211482/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTCAGGGTTAGCTTTTGCCACCCGCAGTGCTCCGGCAACCGCGCCGCCGGAGCTGACGCCACAGAATATTCCTTCCCGCACCGCCAGTTCGCGCATGGTGTTTTCCGCATCGCGCTGATGAATATCCAGCACCTCATCCACCAGAGAAGCGTTGAAAATCCCCGGCAGATATTCCGTAGGCCAGCGGCGAATGCCGGGAATGCTGCTGCCCTCTTCCGGTTGCAGGCCGACAATGGTCACCGGTTTGGATTGTTCGCGCATAAAGCGTGAGACGCCGGTGATAGTGCCGGTCGTCCC  >  CP009273/2532125‑2532423

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: