Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,769 A→T 6.1% intergenic (+197/‑243) 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,7690AT 6.1% 262.9 / 3.8 100intergenic (+197/‑243)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (39/50);  new base T (3/4);  total (42/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCT  >  NC_000913/349621‑349917
                                                                                                                                                    |                                                                                                                                                    
aCGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg                                                                                                                                                     <  1:2228820/150‑1 (MQ=255)
 cgcgACTTATCAGGCCTACGAGGTGCACTGACGTGTAGGTCTGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGc                                                                                                                                                    <  1:524253/150‑1 (MQ=255)
  gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                                                                                                                   <  1:1959072/150‑1 (MQ=40)
     aCTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                                                                                                                <  1:2230033/150‑1 (MQ=39)
       ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                                                                                                                                              >  1:579028/1‑150 (MQ=39)
          tCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                                                                                                                           <  2:866880/150‑1 (MQ=39)
          tCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTAGGATAAGACGCATAGCGTCGCATcc                                                                                                                                           <  1:1579222/150‑1 (MQ=255)
          tCAGGCCTACGAGGGGCACTGAACTGTAGGTCGGATAAGACGGATGGCGGCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                                                                                                                           <  2:1983452/150‑1 (MQ=255)
           cAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                                                                                                                          >  1:912201/1‑150 (MQ=39)
           cAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGAGACTTATCAGGCCTACGGGGTGCACTGAACTGTCGGGCGGATAGGACGCGTAGCGTCGCCTCCg                                                                                                                                          >  1:1683880/1‑150 (MQ=11)
           cAGGCCTACGAGGGGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                                                                                                                          <  2:2230121/150‑1 (MQ=255)
             ggCCTACGAGGTGCACTGAACTGTAGGTGGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                                                                                                                             <  1:1912191/145‑1 (MQ=255)
             ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATGCGACGCTATGCGTCTTATCCGAc                                                                                                                                        <  1:1077561/150‑1 (MQ=14)
             ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                                                        >  1:1350857/1‑150 (MQ=39)
             ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                                                        <  1:2164744/150‑1 (MQ=39)
             ggCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGACCTACGGGGTGAACTGAAATCTAGGTAGGATAAGACGCATAGCGTCGAATCCGaa                                                                                                                                        >  2:2323114/1‑149 (MQ=11)
             ggCCTACGAGGTGCACTGAACTGTAGGTAGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                                                                                                                             >  2:1912191/1‑145 (MQ=255)
               ccTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATGCGACGCTATGCGTCTTATCCGAcct                                                                                                                                      <  1:2190110/150‑3 (MQ=14)
                 taggAGGTGCTGTGATGTGTAGGTCTGATAAGACGGATGGCTTCGCATCTTACAACCGATGGCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGGACTGAACTGTAGGTCGGATAAGTCGCATAGCGTCGCATCCGACAAcc                                                                                                                                    <  2:912201/147‑1 (MQ=11)
                 tACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                                                                                                                    <  1:1962120/150‑1 (MQ=39)
                     aGGTGCACTGAACTGTAGGTCTGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTGGCATCCGACAACCGATg                                                                                                                                <  1:693667/150‑1 (MQ=255)
                     aGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                                                                                                                                <  1:74166/150‑1 (MQ=38)
                     aGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                                                                                                                                <  2:951553/150‑1 (MQ=38)
                      ggTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                                               <  2:1598610/150‑1 (MQ=37)
                      ggTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                                               <  2:353523/150‑1 (MQ=37)
                      ggTGCACTGAACTGTAGGTCGCATAAGACGGATGGCGTCGCATCCGACAACCTATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                                               <  1:2079803/150‑1 (MQ=255)
                       gTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                                                                                                                              <  2:1482228/150‑1 (MQ=37)
                         gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                                                                                                                            <  1:292181/150‑1 (MQ=37)
                         gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                                                                                                                            >  1:103325/1‑150 (MQ=37)
                          cACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                                                                                                                           >  2:1495937/1‑150 (MQ=37)
                           aCTGAACTGTAGGTCGTATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                                                                                                                          <  1:2420280/150‑1 (MQ=37)
                           aCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                                                                                                                          <  2:314626/150‑1 (MQ=37)
                            cTGAACTGTAGGTCGTATAAGACGGATGGCGTCGCATCCGACAACAGAATCCTGATGCGAAGCTCAAGGGACTTATCAGGCCTAGGTGGTGCACTGAAATGTAGGTCGGATAAGAAGCATAGAGTAGCATCCGACAACCGAGGCCTGATg                                                                                                                         >  2:748254/1‑150 (MQ=255)
                            cTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACATCCGATGCCTGAta                                                                                                                         >  1:1830877/1‑149 (MQ=37)
                            cTGAACTGTAGCTCGCATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATg                                                                                                                         <  1:339312/150‑1 (MQ=12)
                             tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGAATCCTGAATCGACGCTGACAAGACTTATAAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATACGACAACCGATGCCTGaata                                                                                                                        >  2:1267718/1‑147 (MQ=11)
                                aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgcc                                                                                                                     <  2:120513/150‑1 (MQ=37)
                                     aGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                                                <  2:1233740/150‑1 (MQ=37)
                                     aGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                                                >  2:876889/1‑150 (MQ=37)
                                         cGGATAAGACGGATGGCTTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg                                                                                                            <  2:1773594/150‑1 (MQ=37)
                                           taTAAGCCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc                                                                                                          <  1:798514/149‑1 (MQ=37)
                                           taTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc                                                                                                          <  2:829283/149‑1 (MQ=37)
                                           gATAAGAGGGATGGCGTCGCATCCTACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATCAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc                                                                                                          <  1:1125999/150‑1 (MQ=255)
                                           gATAAGACGGATGGCGTCGCATCCGACAACCGATTCCTCAAGCGACGCTGACGAGACTAATCAGGACTACGGGGTGAAATGAACTGTAGGTCGGATAAGACGCATAGAGTCGCATAAGACAACAGATGCCAGAGGCGCCGCTgcggcggc                                                                                                          >  2:1088532/1‑143 (MQ=11)
                                           gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc                                                                                                          <  1:2145709/150‑1 (MQ=21)
                                           gATAAGACGGATGGCGGCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc                                                                                                          <  1:499617/150‑1 (MQ=37)
                                           gATAAGACGCAGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc                                                                                                          <  1:2323114/150‑1 (MQ=17)
                                            aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACt                                                                                                         <  2:579028/150‑1 (MQ=37)
                                            aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACt                                                                                                         <  1:1800777/150‑1 (MQ=37)
                                                 aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                                                    <  1:2164782/150‑1 (MQ=37)
                                                 aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                                                    <  1:940755/150‑1 (MQ=37)
                                                     aTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc                                                                                                <  2:386942/150‑1 (MQ=34)
                                                               aTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                                                                                                                        >  2:1577926/1‑100 (MQ=33)
                                                               aTCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGc                                                                                      <  1:793387/150‑1 (MQ=11)
                                                                  cGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                                                                                   <  2:1913935/150‑1 (MQ=34)
                                                                   gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCTTATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                                                  <  2:1984415/150‑1 (MQ=21)
                                                                   gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAAGCGCCCATGAAGCGACTTATCAGGCCAACGCGGAGCACTg                                                                                  >  2:875973/1‑150 (MQ=11)
                                                                       aCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACt                                                                              >  2:1936936/1‑150 (MQ=34)
                                                                       aCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTCATGCGCCGCTCAAGCGAATAATAAGGCATAAGCGGTGCACAGAACt                                                                              >  1:599159/1‑150 (MQ=11)
                                                                          gATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                           >  1:1795349/1‑150 (MQ=34)
                                                                              ccTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTc                                                                       >  1:1932093/1‑150 (MQ=34)
                                                                                tGATGCGCCGCTGACGTTACTTATCAGGCCTACGGTGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                                     <  2:142500/150‑1 (MQ=21)
                                                                                   tGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCAGAAGCGCAGCTCAAGCGACTTATCAGGCCTACGGGGGGCACGGACCTGAAGGTCGGATa                                                                  >  2:819785/1‑150 (MQ=11)
                                                                                       ccgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACACAATCCTGATGCGCCGATCACGCGAATTATCAGGCCTACTGAGTGCCCTGAACAGTAGGACGTATAAGAc                                                              >  2:2420818/1‑150 (MQ=11)
                                                                                        cgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGTATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg                                                             <  2:1760509/150‑1 (MQ=21)
                                                                                         gcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                                                                                                                   >  1:1596126/1‑63 (MQ=35)
                                                                                            gACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc                                                                                                >  1:738095/1‑111 (MQ=32)
                                                                                            gACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc                                                                                                <  2:738095/111‑1 (MQ=32)
                                                                                             aCGTGACTTATCAGTCCTACGGGGTGCACTGAACTGTAGTTCGTATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg                                                        <  2:1814705/150‑1 (MQ=255)
                                                                                               gcgACTTATCAGGCCTACGGGTTTCACTCATCTTTAGGTCGTATAAGACGTATAGCGTCGTATCCGACATCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                      <  2:2092757/148‑1 (MQ=11)
                                                                                               gTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCg                                                      >  2:1800507/1‑150 (MQ=33)
                                                                                                  aCTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                   <  1:927238/150‑1 (MQ=17)
                                                                                                   cTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATTCCTGATGCGCCGATGACGCAACTAATCAGGCCTACTGGGTGCACTCAACTGAAGGTCGGATAAGACGCATAGCGTCga                                                  >  2:464506/1‑149 (MQ=11)
                                                                                                   cTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                   >  2:859536/1‑149 (MQ=32)
                                                                                                   cTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                   <  1:859536/149‑1 (MQ=32)
                                                                                                   cTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc                                                  >  1:886927/1‑150 (MQ=32)
                                                                                                      aTCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTCGGTCGGATAAGACGCGTAGCGTCGCATc                                               >  1:672059/1‑150 (MQ=17)
                                                                                                       tCAGGCCTACGGTTTGCTATGAACTGGAGCTCTTATAATAGGCATTGGTTCGCATCCGACAACCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                              <  2:2186908/150‑1 (MQ=11)
                                                                                                          ggCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                           >  1:1694936/1‑150 (MQ=32)
                                                                                                          ggCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                           >  1:628167/1‑150 (MQ=32)
                                                                                                            ccTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa                                         >  1:616217/1‑150 (MQ=32)
                                                                                                            ccTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa                                         <  1:1800507/150‑1 (MQ=32)
                                                                                                             cTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                                        >  1:577634/1‑150 (MQ=32)
                                                                                                                 ggggTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                    >  2:941084/1‑150 (MQ=35)
                                                                                                                    ggccACTGAACTGTAGGCCGGATAAGACGCATCGCTTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                                 <  1:1034392/147‑1 (MQ=255)
                                                                                                                      gCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                               >  2:2318612/1‑150 (MQ=37)
                                                                                                                        aCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                             <  1:208796/150‑1 (MQ=37)
                                                                                                                           gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcg                          >  1:875335/1‑150 (MQ=37)
                                                                                                                           gAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGAGACGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcc                          >  1:1252702/1‑149 (MQ=255)
                                                                                                                                  aGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                   >  2:1020166/1‑150 (MQ=37)
                                                                                                                                   ggTCGTATTAGACGCATTTCTTCGCATTCGACAACCTATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATTGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc                  <  2:1076405/150‑1 (MQ=11)
                                                                                                                                    gTCGGATGCGACGCTATGCGTCTTATCCGACAACCGATGCCTGATGCGCCGCTGCCGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcg                 >  1:2329073/1‑150 (MQ=11)
                                                                                                                                     tCGGATGCGACGCTATGCGTCTTATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAGATGCCTGATGCGCCGCTGAagc                >  2:425821/1‑150 (MQ=11)
                                                                                                                                     tCGGATGCGACGCTATGCGTCTTATCCGACAACCGATGCCTGATGCCACGCTGACCCCACTAATAAGGCCTACGCGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGATGCGCAGCTGAcga                >  2:1669728/1‑149 (MQ=11)
                                                                                                                                      ctgATAAGACACATAGCGTCGAATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGTGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg               <  2:1805367/148‑1 (MQ=255)
                                                                                                                                             gACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc        <  2:1480399/150‑1 (MQ=37)
                                                                                                                                                gCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATACGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg     >  2:1747794/1‑150 (MQ=25)
                                                                                                                                                 cATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                               <  1:569149/91‑1 (MQ=14)
                                                                                                                                                 cATCGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGc    <  1:603439/150‑1 (MQ=17)
                                                                                                                                                   tAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCt  >  2:2388712/1‑150 (MQ=33)
                                                                                                                                                    |                                                                                                                                                    
ACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCT  >  NC_000913/349621‑349917

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: