Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,760 A→G 5.0% intergenic (+188/‑252) 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,7600AG 5.0% 279.9 / 4.6 102intergenic (+188/‑252)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (40/55);  new base G (3/2);  total (43/59)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 104 total aligned reads displayed.

TGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCG  >  NC_000913/349611‑349904
                                                                                                                                                     |                                                                                                                                                
tGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa                                                                                                                                                  <  2:695324/150‑1 (MQ=255)
tGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTCTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa                                                                                                                                                  <  1:1600933/150‑1 (MQ=255)
   atCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACAGATCAGGCCTACGGGGTGAACTGAACTGTAGGTAGGATAAGa                                                                                                                                               >  1:1623650/2‑150 (MQ=255)
      gatGACGAGATTTATCATGCCTACGAGTTCCACTGTACTTTAGGTTGTAAAATAGGGATTTCTTCGAATCCGACATCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGc                                                                                                                                            <  2:659418/148‑1 (MQ=11)
        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)
                     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)
                       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)
                                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:2164782/150‑1 (MQ=37)
                                                           aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                                       <  1:940755/150‑1 (MQ=37)
                                                              gATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                                                               >  2:569149/1‑91 (MQ=38)
                                                               aTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGcc                                                                                   <  2:386942/150‑1 (MQ=34)
                                                                       gCATCCGACAACCGATGCCTGATGCGCCGCTGTCGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                                                                                                              >  1:2021837/1‑83 (MQ=18)
                                                                       gCATCCGACAACCGATGCCTGATGCGCCGCTGTCGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                                                                                                              <  2:2021837/83‑1 (MQ=18)
                                                                         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)
                                                                                                                                               tCGGATGCGACGCTATGCGTCTTATCCGACAACCGATGCCTGATGCCACGCTGACCCCACTAATAAGGCCTACGCGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGATGCGCAGCTGAcg  >  2:1669728/1‑149 (MQ=11)
                                                                                                                                                ctgATAAGACACATAGCGTCGAATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGTGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg  <  2:1805367/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCG  >  NC_000913/349611‑349904

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: