Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,740 C→A 5.2% intergenic (+168/‑272) prpB → / → prpC 2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,7400CA 5.2% 542.6 / 9.9 167intergenic (+168/‑272)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base C (56/100);  new base A (8/3);  total (64/103)
Fisher's exact test for biased strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

TTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTG  >  NC_000913/349592‑349888
                                                                                                                                                    |                                                                                                                                                    
ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCa                                                                                                                                                     <  2:2774456/150‑1 (MQ=255)
         cGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                                                                                            <  2:1944377/150‑1 (MQ=255)
          gATGCCTGATGCGTCGCTGCCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                                                                                           <  1:116760/150‑1 (MQ=255)
          gATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                                                                                           <  1:116939/150‑1 (MQ=255)
          gATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                                                                                           <  1:656658/150‑1 (MQ=255)
             gCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                                                                                                                        <  2:1555042/150‑1 (MQ=255)
              ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTc                                                                                                                                       <  1:2392444/150‑1 (MQ=255)
               cTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCg                                                                                                                                      >  2:777020/1‑150 (MQ=255)
               cTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGAAGCCTGATGCGCCGCTGACGAGACTTATAAGGCCTACGGGGTGAACTGAACTGTAGGTCg                                                                                                                                      >  2:175163/1‑150 (MQ=255)
                 gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                                                                                                    <  1:721066/150‑1 (MQ=255)
                 gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                                                                                                    <  1:1642681/150‑1 (MQ=255)
                 gATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                                                                                                    >  2:64790/1‑150 (MQ=255)
                      gTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                                               <  1:2789725/150‑1 (MQ=255)
                      gTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                                               <  1:1753619/150‑1 (MQ=255)
                         gCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGc                                                                                                                            <  2:1698001/150‑1 (MQ=255)
                         gCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGc                                                                                                                            >  1:452534/1‑150 (MQ=255)
                          cTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                                                                                                                           <  2:816975/150‑1 (MQ=255)
                             aCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg                                                                                                                        <  1:1794506/150‑1 (MQ=255)
                             aCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg                                                                                                                        <  1:1755422/150‑1 (MQ=255)
                               gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                                                                                      >  1:1131338/1‑150 (MQ=40)
                                cgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                                                                                                     <  2:864173/150‑1 (MQ=40)
                                cgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGt                                                                                                                     <  2:51225/150‑1 (MQ=37)
                                 gACTTATCCGGCCTACGAGGTGCACTGAACTGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACACCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    <  2:683924/150‑1 (MQ=255)
                                 gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    <  2:2606917/150‑1 (MQ=39)
                                 gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    >  1:582908/1‑150 (MQ=39)
                                 gACTGATCAGGCCTCCGCGGTGCACTGACCTGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGCTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTCCGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAGAGCGTc                                                                                                                    <  2:3059656/150‑1 (MQ=11)
                                  aCTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGTCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                                                                                   >  2:1730219/1‑150 (MQ=255)
                                      aTCAGGCCAACGAGGTGCACTGCACTGTAGGTCGGATAAGACGGAGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc                                                                                                               <  1:1284093/150‑1 (MQ=255)
                                       tCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                                                                                              <  2:2262572/150‑1 (MQ=39)
                                         aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                                                                                            <  2:2578302/150‑1 (MQ=39)
                                           gCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACa                                                                                                          <  2:568806/150‑1 (MQ=39)
                                             cTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                                                                                                        <  1:1955326/150‑1 (MQ=39)
                                                cGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACCACCGa                                                                                                     >  1:2840775/1‑150 (MQ=255)
                                                cGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGa                                                                                                     <  2:2657527/150‑1 (MQ=39)
                                                cGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTAGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGATCTGTAGGTCGGATAAGACTCATATCGTCGAATCCGActcccca                                                                                                     <  2:2657505/150‑7 (MQ=14)
                                                 gAGGTGCACTGAACTGTCGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                                                                                    <  1:1866107/150‑1 (MQ=38)
                                                  aGGTGCACTGAACTGGAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                                                                                                   <  2:1008208/150‑1 (MQ=38)
                                                  aGGTGAACTGAATTGTTGGTCGGATAAGAGGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGGTGACGTGACTTATCTGGCCTAGGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                                                                                                   <  1:2665978/150‑1 (MQ=11)
                                                   ggTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTTGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                  <  1:2213159/150‑1 (MQ=37)
                                                   ggTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                  <  2:1102580/150‑1 (MQ=37)
                                                   ggTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                  <  2:219829/150‑1 (MQ=37)
                                                     tGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                                                                                                <  1:583322/150‑1 (MQ=37)
                                                          tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATc                                                                                                               >  2:45505/1‑130 (MQ=38)
                                                          tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                                                                                              >  2:924775/1‑147 (MQ=37)
                                                          tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                                                                                              <  1:924775/147‑1 (MQ=37)
                                                          tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATGCCACATgc                                                                                           >  2:74185/1‑150 (MQ=255)
                                                            aaCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc                                                                                         <  2:2993866/150‑1 (MQ=37)
                                                                 tAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                                                                                    >  1:1733395/1‑150 (MQ=37)
                                                                 tAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTg                                                                                    >  1:1903251/1‑150 (MQ=37)
                                                                       ggATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                              <  1:1534325/150‑1 (MQ=37)
                                                                          tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt                                                                           >  1:1475093/1‑150 (MQ=37)
                                                                          tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt                                                                           >  2:2580013/1‑150 (MQ=37)
                                                                              aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                       >  2:2305890/1‑150 (MQ=37)
                                                                              aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                       <  2:386522/150‑1 (MQ=37)
                                                                              aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGAAGCGCCGCTGACGCGACTTATCa                                                                       >  2:1838723/1‑150 (MQ=255)
                                                                              aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGAAGCGCCGCGTACACGAATTATct                                                                       >  2:1839003/1‑149 (MQ=11)
                                                                                ggATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCCACGCCTCAGGCGCCCCTGACCCGACGTATCAgc                                                                     >  1:2671484/1‑149 (MQ=11)
                                                                                  aTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    >  2:1880673/1‑101 (MQ=17)
                                                                                  aTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    <  1:1880673/101‑1 (MQ=17)
                                                                                  aTGGCGTCGCATCCGACAACCGATGCCTGATGCGACGCTGACGTGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc                                                                   <  2:2886693/150‑1 (MQ=21)
                                                                                     gCGTCGCATCCGACACCCGATGCCTGATGCGCCGCTGAGGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATTGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTAc                                                                <  2:2851966/150‑1 (MQ=12)
                                                                                      cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACg                                                               <  1:1829418/150‑1 (MQ=32)
                                                                                       gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATGCCTGATGCGACGCTGACGCGACTACGAAGGCCTACgg                                                              >  2:1914868/1‑150 (MQ=11)
                                                                                        tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCTGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                                             <  2:1655956/150‑1 (MQ=17)
                                                                                        tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                                             <  2:1812896/150‑1 (MQ=32)
                                                                                        tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                                             <  2:1655928/150‑1 (MQ=32)
                                                                                          gCATCCGACAACGGGTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                           <  2:2452147/150‑1 (MQ=17)
                                                                                          gCATCCGACAACCGGTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                           <  2:1110964/150‑1 (MQ=17)
                                                                                          gCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCATGAATACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                           <  1:1538221/150‑1 (MQ=12)
                                                                                           cATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTg                                                          <  1:2008423/150‑1 (MQ=32)
                                                                                                gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                     >  2:138938/1‑150 (MQ=34)
                                                                                                 aCAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc                                                                                         >  1:2613940/1‑113 (MQ=32)
                                                                                                  cAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                                                   >  1:2658871/1‑150 (MQ=34)
                                                                                                  cAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGaa                                                   >  1:2541154/1‑150 (MQ=34)
                                                                                                    aCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATAAGGCCTACGGGGTGCTCTGAACt                                                 <  2:2360121/150‑1 (MQ=255)
                                                                                                          gCCTGCTGCGCCGCCGCCGTGACTTATCAGGCCTACGGGGTGCACTGCACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAggc                                           <  1:1080108/150‑2 (MQ=255)
                                                                                                          gCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                           <  1:357962/150‑1 (MQ=34)
                                                                                                          gCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                           <  1:1080229/150‑1 (MQ=34)
                                                                                                             tGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGGAGGTCgg                                        >  2:1290004/1‑150 (MQ=21)
                                                                                                             cgATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCTTCGCATCCGACTTCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                        <  2:2408669/149‑1 (MQ=255)
                                                                                                              gATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGACCTGTAGGTCGTATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                       <  2:939953/150‑1 (MQ=21)
                                                                                                              gATGCGCCGCTGAAGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTATGATAAGACGCATTGCGTCGCATCCGACACCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                       <  2:950704/150‑1 (MQ=11)
                                                                                                               aTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATTGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                      <  2:1131338/150‑1 (MQ=21)
                                                                                                                 gcgcCGCTTACTTTTCTTATCTGGCCTTGGGGGTGCACTTAACTGTAGGTCTGATAATACGCAGAGCGTCGGATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATaa                                    <  2:2416590/150‑1 (MQ=11)
                                                                                                                      gcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGc                               >  1:244753/1‑150 (MQ=34)
                                                                                                                       cTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                              >  1:3040951/1‑150 (MQ=34)
                                                                                                                       cTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                              >  1:620054/1‑150 (MQ=34)
                                                                                                                             tGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                        >  1:2279090/1‑150 (MQ=32)
                                                                                                                               aCTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                      <  2:452135/150‑1 (MQ=32)
                                                                                                                                 ttATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                    >  1:732018/1‑150 (MQ=32)
                                                                                                                                 ttATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                    >  2:1238640/1‑150 (MQ=32)
                                                                                                                                 ttATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                    >  2:1373638/1‑150 (MQ=32)
                                                                                                                                    tCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                 >  2:1428830/1‑150 (MQ=32)
                                                                                                                                    tCAGGCCGACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATATCGTCGCATCCGACATCCGATGCCTGATGCGCCGCTGACGCGAATTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                 <  1:330553/150‑1 (MQ=255)
                                                                                                                                      aGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa               >  2:615680/1‑150 (MQ=32)
                                                                                                                                      aGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa               >  1:453224/1‑150 (MQ=32)
                                                                                                                                      aGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa               >  1:378843/1‑150 (MQ=32)
                                                                                                                                      aGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa               >  1:2488766/1‑150 (MQ=32)
                                                                                                                                         ccTACGGGGTGCACTGAACTGTAGGACGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATATCGTCGCATCCGACaa            >  1:526750/1‑150 (MQ=17)
                                                                                                                                                   gCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCTTCCGACAACCGATGCCTg  >  1:851646/1‑150 (MQ=37)
                                                                                                                                                    |                                                                                                                                                    
TTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTG  >  NC_000913/349592‑349888

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: