Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 345,384 G→C 5.4% 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.4% 595.3 / 21.4 179intergenic (+101/‑339)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (60/109);  new base C (6/4);  total (66/113)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

ACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGG  >  CP009273/345237‑345532
                                                                                                                                                   |                                                                                                                                                    
 cGAAGAGACGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGAt                                                                                                                                                    <  2:1234789/149‑1 (MQ=255)
  gAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATg                                                                                                                                                   <  2:3480438/149‑1 (MQ=255)
   aaGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATgg                                                                                                                                                  >  1:3359708/1‑149 (MQ=255)
    agagACGCCCGACAACCTGTTTGCCCGAAGCCTGGTGAAAAAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGc                                                                                                                                                 <  2:432171/149‑1 (MQ=25)
    agagAAGCTCGACACCCTGTTGGCCCGTAGCCCGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGc                                                                                                                                                 <  2:431909/149‑1 (MQ=255)
          gCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCa                                                                                                                                           >  1:290702/1‑149 (MQ=255)
          gCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGGTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCa                                                                                                                                           >  1:290244/1‑149 (MQ=255)
             cGACAACCTTTTTGCCCGTCGCCAGGATAAAAAAAACACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATcc                                                                                                                                        <  1:1464149/149‑1 (MQ=21)
                  cccTGTTTGCCCGTCGCCAGGTTAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACaa                                                                                                                                   <  2:958820/148‑1 (MQ=255)
                   ccTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGGCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAAc                                                                                                                                  >  1:3436537/1‑149 (MQ=255)
                      gTTTGCCGTTAGCCTGGTTAATTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGa                                                                                                                               <  1:4328008/149‑1 (MQ=21)
                         tGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGc                                                                                                                            <  2:338042/149‑1 (MQ=255)
                          gCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGcc                                                                                                                           <  2:2041830/149‑1 (MQ=255)
                           cccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCt                                                                                                                          <  2:2615134/149‑1 (MQ=255)
                           cccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCt                                                                                                                          <  2:355127/149‑1 (MQ=255)
                            ccGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTg                                                                                                                         >  1:1015026/1‑149 (MQ=255)
                                 gCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACGTATAAGGCCGAAGAGGTGCACTGAACAGTAGGTCGGATAATACGGATGGCGTCGCAACCGACAACCGATGCCTGATgcg                                                                                                                    >  1:2172809/1‑149 (MQ=21)
                                  ccAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATTGCGTCGCATCCGACAACCGATGCCTGCTgcgc                                                                                                                   >  1:2129846/1‑149 (MQ=255)
                                    aGGTGAAATAAAAAACGCCCGTTGATTGTCTTCGACAGCCGATGCCTGATGCGTCGCTGGCGCGGCTTTTCAGGCCTGCGTGGTGCCCTGAACTGTAGGTGGGATAAGACGGATGGCGTCGCATCCGACACCCGGTGCCTGAGGcgccg                                                                                                                 >  1:2959375/1‑149 (MQ=255)
                                    aGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGcgccg                                                                                                                 >  1:524883/1‑149 (MQ=255)
                                     ggTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                                                                                                                >  1:2129591/1‑149 (MQ=255)
                                      gTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATTCGTCGCGGACGCGACTTATAAGGCCAACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTAGCATCCGACAACCGAAGCCTGATGCGCCGCt                                                                                                               >  1:1183004/1‑149 (MQ=18)
                                      gTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCAGTTGCGTCGCTGACGCGACATATAAGGCCTACGAGGTGCACTGAACTGTATGTCGGATAAGAAGGTTGGCGTCGCATACGACAACCGATGACTGATTCGGCGCt                                                                                                               >  1:1182904/1‑149 (MQ=255)
                                       tGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGAAGCGCCGCTg                                                                                                              >  2:3089652/1‑149 (MQ=255)
                                        taaaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                                             <  1:1412266/148‑1 (MQ=255)
                                        gAAATAAAACACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                                             <  1:2314070/149‑1 (MQ=255)
                                        gAAATAAAAAACGCCCGTTGATTGTATGCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                                             <  1:2313826/149‑1 (MQ=255)
                                          aaTAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACg                                                                                                           >  2:4178348/1‑149 (MQ=255)
                                            tAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGGCGCTGACGCGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGAGAAGACGCATGGCGTCGCGACCGGCAACCGATGCCTGATGCGCCGCTGACGTg                                                                                                         >  2:1753470/1‑149 (MQ=17)
                                               aacaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACt                                                                                                      <  1:909636/146‑1 (MQ=255)
                                               aaaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACt                                                                                                      >  2:4218673/1‑149 (MQ=255)
                                                aaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGTTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCAGATGCGCCGCGGACGTGACtt                                                                                                     >  1:3476676/1‑149 (MQ=255)
                                                 aaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTa                                                                                                    <  2:2756991/149‑1 (MQ=255)
                                                   cGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                                                                  <  1:4218673/149‑1 (MQ=255)
                                                      ccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTAACAGGCCAACGACGTGCACAGAACTGTAGGTCGGATAAGACGGATGGCGTCGCAACCCACAACCGAAGCCAGATGCGCCCCAGACGAGACTAATCAgg                                                                                               >  1:3151989/1‑149 (MQ=17)
                                                       cGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGc                                                                                              <  2:4144348/149‑1 (MQ=255)
                                                       cGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACGTATCAGGCCTACGAGGTGCACAGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCCAAGCCAGAAGCCCCACTGACCAGACTTAGCAGGc                                                                                              >  1:950352/1‑149 (MQ=17)
                                                        gTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgc                                                                                                                     >  2:3461387/1‑125 (MQ=255)
                                                        gTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGcc                                                                                             >  2:2186437/1‑149 (MQ=255)
                                                         ttGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCt                                                                                            <  2:3711158/149‑1 (MQ=255)
                                                          tGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTa                                                                                           <  2:2129846/149‑1 (MQ=255)
                                                           gATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTAc                                                                                          <  2:2162210/149‑1 (MQ=255)
                                                           gATTGTATTCGACAGCCGATGCCTGATGCGTCGCAGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTAc                                                                                          <  2:5330/149‑1 (MQ=255)
                                                                tATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGt                                                                                     <  1:2921362/149‑1 (MQ=255)
                                                                 aTTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTg                                                                                    <  2:3730420/149‑1 (MQ=255)
                                                                      aCAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTg                                                                               <  2:337330/149‑1 (MQ=255)
                                                                      aCAGCCGATGCCTGATGCGACGCTGACGAGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGCCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                                               <  1:1950722/149‑1 (MQ=255)
                                                                        aGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                             <  2:741627/149‑1 (MQ=255)
                                                                         gCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTTATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                            <  2:898526/149‑1 (MQ=255)
                                                                         gCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                            <  1:955394/149‑1 (MQ=255)
                                                                         gCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                            <  2:3585121/149‑1 (MQ=255)
                                                                         gCCGATGCCTGAGGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                            <  2:2645649/149‑1 (MQ=255)
                                                                         gCCGATGCCAGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGAAGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaactgaac                                                                            >  1:3403305/1‑149 (MQ=255)
                                                                          ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACt                                                                           <  2:2377312/149‑1 (MQ=255)
                                                                             aTGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                        >  1:3819648/1‑149 (MQ=255)
                                                                              tGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                       <  2:2481619/149‑1 (MQ=255)
                                                                              tGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                       >  1:1068602/1‑149 (MQ=255)
                                                                              tGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                       >  1:1069287/1‑149 (MQ=255)
                                                                                ccTGCTGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                                                     <  2:2604536/149‑1 (MQ=255)
                                                                                ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                                                     <  1:4178685/149‑1 (MQ=255)
                                                                                   gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                                  <  2:615981/149‑1 (MQ=255)
                                                                                   gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                                  <  1:1550812/149‑1 (MQ=255)
                                                                                   gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGt                                                                         <  2:3591517/142‑1 (MQ=21)
                                                                                   gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGt                                                                         >  1:3591517/1‑142 (MQ=21)
                                                                                   gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTAGCAACCGACAACCGAAGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                                  >  1:3591334/1‑149 (MQ=255)
                                                                                    aTGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                                 <  2:3874334/149‑1 (MQ=255)
                                                                                      gCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa                                                               <  2:1806976/149‑1 (MQ=255)
                                                                                      gCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa                                                               <  2:2975072/149‑1 (MQ=255)
                                                                                          cGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGAc                                                           <  1:86713/149‑1 (MQ=17)
                                                                                             tGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                                                        <  1:180312/149‑1 (MQ=255)
                                                                                              gACGCGACTGATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAt                                                       <  1:2304609/149‑1 (MQ=255)
                                                                                               aCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                                      <  2:290244/149‑1 (MQ=255)
                                                                                               aCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                                      <  2:290453/149‑1 (MQ=255)
                                                                                                 gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGc                                                    <  2:1651805/149‑1 (MQ=40)
                                                                                                    aCTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                 >  2:112224/1‑149 (MQ=17)
                                                                                                                tACGAGGTTCACTGTACTGTAGTTCGTATAAGACGGATGGAGTCGCATACGACAACCGATGCCTGATGCGCCTCTGAAGTGAATTATCAGGCCTACGGTGTGCACTGAAATGTAGTTATGATAAGAGGCATAGCGTCGCATCCTACAAc                                     <  2:3098456/149‑1 (MQ=14)
                                                                                                                 aCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATTGCGTCGCATCCGACAAcc                                    <  1:1143540/149‑1 (MQ=255)
                                                                                                                 aCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                    >  1:1358134/1‑149 (MQ=39)
                                                                                                                 aCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAca                                    >  1:3411461/1‑148 (MQ=39)
                                                                                                                   gAGGGGCACTGAACGGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGa                                  <  2:1851337/149‑1 (MQ=255)
                                                                                                                          acagAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATGCCAGa                           >  2:4235995/4‑149 (MQ=255)
                                                                                                                          aCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                           >  1:1919067/1‑149 (MQ=37)
                                                                                                                          aCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGACTGa                           >  1:4093858/1‑149 (MQ=17)
                                                                                                                          aCTGAACTGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                           <  2:4134412/149‑1 (MQ=37)
                                                                                                                           cTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                          >  2:1519651/1‑149 (MQ=37)
                                                                                                                            tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGTAGCCTGagg                         >  1:1518783/1‑147 (MQ=255)
                                                                                                                             gcaCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACATCCGATGCCTGATGCGACGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATTGCGTCGCATCCGACAACCGATGCCTGATgc                        <  2:1726709/147‑1 (MQ=255)
                                                                                                                                  gTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                   <  2:1936894/149‑1 (MQ=37)
                                                                                                                                   tAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCt                  <  2:1595022/149‑1 (MQ=37)
                                                                                                                                    aGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                 >  1:1216807/1‑149 (MQ=37)
                                                                                                                                    aGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                 <  2:834650/149‑1 (MQ=37)
                                                                                                                                     ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                <  1:2524849/149‑1 (MQ=37)
                                                                                                                                     ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                <  1:215743/149‑1 (MQ=37)
                                                                                                                                     ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                <  1:3718615/149‑1 (MQ=37)
                                                                                                                                          gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa           <  1:2631280/149‑1 (MQ=37)
                                                                                                                                           aTAAGAGGGAGTGCTTCGCATCCGACACCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc          <  2:1183004/149‑1 (MQ=11)
                                                                                                                                            tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGATGCGCCGCTGACTCGACt         >  2:878175/1‑149 (MQ=255)
                                                                                                                                             aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt        <  2:508536/149‑1 (MQ=37)
                                                                                                                                             aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAATCCAGAAGCGCCGCTGACGCGACta        >  2:1432074/1‑148 (MQ=255)
                                                                                                                                                   gATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg  >  2:3282178/1‑149 (MQ=35)
                                                                                                                                                   |                                                                                                                                                    
ACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGG  >  CP009273/345237‑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: