Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 345,384 G→C 5.1% intergenic (+101/‑339) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273345,3840GC 5.1% 1129.5 / 27.7 323intergenic (+101/‑339)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (132/172);  new base C (15/4);  total (147/176)
Fisher's exact test for biased strand distribution p-value = 3.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

GTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGG  >  CP009273/345235‑345532
                                                                                                                                                     |                                                                                                                                                    
  aCGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGAt                                                                                                                                                    <  2:3557881/150‑1 (MQ=255)
   cGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTTAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATg                                                                                                                                                   <  2:4053938/150‑1 (MQ=255)
      agagAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                                                                                                                >  2:4578694/1‑150 (MQ=255)
        agaAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAATACGGATGGCGTc                                                                                                                                              >  1:984101/1‑150 (MQ=255)
           aGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACGGAACTGTAGGTCGGATAAGACGGATGGCGTCGCa                                                                                                                                           >  1:2493924/1‑150 (MQ=255)
              tCGACAACCTGTTTGCCCGTAGCCTGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATcc                                                                                                                                        <  2:3027195/150‑1 (MQ=255)
                gACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGa                                                                                                                                      >  2:3982091/1‑150 (MQ=255)
                cacacCCGGTTTGCCCTTAGCCAGTTTAAAAAAAACACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGa                                                                                                                                      <  1:3708185/149‑1 (MQ=21)
                 aCAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATAACGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCCTCCGAc                                                                                                                                     >  1:4401190/1‑150 (MQ=255)
                      cTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAATACGGATGGCGTCGCATCCGACCACCg                                                                                                                                >  1:1291033/1‑150 (MQ=255)
                      cTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAAATGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCg                                                                                                                                >  1:3110092/1‑150 (MQ=255)
                      cTGTTTGCCCGTACCCAGGTTAAATAAAACACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCg                                                                                                                                <  1:1890299/150‑1 (MQ=255)
                       tgtgcgCCCGTAGCCCGGTGAAAAAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGa                                                                                                                               <  1:2620926/145‑1 (MQ=255)
                       ctcgtGCCCGTACCCAGGTGAAATAAAAATCGCCCGTTGATGGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcaat                                                                                                                               <  2:740839/146‑4 (MQ=14)
                        gTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGAt                                                                                                                              <  1:3935640/150‑1 (MQ=255)
                         tttGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATg                                                                                                                             <  2:1328414/150‑1 (MQ=255)
                          ttGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGc                                                                                                                            <  2:1885505/150‑1 (MQ=255)
                            gCCCGTAGCCAGTTTAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCt                                                                                                                          <  2:105203/150‑1 (MQ=255)
                                  aGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGTTGCGTCGCTTACGCCAACTATCAGGCCGAACAGGTGAGCGGAACTGAACTACGGTGAAGTCGCATCCCGTCCCATGCGACCACCGATGCCGGATgcg                                                                                                                    >  1:4699909/1‑150 (MQ=255)
                                  aGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACAGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATgcg                                                                                                                    >  1:774450/1‑150 (MQ=255)
                                    ccAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCCAAGCCAGACGcgcc                                                                                                                  >  2:3439331/1‑150 (MQ=255)
                                      aGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGACGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                                                                                                                >  2:4038466/1‑150 (MQ=255)
                                       ggTGATATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACAGTACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt                                                                                                               >  1:1969115/1‑150 (MQ=255)
                                       ggTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATAAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATTCCTGATGCGCCGCt                                                                                                               >  2:4298266/1‑150 (MQ=255)
                                        gTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGCGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                                                                                                              <  2:270895/150‑1 (MQ=255)
                                         tGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACATATAAGGCCAACGAGGTGCCCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGCTGCGCCGCTGa                                                                                                             >  1:3000999/1‑150 (MQ=18)
                                          gAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCAGATGCGTCGCTGACGCGACTTATAAGGCCAACGAGGTGCACGGAACATTAGGTCGGATAAGACCGCTGGCGTCGCATCCCACCACCGATGCCTGATGCGCCGCTGAc                                                                                                            >  1:2873709/1‑150 (MQ=18)
                                           aaaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACg                                                                                                           <  2:4147001/150‑1 (MQ=255)
                                           aaaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACATATCAGGCCTACGAGGTGCACCGAACTGTAGGTCGGATAATAAGGATTGCGTCGCATCCCACAACCGATGCCAGGTTCGCCCCTGAcc                                                                                                           >  1:1189953/1‑149 (MQ=255)
                                            aaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGAAGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCAGATGCGCCGCTGACGt                                                                                                          >  1:329677/1‑150 (MQ=255)
                                              tAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGGCGCTGACGCGACTTATCAGGCCTACGCGGTGCCCTGAACTGTAGGTCGGATAAGACGGTTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGa                                                                                                        >  1:1363508/1‑150 (MQ=17)
                                               aaaaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGAc                                                                                                       <  1:3559396/150‑1 (MQ=255)
                                                  aaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                                              >  2:158760/1‑108 (MQ=14)
                                                  aaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                                              <  1:158760/108‑1 (MQ=14)
                                                     cGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                                                 <  2:1688237/150‑1 (MQ=255)
                                                       cccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAgg                                                                                               <  2:3533474/150‑1 (MQ=255)
                                                         cGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCAGACGCGACTTATCAGGCCTACGAGGTGCACAGAACTGTAGGTCGGATAAGACGGATGGCGTCGCAACCGACAACCGATGCCTGATGCGCCGCTGACGTGAATTAACAGGcc                                                                                             >  2:4074896/1‑150 (MQ=17)
                                                          gTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCAGACGTGACTTATCAGGCCt                                                                                            >  1:1065726/1‑150 (MQ=255)
                                                            tGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGTTAATACGGATGGCGTCGCATACGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTAc                                                                                          >  1:4251908/1‑150 (MQ=255)
                                                            tGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTAc                                                                                          <  2:2132769/150‑1 (MQ=255)
                                                             gATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGAAGCCTGATTCGCCGCTGACGTGACATATCAGGCCTaag                                                                                         >  2:3266150/1‑148 (MQ=255)
                                                                  tATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTg                                                                                    <  2:960236/150‑1 (MQ=255)
                                                                    ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGGCTTATCAGGCCTACGGGGTGCa                                                                                  <  1:4382617/150‑1 (MQ=255)
                                                                    ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCa                                                                                  <  1:4369247/150‑1 (MQ=255)
                                                                      cGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCACTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                                <  1:2765374/150‑1 (MQ=255)
                                                                         cAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                             <  2:1947429/150‑1 (MQ=255)
                                                                         cAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACAGATTGCGTCGCATCCGACAACCCAATCCTGAAGCGCCGCTGACGAGACATATCAGGCCTACGGGGTGCACTGaa                                                                             >  1:1930224/1‑150 (MQ=14)
                                                                            ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTg                                                                          >  2:4433720/1‑150 (MQ=255)
                                                                               aTGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                       <  1:1004097/150‑1 (MQ=255)
                                                                               aTGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                       <  2:2878410/150‑1 (MQ=255)
                                                                               aTGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                       >  1:4158159/1‑150 (MQ=255)
                                                                                  ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTc                                                                    <  2:2072818/150‑1 (MQ=255)
                                                                                  ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTc                                                                    <  1:494457/150‑1 (MQ=255)
                                                                                   cTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                                   <  1:598149/150‑1 (MQ=255)
                                                                                       tGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa                                                               <  1:3355668/150‑1 (MQ=255)
                                                                                         cGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                             <  1:71056/150‑1 (MQ=255)
                                                                                           tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                           >  1:2544662/1‑150 (MQ=255)
                                                                                            cGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATCAGACg                                                          >  2:895842/1‑150 (MQ=255)
                                                                                              cTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGCCAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                                                        <  1:633565/150‑1 (MQ=255)
                                                                                                gACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                                      >  1:2223401/1‑150 (MQ=37)
                                                                                                   gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                   <  1:4716780/150‑1 (MQ=40)
                                                                                                   gcgACTGATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                   <  1:3451116/150‑1 (MQ=255)
                                                                                                    cgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                                  >  1:3979571/1‑150 (MQ=40)
                                                                                                    cgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                                  >  2:2335541/1‑150 (MQ=25)
                                                                                                     gACTTATCAGGCCTACGCGGTGCACTGAACAGTAGGTAGGATAATACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                 <  1:2287109/150‑1 (MQ=11)
                                                                                                     gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                 <  1:1915296/150‑1 (MQ=39)
                                                                                                     gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                 <  1:4179215/150‑1 (MQ=39)
                                                                                                      aCTTATCATGCCTACGAGGTGCACTGAACTGTAGTTCGGATAAGACGGATGGCTTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                <  2:225417/150‑1 (MQ=255)
                                                                                                        ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                                              >  2:1823508/1‑150 (MQ=39)
                                                                                                           tCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                           >  2:1280206/1‑150 (MQ=39)
                                                                                                             aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                         >  1:2025236/1‑150 (MQ=37)
                                                                                                                ccTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGAGGGCGTCGCATCCGACATCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa                                      <  2:640099/150‑1 (MQ=255)
                                                                                                                  tACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCTTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                    <  2:4225302/150‑1 (MQ=255)
                                                                                                                     gAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                 <  2:2579849/150‑1 (MQ=38)
                                                                                                                     gAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACAGTAGGTCGGATAAGACGCATAGCGTCTCATCCGACAACAGag                                 >  2:314034/1‑149 (MQ=255)
                                                                                                                      aGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                                <  1:3218645/150‑1 (MQ=38)
                                                                                                                         tGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                             >  2:1423035/1‑150 (MQ=37)
                                                                                                                          gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                            <  1:349128/150‑1 (MQ=37)
                                                                                                                               gAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAGgcg                       >  2:4561607/1‑150 (MQ=37)
                                                                                                                                aaCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc                      >  1:3583576/1‑150 (MQ=37)
                                                                                                                                 aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc                     <  2:3818614/150‑1 (MQ=37)
                                                                                                                                  cTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg                    <  2:1969115/150‑1 (MQ=37)
                                                                                                                                   tGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                   >  1:1829786/1‑150 (MQ=37)
                                                                                                                                    gTAGGTCGTATAAGACGGATGGCGTCGCATTCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt                  <  2:1164627/150‑1 (MQ=255)
                                                                                                                                    gTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAAca                                                                                                                                 >  1:3803632/1‑38 (MQ=255)
                                                                                                                                     tAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                 <  1:3638375/150‑1 (MQ=37)
                                                                                                                                     tAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                 >  2:4686901/1‑150 (MQ=37)
                                                                                                                                     tAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                 >  2:1740838/1‑150 (MQ=37)
                                                                                                                                       ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc               <  2:426940/150‑1 (MQ=37)
                                                                                                                                           ggATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa           >  2:1493363/1‑150 (MQ=37)
                                                                                                                                            taTAAGACGGCTGGGGTCGCATCCGACACCCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc          <  2:2291483/149‑1 (MQ=255)
                                                                                                                                            gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc          >  1:4520655/1‑150 (MQ=37)
                                                                                                                                              tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt        <  1:3078533/150‑1 (MQ=37)
                                                                                                                                              tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt        <  1:1646370/150‑1 (MQ=37)
                                                                                                                                              tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt        <  1:4442631/150‑1 (MQ=37)
                                                                                                                                               aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTa       >  1:3254731/1‑150 (MQ=37)
                                                                                                                                               aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCAGCGGCCCCGCTCACGCGACtga       >  2:3347337/1‑148 (MQ=11)
                                                                                                                                               aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACt                                                                                                      >  1:387680/1‑55 (MQ=255)
                                                                                                                                                aGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt      >  2:3757903/1‑150 (MQ=37)
                                                                                                                                                 gACGGATGGGGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc     <  2:3707013/150‑1 (MQ=37)
                                                                                                                                                     |                                                                                                                                                    
GTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGG  >  CP009273/345235‑345532

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: