Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,961,295 G→T 48.0% R167S (CGC→AGC)  ptsP ← fused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,961,2950GT48.0% 196.7 / inf 315R167S (CGC→AGC) ptsPfused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (91/73);  new base T (76/75);  total (167/148)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 317 total aligned reads displayed.

AGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTT  >  CP009273/2961147‑2961444
                                                                                                                                                    |                                                                                                                                                     
   gcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTatat                                                                                                                                                   >  1:2536658/1‑150 (MQ=255)
   gcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGatat                                                                                                                                                   >  2:658048/1‑150 (MQ=255)
   gcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGatat                                                                                                                                                   >  2:1519804/1‑150 (MQ=255)
   gcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGatat                                                                                                                                                   >  1:2306309/1‑150 (MQ=255)
    ccccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATAtt                                                                                                                                                  >  2:3664527/1‑150 (MQ=255)
    ccccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATAtt                                                                                                                                                  <  2:1688789/150‑1 (MQ=255)
     cccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTg                                                                                                                                                 <  2:2401190/150‑1 (MQ=255)
     cccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTg                                                                                                                                                 >  2:659288/1‑150 (MQ=255)
      ccGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGc                                                                                                                                                <  2:3318405/150‑1 (MQ=255)
        ggTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGccc                                                                                                                                              <  1:189848/150‑1 (MQ=255)
        ggTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                                                                                                                              >  2:208888/1‑150 (MQ=255)
               cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCTCGTCTGGCGATATTGCCCAAACAAg                                                                                                                                       >  2:3336601/1‑150 (MQ=255)
               cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAg                                                                                                                                       <  2:2124850/150‑1 (MQ=255)
               cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAg                                                                                                                                       >  2:1013193/1‑150 (MQ=255)
               cGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAg                                                                                                                                       >  1:790100/1‑150 (MQ=255)
                 ttCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGc                                                                                                                                     <  1:1775472/150‑1 (MQ=255)
                  tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCa                                                                                                                                    >  1:11310/1‑150 (MQ=255)
                  tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCa                                                                                                                                    >  1:4277701/1‑150 (MQ=255)
                  tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCa                                                                                                                                    <  2:3560714/150‑1 (MQ=255)
                  tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCa                                                                                                                                    >  1:3830560/1‑150 (MQ=255)
                   cgcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAg                                                                                                                                   <  2:4408047/150‑1 (MQ=255)
                        tCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAAc                                                                                                                              >  2:3350841/1‑150 (MQ=255)
                           agagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGc                                                                                                                           >  2:4437946/1‑150 (MQ=255)
                             agCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGa                                                                                                                         <  2:1503416/150‑1 (MQ=255)
                                   aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga                                                                                                                   <  2:1443466/150‑1 (MQ=255)
                                   aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga                                                                                                                   <  2:3774618/150‑1 (MQ=255)
                                      cAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaag                                                                                                                <  1:3240541/150‑1 (MQ=255)
                                      cAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaag                                                                                                                <  1:4489773/150‑1 (MQ=255)
                                         cGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAt                                                                                                             <  2:2556940/150‑1 (MQ=255)
                                         cGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAt                                                                                                             <  1:355814/150‑1 (MQ=255)
                                          gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATa                                                                                                            >  2:4395206/1‑150 (MQ=255)
                                                gCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGcc                                                                                                      <  1:127003/150‑1 (MQ=255)
                                                 ccTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCa                                                                                                     >  2:3301607/1‑150 (MQ=255)
                                                 ccTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCa                                                                                                     >  2:4245365/1‑150 (MQ=255)
                                                    gATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCt                                                                                                  <  2:1234642/150‑1 (MQ=255)
                                                    gATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCt                                                                                                  >  2:199370/1‑150 (MQ=255)
                                                    gATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCt                                                                                                  <  2:2186809/150‑1 (MQ=255)
                                                          ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTgg                                                                                            <  1:1664955/150‑1 (MQ=255)
                                                           cTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGc                                                                                           >  2:2847111/1‑150 (MQ=255)
                                                             gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAACCAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCaa                                                                                         <  1:2687571/150‑1 (MQ=255)
                                                              ttCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAg                                                                                        >  1:3125516/1‑150 (MQ=255)
                                                              ttCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAg                                                                                        >  1:543047/1‑150 (MQ=255)
                                                                 cATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGt                                                                                     <  1:1042390/150‑1 (MQ=255)
                                                                  aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTc                                                                                    <  1:4136838/150‑1 (MQ=255)
                                                                  aTTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTc                                                                                    <  1:4136780/150‑1 (MQ=255)
                                                                     aaaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCAcc                                                                                 <  2:2055728/150‑1 (MQ=255)
                                                                      aaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCa                                                                                >  2:1860706/1‑150 (MQ=255)
                                                                      aaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGCGGCAAGCGTCACCa                                                                                >  2:3178631/1‑150 (MQ=255)
                                                                      aaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCAGCCATCTGGGTGGCAAGCGACACCa                                                                                >  1:2469151/1‑150 (MQ=255)
                                                                       aGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCag                                                                               >  1:3228540/1‑150 (MQ=255)
                                                                         gTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCagga                                                                             <  1:1678290/150‑1 (MQ=255)
                                                                          tAACGTGGCCTCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                            >  2:1782809/1‑150 (MQ=255)
                                                                          tAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCCGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaa                                                                            <  2:3830452/150‑1 (MQ=255)
                                                                            aCGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagg                                                                          <  1:28914/150‑1 (MQ=255)
                                                                            aCGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagg                                                                          <  1:676878/150‑1 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagga                                                                         <  2:4569651/150‑1 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagga                                                                         <  1:1162284/150‑1 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagga                                                                         <  1:2212750/150‑1 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagga                                                                         >  2:4265307/1‑150 (MQ=255)
                                                                              gTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAt                                                                        >  2:602526/1‑150 (MQ=255)
                                                                                   aTCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAttcttc                                                                   <  1:3073460/150‑1 (MQ=255)
                                                                                       tGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACtt                                                               <  1:4001519/150‑1 (MQ=255)
                                                                                         ccAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTc                                                             >  1:209068/1‑150 (MQ=255)
                                                                                         ccAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTc                                                             >  1:2346161/1‑150 (MQ=255)
                                                                                          cAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACtttag                                                            <  1:424090/150‑3 (MQ=255)
                                                                                            gCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTc                                                          <  1:4437946/150‑1 (MQ=255)
                                                                                            gCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTc                                                          >  2:2972532/1‑150 (MQ=255)
                                                                                             ccTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGAGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGTCTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCa                                                         <  1:1724353/150‑1 (MQ=255)
                                                                                              cTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCAt                                                        <  2:13521/150‑1 (MQ=255)
                                                                                                tCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATAc                                                      <  1:954711/150‑1 (MQ=255)
                                                                                                 cGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACt                                                     <  2:3100468/150‑1 (MQ=255)
                                                                                                 cGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACt                                                     >  1:621200/1‑150 (MQ=255)
                                                                                                    cAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGgc                                                  >  1:3652285/1‑150 (MQ=255)
                                                                                                    cAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGgc                                                  >  1:3614366/1‑150 (MQ=255)
                                                                                                     aaTCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGgcg                                                 <  2:1629195/150‑1 (MQ=255)
                                                                                                       tCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCa                                               <  2:385742/150‑1 (MQ=255)
                                                                                                         gCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAAc                                             >  1:1337219/1‑150 (MQ=255)
                                                                                                         gCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAAc                                             >  2:4031441/1‑150 (MQ=255)
                                                                                                            acacCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAActct                                          <  1:2621766/150‑1 (MQ=255)
                                                                                                            acacCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAActct                                          >  1:3369865/1‑150 (MQ=255)
                                                                                                            acacCAGGTGCTGCCGGTAATGCGAGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAActct                                          <  1:2037486/150‑1 (MQ=255)
                                                                                                             cacCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACtctc                                         >  1:4560541/1‑150 (MQ=255)
                                                                                                               ccAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGt                                       <  1:3283238/150‑1 (MQ=255)
                                                                                                                cAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGtt                                      >  2:769646/1‑150 (MQ=255)
                                                                                                                 aGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTg                                     <  2:3044354/150‑1 (MQ=255)
                                                                                                                  ggTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGc                                    >  2:1910931/1‑150 (MQ=255)
                                                                                                                  ggTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGc                                    >  2:4241660/1‑150 (MQ=255)
                                                                                                                   gTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCt                                   <  2:2920874/150‑1 (MQ=255)
                                                                                                                       tgcCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                               >  1:3434825/1‑150 (MQ=255)
                                                                                                                           ggTAATGCGCGGATTCGCGTCTGGCTATAATGCCCAAACAAGGCAGACAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAATCTTCACAATCGTCATACTGGCGCAACTCTCGTTGCTGTaccacc                           >  1:375958/1‑150 (MQ=255)
                                                                                                                           ggTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTaccacc                           >  1:303735/1‑150 (MQ=255)
                                                                                                                            gTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTAccacca                          >  1:879681/1‑150 (MQ=255)
                                                                                                                            gTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTAccacca                          <  2:2225170/150‑1 (MQ=255)
                                                                                                                              aaTGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGt                        <  2:11310/150‑1 (MQ=255)
                                                                                                                               aTGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                       >  1:46283/1‑150 (MQ=255)
                                                                                                                               aTGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                       >  1:1527998/1‑150 (MQ=255)
                                                                                                                               aTGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                       >  1:1594487/1‑150 (MQ=255)
                                                                                                                                    cgGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCa                  >  1:2362401/1‑150 (MQ=255)
                                                                                                                                     ggATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCaa                 <  2:2306303/150‑1 (MQ=255)
                                                                                                                                                 ggCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCCGTACACCAAGCAACTGGCGAc     >  1:2798580/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTT  >  CP009273/2961147‑2961444

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: