Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,767 A→T 5.1% intergenic (+195/‑245) 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,7670AT 5.1% 275.9 / 3.5 100intergenic (+195/‑245)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (37/57);  new base T (3/2);  total (41/59)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGC  >  NC_000913/349619‑349915
                                                                                                                                                    |                                                                                                                                                    
tGACGCGACTTATCTGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATAAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAt                                                                                                                                                     <  1:2031630/150‑1 (MQ=255)
  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:353523/150‑1 (MQ=37)
                        ggTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                                             <  2:1598610/150‑1 (MQ=37)
                        ggTGCACTGAACTGTAGGTCGCATAAGACGGATGGCGTCGCATCCGACAACCTATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                                             <  1:2079803/150‑1 (MQ=255)
                         gTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                                                                                                                            <  2:1482228/150‑1 (MQ=37)
                           gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                                                                                                                          >  1:103325/1‑150 (MQ=37)
                           gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTg                                                                                                                          <  1:292181/150‑1 (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:940755/150‑1 (MQ=37)
                                                   aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                                                  <  1:2164782/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                                                 <  1:859536/149‑1 (MQ=32)
                                                                                                     cTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCg                                                 >  2:859536/1‑149 (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)
                                                                                                                                                    |                                                                                                                                                    
TGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGC  >  NC_000913/349619‑349915

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: