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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913349,7400CA 5.0% 318.5 / 4.3 135intergenic (+168/‑272)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base C (52/72);  new base A (5/2);  total (60/75)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

TTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTG  >  NC_000913/349592‑349888
                                                                                                                                                    |                                                                                                                                                    
ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGACTTAGTCGCCTCTGACGTGACTTATCAGGCCTACGGGGGGc                                                                                                                                                      >  1:1689834/1‑149 (MQ=17)
  cGACAGCCGTTGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCAc                                                                                                                                                    <  1:891117/149‑1 (MQ=255)
        ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGCCGTCCAATCCGACAACCGATGCCTGCTTCCACGCTGACACAACGAATCAGGCCGACCGGGGGCGCGGAAca                                                                                                                                              >  2:646950/1‑148 (MQ=11)
          gaggCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                                                                                            <  1:24045/146‑1 (MQ=255)
          gATGCCTGAGGCGCCGCTGACGCGACTTATCATGCCTACGAGGTGCACTAAACTGTAGGTCGGATAAGACGGATGGCGCCGCATCCGACAACCGATGCCGTATGCGCCGCTGACGTGACTTCTCAGGCCTACGGGGTGCACTGAACTGt                                                                                                                                            <  2:1421835/149‑1 (MQ=14)
          gATGCCGTATGCGTCGCTGACGCGACTTATCAGGCCTGCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATATCGCGTTCGGGATGAAATGAACTGt                                                                                                                                            <  2:400094/149‑1 (MQ=14)
            tGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                                                                                                          <  2:1339486/149‑1 (MQ=255)
              ccTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                                                                                                                        >  1:889483/1‑149 (MQ=255)
                tGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGGCATGAACTGTAGGTCg                                                                                                                                      >  2:1261281/1‑149 (MQ=255)
                tGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                                                                                                      <  1:969410/149‑1 (MQ=255)
                tGATGCGTCGCTGACGCGACTTAGCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                                                                                                      <  1:682805/149‑1 (MQ=255)
                tGATGCGTCGCTGACGCGACTTAACAGGCCTACGAGGTCCACTGAAGTGTAGGTCGGATAAGACGAATGGCGGTGCATCCGACAACCGATGACTGATGCGCCGCTGACGTGGCTTATAAGCCAAAGGTGTTGCACTGAGCTATGGGtgg                                                                                                                                      >  2:1355568/1‑147 (MQ=255)
                agATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGCCGGATAAGCCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                                                                                                      <  1:705460/148‑1 (MQ=255)
                      gTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                                                                                                <  1:1405665/149‑1 (MQ=255)
                      gTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACCAAGCAGGCCCACGGGGGGCACGGAACCGTAGGGGGGGTAAg                                                                                                                                >  1:126270/1‑149 (MQ=14)
                       tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGa                                                                                                                               >  1:1730241/1‑149 (MQ=255)
                       tCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                                                                                               >  1:303590/1‑149 (MQ=38)
                         gCTGGCGCGACTTATCAGGCCTACGAGGGGCCCTGCACAGTAGGTCGGATGAGACGGATGGCGTCGCATCAGACAACCGAGGCCGGATTCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg                                                                                                                             <  2:1645217/149‑1 (MQ=14)
                         gCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACg                                                                                                                             <  2:537838/149‑1 (MQ=38)
                         gCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACg                                                                                                                             >  2:1801309/1‑149 (MQ=17)
                          cTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGc                                                                                                                            >  2:1622349/1‑149 (MQ=17)
                             aCGCGACTGATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATa                                                                                                                         <  2:1812864/149‑1 (MQ=255)
                                 gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCATGAACTGTCGGTCGGATAGAACGCATAGCGt                                                                                                                     >  2:1133595/1‑149 (MQ=11)
                                  tctgaTCAGGCCTACGAGGGGCACTGAACGGTATGTCGGATAAGACGGATGGCGTCGCATCCGACAACCTATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTTCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    <  2:1463787/145‑1 (MQ=14)
                                  aCTTATCGGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    <  2:1665353/149‑1 (MQ=255)
                                  aCTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCAACGGGGGGCACTGAACGGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    >  2:474770/1‑149 (MQ=255)
                                    ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATCGCGTccc                                                                                                                  >  2:128174/1‑147 (MQ=255)
                                        cAGGCCTACGAGTTGCACTGGACCGTAGGGCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                                                                                              <  1:490179/149‑1 (MQ=255)
                                         aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGAGTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTTGGATAAGACGCATAGCGTCGCATCCg                                                                                                             >  1:10695/1‑149 (MQ=255)
                                         aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                                                                                             >  2:850035/1‑149 (MQ=255)
                                         aGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTAGCAGGCCTACGGGGTGCACCGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCg                                                                                                             >  1:70706/1‑149 (MQ=255)
                                          ggCCTACGAGGTGCCCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATTCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGa                                                                                                            <  2:942326/149‑1 (MQ=255)
                                              tACGAGGTGCACTTAACTGGCGGTCGAGTAAGCCGCTGTGGGTCGCCTCCGACAGCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAc                                                                                                        <  1:755505/149‑1 (MQ=11)
                                              tACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGGGCACTGAACTGTAGGCCGGCTAACCAGAATAGCGTCGCATCCGACAAc                                                                                                        >  1:328396/1‑149 (MQ=14)
                                               aCGAGGTGCACTAAACTGTAGGGCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                                                                                                       <  1:642418/149‑1 (MQ=255)
                                                 gAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATACGACGCATAGCGTCGCATCCGGCAAccca                                                                                                     >  2:1418001/1‑147 (MQ=255)
                                                  gggggCACTGACCTGGAGGTCGGATAAGACGGATGGCGTCGGATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAt                                                                                                    <  2:66709/145‑1 (MQ=255)
                                                    gTGCACTGAACTGTAGGGCGGCTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGc                                                                                                  <  1:459034/149‑1 (MQ=255)
                                                     tGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                                                                                                 >  1:1203479/1‑149 (MQ=39)
                                                      gCACTGACCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                                                                                                <  1:974344/149‑1 (MQ=38)
                                                           gAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgc                                                                                           <  2:1714607/149‑1 (MQ=38)
                                                             gcTGGAGGTCGGATAAGACGCATGGCGTCGTATTCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgcgc                                                                                         <  2:1573239/148‑1 (MQ=255)
                                                                gTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATACGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc                                                                                      >  1:1669712/1‑149 (MQ=38)
                                                                gTAGGTCGGATAAGACGGATGGCAACGGATGTGTCAACCGATGCCTGATGCGCCGCGGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAAACGCATAGCGTCGCATCCGACAACCGATGCCTGATGgcccgc                                                                                      >  2:1749941/1‑149 (MQ=11)
                                                                   ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                   <  1:459851/149‑1 (MQ=37)
                                                                   ggTCGGATAAGACGGAGGGCGTCGCATCCGCCAACCGATGCCTGAGGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                   <  1:82189/149‑1 (MQ=11)
                                                                    gTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGTACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAc                                                                                  <  2:755565/149‑1 (MQ=38)
                                                                     tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCAATGCCTGATGCCCCGCGGAcg                                                                                 >  2:680347/1‑149 (MQ=255)
                                                                      cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgc                                                                                >  1:1051652/1‑149 (MQ=39)
                                                                        ttacaGACGGGTGGCGCTGCATCCGACCACCGATGGCTGATGCGCCGCTGACGTGACTTATGAGGCCTACGGTGGGCACTGAACTGTAGGTGGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                              <  2:543704/145‑1 (MQ=11)
                                                                        gATACGACGCATGGGGTGGCATCCGACAACCGATGCCGGATGCGCCGCTGACGTGACATATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGTATAGCGTTGCATCCGACAACCGCTGCCTGATGCGCCGCTGACGCGa                                                                              <  2:1355298/149‑1 (MQ=11)
                                                                        aaTAAGCCGGAGGGCGTCGCAGCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                              <  2:117128/148‑1 (MQ=255)
                                                                         aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAAGCCTGATGCGCCGCTGACGCGAc                                                                             <  1:1001703/149‑1 (MQ=37)
                                                                           aaGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt                                                                           <  1:1204792/149‑1 (MQ=38)
                                                                             gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                         <  2:889483/149‑1 (MQ=37)
                                                                             gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                         <  1:1610724/149‑1 (MQ=37)
                                                                             gACGGATGGCGCCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACACCCGCTGCCTGATGCGCCGCTGCCGCGACtcac                                                                         >  2:249869/1‑146 (MQ=11)
                                                                              aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCCATGCCTGACGCGCCGCTGACACGACttcta                                                                        >  2:1159436/1‑146 (MQ=255)
                                                                              aCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCACAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATc                                                                        >  1:415212/1‑149 (MQ=38)
                                                                                  aTGTCGGCGCATCCGCTAAACGAGACCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGc                                                                    <  2:372529/149‑1 (MQ=11)
                                                                                  aTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGc                                                                    >  1:94755/1‑149 (MQ=35)
                                                                                  aTGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGc                                                                    <  1:1599702/149‑1 (MQ=35)
                                                                                     gcggcgCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTa                                                                 <  1:171955/149‑1 (MQ=17)
                                                                                       gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGCCTTCTCAGGCCTCCg                                                               >  2:870768/1‑149 (MQ=12)
                                                                                       gTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACg                                                               <  1:7834/149‑1 (MQ=32)
                                                                                        tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGAAGCTGACGCCGCCTCTCAGGCCTACgg                                                              >  2:1405272/1‑149 (MQ=11)
                                                                                        tCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGGAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGAGGCCAGCTGCGCCGCTCACGCGCCTACGCAGGccaccgg                                                              >  2:1777132/1‑144 (MQ=11)
                                                                                           cctCCGACAACCGATGCCTGATGCGCCGCCGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAGGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                           <  2:1777512/147‑1 (MQ=12)
                                                                                            aTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTCATGACCAACGGACGCGACGTATCAGGCCAACGGGGTg                                                          >  2:299411/1‑149 (MQ=11)
                                                                                                gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACt                                                      >  1:586642/1‑149 (MQ=35)
                                                                                                 aCAACCGATGCCTGGTGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCCCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                     <  2:1785970/149‑1 (MQ=25)
                                                                                                 aCAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCTCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                     <  1:913418/149‑1 (MQ=25)
                                                                                                 aCAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCAGCTTCCCAGCGGAAGCAACGTAGCGGGCCTCCGGGGCGCACta                                                     >  1:1290276/1‑148 (MQ=11)
                                                                                                     gcgctGCCTGATGCGCTGCTGACGGGACTTATCAGGCCCACGGGGTGCACTGACCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGCACt                                                 <  2:1379268/145‑1 (MQ=11)
                                                                                                      cGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                                                                                                                    >  2:460494/1‑81 (MQ=35)
                                                                                                        aTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGGCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCGTGAACTGTa                                              >  1:20578/1‑149 (MQ=255)
                                                                                                         tGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCACAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                             <  1:1115486/149‑1 (MQ=25)
                                                                                                           ccTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACg                                                               <  1:626785/129‑1 (MQ=32)
                                                                                                              gATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                        >  2:579159/1‑149 (MQ=35)
                                                                                                                tGCCCCGCCGACGGGACTTATCAGGCCCACGGGGGGCACTGGACTGGAGGTAGGATAAGACGCAGAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                      <  2:592460/149‑1 (MQ=11)
                                                                                                                  cgccgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCTGATaa                                    <  2:1203479/149‑1 (MQ=25)
                                                                                                                   gccgcTGACGTGACTTAGCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                   <  1:350641/149‑1 (MQ=25)
                                                                                                                    ccgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGCGCACTGAACTGTAGGTCGGATAAGa                                  <  1:1014176/149‑1 (MQ=25)
                                                                                                                    ccgcTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                  <  2:1761111/149‑1 (MQ=17)
                                                                                                                    ccgcAGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGCGCACTGAACTGTAGGTCGGATAAGa                                  <  1:1014177/149‑1 (MQ=12)
                                                                                                                     cgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGCCTTATCAGGCAGCCGGGGTGGCCTGAACTGTAGACCGGATAAGAc                                 >  2:485059/1‑149 (MQ=11)
                                                                                                                     cgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGACGCGGACGCGACTTCTAAGGCCTACGGGGTGCACTGAACTGAAGGCCGGATAAGAc                                 >  2:29246/1‑149 (MQ=11)
                                                                                                                            gTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATTCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGc                          >  1:195432/1‑149 (MQ=18)
                                                                                                                            gTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGc                          <  1:1358235/149‑1 (MQ=33)
                                                                                                                             tGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGCGCACTGAACTGCAGGTCGGATAAGACGCATAGCg                         >  2:1346139/1‑149 (MQ=17)
                                                                                                                              gACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGt                        >  1:406649/1‑149 (MQ=32)
                                                                                                                               aattATCAGGCCTCGGGGGTGCACTGAACTGTAGGCCGGATAAGCCGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                       <  2:206057/147‑1 (MQ=11)
                                                                                                                               aCTGATCGGGCCTCCGGGGGGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTc                       <  2:1512900/149‑1 (MQ=255)
                                                                                                                                 ttATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGc                     <  1:1607776/149‑1 (MQ=32)
                                                                                                                                   aTCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCCTAGCGTCGCAt                   >  2:1570328/1‑149 (MQ=17)
                                                                                                                                        gCCTACGGGGGGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc              <  2:554265/149‑1 (MQ=17)
                                                                                                                                          cTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa            >  1:1051202/1‑149 (MQ=32)
                                                                                                                                           tACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGACCTACGGGGTGCACTGAACTGTAGGCCGGATAAGACGCATAGCGTCGCATCCGACAAc           >  2:1309806/1‑149 (MQ=17)
                                                                                                                                             cGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCg         >  1:1765499/1‑149 (MQ=33)
                                                                                                                                                  tGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCCATGcc    >  1:1404501/1‑149 (MQ=25)
                                                                                                                                                    |                                                                                                                                                    
TTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTG  >  NC_000913/349592‑349888

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