Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,961,295 G→T 21.4% 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,2950GT21.4% 687.3 / 131.7 248R167S (CGC→AGC) ptsPfused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (89/106);  new base T (25/28);  total (114/134)
Fisher's exact test for biased strand distribution p-value = 8.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 253 total aligned reads displayed.

AGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGTT  >  CP009273/2961147‑2961444
                                                                                                                                                    |                                                                                                                                                     
aGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGa                                                                                                                                                      >  2:853253/1‑150 (MQ=255)
aGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGa                                                                                                                                                      >  1:3108669/1‑150 (MQ=255)
 gcgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGat                                                                                                                                                     <  2:1750272/150‑1 (MQ=255)
  cgcCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGata                                                                                                                                                    <  1:3353032/150‑1 (MQ=255)
        ggTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGccc                                                                                                                                              <  2:609205/150‑1 (MQ=255)
          tCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCaa                                                                                                                                            >  2:2006414/1‑150 (MQ=255)
            aGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAAc                                                                                                                                          >  1:1210333/1‑150 (MQ=255)
            aGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAAc                                                                                                                                          >  2:3855515/1‑150 (MQ=255)
             gTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACa                                                                                                                                         <  2:2748481/150‑1 (MQ=255)
                  tCGCGTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCa                                                                                                                                    >  2:500842/1‑150 (MQ=255)
                    gcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGt                                                                                                                                  <  1:270368/150‑1 (MQ=255)
                    gcgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGt                                                                                                                                  <  2:4577952/150‑1 (MQ=255)
                     cgTTCCAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTc                                                                                                                                 >  1:4011006/1‑150 (MQ=255)
                         ccAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACt                                                                                                                             <  1:511698/150‑1 (MQ=255)
                          cAGAGCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTg                                                                                                                            >  1:1030472/1‑150 (MQ=255)
                           agagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGc                                                                                                                           <  1:226329/150‑1 (MQ=255)
                           agagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGc                                                                                                                           >  1:1693253/1‑150 (MQ=255)
                            gagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCg                                                                                                                          <  2:4012881/150‑1 (MQ=255)
                            gagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCg                                                                                                                          <  2:1669774/150‑1 (MQ=255)
                            gagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCg                                                                                                                          <  1:114242/150‑1 (MQ=255)
                            gagCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCg                                                                                                                          >  1:3885961/1‑150 (MQ=255)
                             agCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGa                                                                                                                         <  2:2588100/150‑1 (MQ=255)
                             agCCGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCAGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGa                                                                                                                         >  2:1885876/1‑150 (MQ=255)
                               ccGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACt                                                                                                                       <  2:538562/150‑1 (MQ=255)
                                cGGATCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTg                                                                                                                      <  2:1739311/150‑1 (MQ=255)
                                   aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga                                                                                                                   >  2:362699/1‑150 (MQ=255)
                                   aTCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTgaga                                                                                                                   <  1:3748593/150‑1 (MQ=255)
                                    tCCAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaa                                                                                                                  >  2:4264879/1‑150 (MQ=255)
                                      cAGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaag                                                                                                                <  1:1592410/150‑1 (MQ=255)
                                       aGCGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGagaaaga                                                                                                               <  2:4545125/150‑1 (MQ=255)
                                         cGTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAAt                                                                                                             >  1:159466/1‑150 (MQ=255)
                                          gTTGATGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATa                                                                                                            >  2:2647184/1‑150 (MQ=255)
                                              aTGCCCGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTg                                                                                                        >  2:2144679/1‑150 (MQ=255)
                                               tGCCTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGc                                                                                                       >  1:2854072/1‑150 (MQ=255)
                                                  cTGATACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCAt                                                                                                    <  1:7128/150‑1 (MQ=255)
                                                      tACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTgg                                                                                                >  1:1063057/1‑150 (MQ=255)
                                                      tACACCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTgg                                                                                                <  1:301335/150‑1 (MQ=255)
                                                         acCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTg                                                                                             <  2:4116960/150‑1 (MQ=255)
                                                         acCTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTg                                                                                             <  1:4135989/150‑1 (MQ=255)
                                                          ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTgg                                                                                            <  1:3383835/150‑1 (MQ=255)
                                                          ccTGTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTgg                                                                                            >  1:2068063/1‑150 (MQ=255)
                                                             gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCaa                                                                                         >  1:4067742/1‑150 (MQ=255)
                                                             gTTCCATTAAAGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCaa                                                                                         <  1:3662915/150‑1 (MQ=255)
                                                                     aaaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCAcc                                                                                 >  2:749672/1‑150 (MQ=255)
                                                                     aaaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCAcc                                                                                 <  1:1905953/150‑1 (MQ=255)
                                                                      aaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGTCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCa                                                                                <  2:1604647/150‑1 (MQ=255)
                                                                      aaGGTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCa                                                                                <  2:2933308/150‑1 (MQ=255)
                                                                        ggTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCagg                                                                              >  1:3237314/1‑150 (MQ=255)
                                                                        ggTAACGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACCAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCagg                                                                              <  1:2836809/150‑1 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagga                                                                         >  2:915556/1‑150 (MQ=255)
                                                                             cGTGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCaggaagga                                                                         >  2:32344/1‑150 (MQ=255)
                                                                               tGGCATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAtt                                                                       <  2:4024655/150‑1 (MQ=255)
                                                                                  cATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAttctt                                                                    >  2:713328/1‑150 (MQ=255)
                                                                                  cATCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAttctt                                                                    <  1:800432/150‑1 (MQ=255)
                                                                                   aTCCTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGAttcttc                                                                   <  1:195777/150‑1 (MQ=255)
                                                                                      cTGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACt                                                                >  2:2731026/1‑150 (MQ=255)
                                                                                       tGCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACtt                                                               >  1:2476598/1‑150 (MQ=255)
                                                                                        gCCAGCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACttt                                                              <  1:1091690/150‑1 (MQ=255)
                                                                                            gCCTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTc                                                          <  1:633122/150‑1 (MQ=255)
                                                                                             ccTTCGGCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCa                                                         <  2:3306618/150‑1 (MQ=255)
                                                                                                   gCAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTgg                                                   <  1:428592/150‑1 (MQ=255)
                                                                                                    cAATCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGgc                                                  <  2:4599207/150‑1 (MQ=255)
                                                                                                     aaTCGCCACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGgcg                                                 >  2:2228725/1‑150 (MQ=255)
                                                                                                          ccACACCAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAAct                                            <  2:3730138/150‑1 (MQ=255)
                                                                                                               ccAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGt                                       <  2:1391521/150‑1 (MQ=255)
                                                                                                               ccAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGt                                       >  2:326364/1‑150 (MQ=255)
                                                                                                                cAGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGtt                                      <  1:2024833/150‑1 (MQ=255)
                                                                                                                 aGGTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTg                                     >  1:2307420/1‑150 (MQ=255)
                                                                                                                   gTGCTGCCGGTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCt                                   <  1:2939428/150‑1 (MQ=255)
                                                                                                                       tgcCGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTac                               >  2:3092839/1‑150 (MQ=255)
                                                                                                                         ccGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTacca                             <  2:524410/150‑1 (MQ=255)
                                                                                                                         ccGGTAATGCGCGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTacca                             <  2:3638428/150‑1 (MQ=255)
                                                                                                                           ggTAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTaccacc                           >  2:2955333/1‑150 (MQ=255)
                                                                                                                             tAATGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAg                         >  2:2562796/1‑150 (MQ=255)
                                                                                                                              aaTGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGt                        <  2:1722358/150‑1 (MQ=255)
                                                                                                                               aTGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTa                       <  1:565290/150‑1 (MQ=255)
                                                                                                                                tGCGCGGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTac                      <  1:75994/150‑1 (MQ=255)
                                                                                                                                 gcgcGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaca                     <  1:1619669/150‑1 (MQ=255)
                                                                                                                                 gcgcGGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTaca                     <  1:1885876/150‑1 (MQ=255)
                                                                                                                                  cgcgGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTacac                    >  1:999840/1‑150 (MQ=255)
                                                                                                                                  cgcgGATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTacac                    >  1:1347749/1‑150 (MQ=255)
                                                                                                                                   gcgGATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACAcc                   <  2:2303872/150‑1 (MQ=255)
                                                                                                                                     ggATTCGCGTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCaa                 >  1:1092286/1‑150 (MQ=255)
                                                                                                                                     ggATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCaa                 <  1:2285950/150‑1 (MQ=255)
                                                                                                                                      gATTCGCGTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAg                <  2:4402628/150‑1 (MQ=255)
                                                                                                                                          cgcgTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAAc            <  1:3481366/150‑1 (MQ=255)
                                                                                                                                            cgTCTGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTg          >  1:1813535/1‑150 (MQ=255)
                                                                                                                                             gTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTgg         >  1:4190157/1‑150 (MQ=255)
                                                                                                                                             gTCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTgg         >  2:510519/1‑150 (MQ=255)
                                                                                                                                              tCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGc        <  2:3466317/150‑1 (MQ=255)
                                                                                                                                              tCTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGc        <  2:1312899/150‑1 (MQ=255)
                                                                                                                                               cTGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCg       >  2:3932229/1‑150 (MQ=255)
                                                                                                                                                tGGCTATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGa      >  1:1726348/1‑150 (MQ=255)
                                                                                                                                                tGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGa      >  2:2081111/1‑150 (MQ=255)
                                                                                                                                                tGGCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGa      >  2:3618004/1‑150 (MQ=255)
                                                                                                                                                 ggCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGAc     <  2:3439400/150‑1 (MQ=255)
                                                                                                                                                  gCGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACg    >  2:2636079/1‑150 (MQ=255)
                                                                                                                                                   cGATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGt   >  2:4378030/1‑150 (MQ=255)
                                                                                                                                                    gATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGtt  >  2:2778125/1‑150 (MQ=255)
                                                                                                                                                    gATATTGCCCAAACAAGGCAGTCAACTGCGACTGAGAAAGAATAGCTGCCATCTGGGTGGCAAGCGTCACCAGGAAGGATTCTTCACTTTCGTCATACTGGCGCAACTCTCGTTGCTGTACCACCAGTACACCAAGCAACTGGCGACGtt  <  1:1251756/150‑1 (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: