Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 349,717 T→C 55.0% intergenic (+145/‑295) 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,7170TC55.0% ‑2.5 / 101.1 89intergenic (+145/‑295)prpB/prpC2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (25/23);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAC  >  NC_000913/349618‑349813
                                                                                                   |                                                                                                
cTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGCCGt                                                                                                  <  2:1539063/100‑1 (MQ=255)
  gACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGTCTGATGCGCCGCTGACGTGa                                                                                                <  2:1352578/100‑1 (MQ=255)
  gACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                                                <  2:1257823/100‑1 (MQ=255)
  gACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                                                >  1:98572/1‑100 (MQ=255)
     gcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt                                                                                             >  1:1983914/1‑100 (MQ=255)
     gcgACGTATCACGCCTACAAGGTGCACTGAATTGTAGGGCAGATAATTCGGATGTCGTCGCATCCGACAACCGTTGCCTGATGCGCCGCTGACGCGAc                                                                                               <  2:127300/98‑1 (MQ=11)
      cgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACtt                                                                                             >  1:242812/1‑99 (MQ=255)
       gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                           >  2:361533/1‑100 (MQ=25)
         cTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                                         >  1:546424/1‑100 (MQ=255)
          ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAg                                                                                        <  2:546425/100‑1 (MQ=255)
           tATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAgg                                                                                       >  2:1326009/1‑100 (MQ=255)
           tATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg                                                                                       >  2:1509392/1‑100 (MQ=255)
           tATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAgg                                                                                       >  1:1346190/1‑100 (MQ=255)
           tATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAgg                                                                                       >  2:2065259/1‑100 (MQ=38)
            atcaggcCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGc                                                                                      >  2:153126‑M2/8‑100 (MQ=255)
            atcaggcCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGc                                                                                      >  1:304569‑M2/8‑100 (MQ=255)
             tcaggcCTACGAGGTGCACTGAACTGTAGATCTGATAAGACGGATGGCGTCGCATTCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGc                                                                                      <  2:2139834‑M2/93‑1 (MQ=255)
              caggcCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCt                                                                                    >  1:1897691‑M2/6‑100 (MQ=255)
                ggcCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTAc                                                                                  >  1:1950433‑M2/4‑100 (MQ=255)
                  cCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGCTGCGCCGCTGACGCGACTTATCAGGCCTACgg                                                                                >  1:373571‑M2/2‑100 (MQ=255)
                  cCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgg                                                                                >  1:1711286‑M2/2‑100 (MQ=255)
                    tACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                               >  2:1542601/1‑99 (MQ=255)
                    tACGAGGTGCACTGAACTGTAGGTCGGATAAGAAGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACgggg                                                                              <  1:38183/100‑1 (MQ=255)
                     aCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                                             >  1:1231404/1‑100 (MQ=255)
                       gAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGc                                                                           <  2:302996/100‑1 (MQ=25)
                         ggTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCAc                                                                         <  1:1571275/100‑1 (MQ=255)
                           tGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTg                                                                       >  2:1975350/1‑100 (MQ=17)
                            gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGa                                                                      >  2:1948409/1‑100 (MQ=18)
                            gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGa                                                                      >  1:762110/1‑100 (MQ=25)
                            gCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGa                                                                      >  2:383235/1‑100 (MQ=14)
                             cACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACTGGGTGCACTGaa                                                                     >  1:1954855/1‑100 (MQ=255)
                              aCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                    <  1:212846/100‑1 (MQ=17)
                              aCTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGAAAACCAATGCCTGATGAACCGATGCCAAGCATTATAAGACCTACGGAGTGAACTGaaa                                                                    >  1:1874474/1‑99 (MQ=2)
                               cTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCCCTGa                                                                      >  2:199483/1‑97 (MQ=12)
                                tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTg                                                                  <  2:592601/100‑1 (MQ=40)
                                tGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTg                                                                  <  1:420054/100‑1 (MQ=17)
                                 gAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGt                                                                 <  1:2243024/100‑1 (MQ=255)
                                 gAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGAAAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAAc                                                                    <  2:1750327/97‑1 (MQ=255)
                                   aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                >  2:205059/1‑99 (MQ=255)
                                   aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                               <  1:1717232/100‑1 (MQ=40)
                                   aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                               <  1:1150873/100‑1 (MQ=14)
                                   aCTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAg                                                               <  2:1231399/100‑1 (MQ=18)
                                   aCTGTAGGTCGAATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                <  2:1615086/99‑1 (MQ=9)
                                    cTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTa                                                                <  2:1723981/98‑1 (MQ=17)
                                    cTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAgg                                                              <  1:199213/100‑1 (MQ=14)
                                    cTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAgg                                                              <  2:304570/100‑1 (MQ=17)
                                     tGTAGGTCGGATAAGACGTATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAgg                                                              <  1:588247/99‑1 (MQ=38)
                                         ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCg                                                           <  2:2013995/98‑1 (MQ=14)
                                         ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                         <  2:465787/100‑1 (MQ=14)
                                         ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                         <  1:2051589/100‑1 (MQ=18)
                                         ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGa                                                         >  1:760702/1‑100 (MQ=18)
                                          gTCGAATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACTCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                        <  2:242812/100‑1 (MQ=255)
                                           tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCgg                                                          <  1:60696/97‑1 (MQ=14)
                                           tCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATa                                                       <  1:814596/100‑1 (MQ=17)
                                            cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATaa                                                      <  1:1271378/100‑1 (MQ=18)
                                            cGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATaa                                                      >  1:460482/1‑100 (MQ=18)
                                              gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                    >  1:2055951/1‑100 (MQ=255)
                                              gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGa                                                    >  1:1107907/1‑100 (MQ=18)
                                               aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAc                                                   >  1:754797/1‑100 (MQ=255)
                                               aTAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                     <  2:1636508/98‑1 (MQ=14)
                                                tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCTGATAAGACg                                                  >  2:1063574/1‑100 (MQ=255)
                                                tAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACg                                                  <  2:1008283/100‑1 (MQ=255)
                                                 gatctggATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGc                                                 <  1:834768/95‑1 (MQ=2)
                                                 aaGACGGATGGCGTCGCATCCGACAACCGACGCCGAACGCGAGACCGACGTGACGTAAAAGACCTAGGCGCTGAATAGAACTCGAGGGCGGATataac                                                   >  2:1723662/1‑94 (MQ=255)
                                                  aGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                                                >  2:1267206/1‑100 (MQ=18)
                                                   gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAt                                               >  2:444840/1‑100 (MQ=17)
                                                   gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCAt                                               >  1:1537182/1‑100 (MQ=17)
                                                   gACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAGAAGACg                                                  >  2:1558075/1‑97 (MQ=12)
                                                     cGGATGGCGTCGCATCCGACAACCGATGCCTGTTGCGCCGCTCACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAg                                             <  2:1574491/100‑1 (MQ=255)
                                                      ggATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCagatc                                            >  2:834768/1‑96 (MQ=2)
                                                            cGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                                      >  1:1146575/1‑100 (MQ=33)
                                                                  aTCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGAc                                >  2:1292359/1‑100 (MQ=33)
                                                                      gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACa                               >  2:973612/1‑97 (MQ=32)
                                                                      gACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcc                            <  1:370782/100‑1 (MQ=32)
                                                                          aCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAcaaa                          >  2:1760357/1‑95 (MQ=18)
                                                                           ccGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATg                        <  2:2208558/99‑1 (MQ=32)
                                                                            cGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGcc                      <  1:119608/100‑1 (MQ=32)
                                                                               tGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                   <  1:205059/100‑1 (MQ=33)
                                                                               tGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGa                   <  1:153126/100‑1 (MQ=33)
                                                                                gCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                  <  2:1388303/100‑1 (MQ=32)
                                                                                cgatcttGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAt                  <  1:1387570/94‑1 (MQ=16)
                                                                                 ccTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCt                     >  1:987741/1‑96 (MQ=33)
                                                                                  cTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATgc                >  2:1584594/1‑100 (MQ=33)
                                                                                      tGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAtagatcg            >  2:1387564/1‑94 (MQ=16)
                                                                                      tGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg            <  1:840558/100‑1 (MQ=33)
                                                                                      tGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccg            >  2:1941438/1‑100 (MQ=33)
                                                                                       gcgcCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGcgccgc           <  1:1412430/100‑1 (MQ=33)
                                                                                            gcTGACGTGACTTATCGGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcg      >  1:835324/1‑100 (MQ=18)
                                                                                                aCGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCTTCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc  <  1:1584605/100‑1 (MQ=17)
                                                                                                aCGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAc  <  1:2002688/100‑1 (MQ=32)
                                                                                                   |                                                                                                
CTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGAC  >  NC_000913/349618‑349813

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