Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 649,964 +G coding (527/909 nt) citE ← citrate lyase, citryl‑ACP lyase (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913649,9601.G97.6% 435.2 / 2.2 125coding (531/909 nt)citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (66/56);  total (68/57)
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

Only 100 of 152 total aligned reads displayed.

GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑G‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649861‑650054
                                                                                                    |                                                                                                 
gttgttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCATGCGAACAGCAGTTCAGttccttcc                                                                                                    >  1:2688648/1‑100 (MQ=255)
   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g                                                                                              <  1:1101039/100‑1 (MQ=255)
    ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g‑g                                                                                            <  2:1839685/100‑1 (MQ=255)
    ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g‑g                                                                                            <  1:2013577/100‑1 (MQ=255)
     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g‑gg                                                                                           >  1:2740726/1‑99 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g‑ggagcg                                                                                       >  2:1303034/1‑95 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g‑ggagcg                                                                                       >  1:817054/1‑95 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g‑ggagcg                                                                                       >  1:3034679/1‑95 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑g‑g‑g‑ggagcg                                                                                       >  1:1438394/1‑95 (MQ=255)
          tCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGt                                                                                       >  1:2675700/1‑100 (MQ=255)
           cGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGtt                                                                                      >  2:894070/1‑100 (MQ=255)
             gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCt                                                                                    >  1:1018653/1‑100 (MQ=255)
             gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCt                                                                                    <  2:1092521/100‑1 (MQ=255)
             gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCt                                                                                    >  2:922345/1‑100 (MQ=255)
              gatctaCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTg                                                                                   <  1:1475834/95‑1 (MQ=255)
              aaccggcGGTATCGAACGACGGAATACCCGCTGCGCGCGCGGCCTGCAAAATGGAACAGGCTGCGAACAGGAGTTCAGGTCCTGCCGG‑G‑G‑GAGCGTTCTg                                                                                   <  1:3031385/94‑1 (MQ=255)
              aaTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTg                                                                                   <  1:2211789/100‑1 (MQ=255)
               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGt                                                                                  >  2:2196004/1‑100 (MQ=255)
               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGt                                                                                  <  1:2805950/100‑1 (MQ=255)
               aTAGACGGTATCGAACGCATGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGt                                                                                  >  2:3063176/1‑100 (MQ=255)
                tAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTa                                                                                 <  2:642875/100‑1 (MQ=255)
                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTAc                                                                                <  2:3047061/100‑1 (MQ=255)
                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTAc                                                                                <  2:3198467/100‑1 (MQ=255)
                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTAc                                                                                <  1:1613590/100‑1 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTAc                                                                                <  2:591215/99‑1 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACg                                                                               >  1:910509/1‑100 (MQ=255)
                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTGCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGc                                                                              <  1:2296277/100‑1 (MQ=255)
                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCa                                                                             <  1:1808397/100‑1 (MQ=255)
                      gTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCa                                                                             <  1:2619940/98‑1 (MQ=255)
                       tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGt                                                                          <  2:231022/100‑1 (MQ=255)
                         tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTg                                                                        <  2:2385995/100‑1 (MQ=255)
                         tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTg                                                                        <  2:131307/100‑1 (MQ=255)
                          cGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTgc                                                                       >  1:1491632/1‑100 (MQ=255)
                           gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTgcg                                                                      >  2:2528437/1‑100 (MQ=255)
                            aaCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTgcgc                                                                     <  2:1200376/100‑1 (MQ=255)
                             aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGca                                                                    <  1:19054/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACTGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGcac                                                                   >  1:1906582/1‑100 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GGGCGTTCTGTACGCAGGTTGCGcac                                                                   <  1:97784/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGcac                                                                   >  1:2688937/1‑100 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGcac                                                                   <  1:1349275/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATCGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGcac                                                                   >  2:2442338/1‑100 (MQ=255)
                               gCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGcaca                                                                  >  2:2964772/1‑100 (MQ=255)
                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  2:2909123/1‑100 (MQ=255)
                                 cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATa                                                                >  1:1103359/1‑100 (MQ=255)
                                 cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATa                                                                >  1:753594/1‑100 (MQ=255)
                                  tGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAg                                                               >  1:266168/1‑100 (MQ=255)
                                  tGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAg                                                               <  2:2642790/100‑1 (MQ=255)
                                   gAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGt                                                              <  1:612239/100‑1 (MQ=255)
                                       aCCCGCATAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  1:1831995/1‑100 (MQ=255)
                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          <  2:2688979/100‑1 (MQ=255)
                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACACAGtcttc                                                          <  2:1825684/100‑1 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:2605451/1‑100 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:1331986/1‑100 (MQ=255)
                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        <  2:2801315/100‑1 (MQ=255)
                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        <  2:1420651/100‑1 (MQ=255)
                                          cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  1:479301/1‑100 (MQ=255)
                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCa                                                      <  1:922357/100‑1 (MQ=255)
                                             agagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCAcc                                                    <  2:3222584/100‑1 (MQ=255)
                                              gagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCg                                                   >  1:344062/1‑100 (MQ=255)
                                                gcgcgcGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAg                                                 <  1:2877902/100‑1 (MQ=255)
                                                       gCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATa                                          >  1:336070/1‑100 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATTGTCTTCTGCACCGAGGGCGATAcc                                        >  1:348536/1‑100 (MQ=255)
                                                          tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                       >  1:1733878/1‑100 (MQ=255)
                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  1:435884/1‑100 (MQ=255)
                                                            cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑GAG‑G‑GAGCGTTCTGTACGAAGGTTCCGCACATAGGCTTCTGCACCGAGGGCGATACCGAt                                     >  1:2063414/1‑100 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  2:1491654/100‑1 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  1:1158535/100‑1 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    >  1:37776/1‑100 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  1:1273807/100‑1 (MQ=255)
                                                               aaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  >  2:1786250/1‑100 (MQ=255)
                                                                 tGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                >  1:1438156/1‑100 (MQ=255)
                                                                 tGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑G‑GTG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                <  1:2952075/100‑1 (MQ=255)
                                                                   gAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACg                               >  2:926026/1‑99 (MQ=255)
                                                                   gAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGc                              >  2:1732658/1‑100 (MQ=255)
                                                                     aCAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            <  1:1974952/100‑1 (MQ=255)
                                                                     aCAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            >  1:600304/1‑100 (MQ=255)
                                                                      cAGCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCg                           <  2:377127/100‑1 (MQ=255)
                                                                      cAGCGTGCGAACAACAGTTCAGGTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGCTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  <  1:2076917/93‑1 (MQ=255)
                                                                        gCGTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                         >  2:3175471/1‑100 (MQ=255)
                                                                          gTGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAg                       >  2:110707/1‑100 (MQ=255)
                                                                           tgtgAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                      >  2:1485181/4‑100 (MQ=255)
                                                                           tGCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                      >  2:2697444/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  1:2203969/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  1:337551/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     <  1:715946/100‑1 (MQ=255)
                                                                              gAACAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTg                   <  2:512433/100‑1 (MQ=255)
                                                                                aCAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAg                 <  2:2060975/100‑1 (MQ=255)
                                                                                aCAGCAGTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAg                 <  2:1945477/100‑1 (MQ=255)
                                                                                  agcagTTCAGTTCCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  2:2068228/100‑1 (MQ=255)
                                                                                      gttcagttcCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAttt           >  2:1181885/1‑100 (MQ=255)
                                                                                        tcagttcCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTcc         <  2:3406770/100‑1 (MQ=255)
                                                                                         cagttcCTTCC‑G‑G‑GCGAGCGCTCTGTACGCAGGTTGCGCACATAGTCTTCTGCCCCGAGCGCGATCCCGATCACCCGCTCCGAAGCGTGAGCGATTTCCa        >  1:3279144/1‑100 (MQ=255)
                                                                                          agttcCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGCGCGATCCCGATCAAACGCCCGGAAGCGTCAGCGATTTCCAc       >  1:561777/1‑100 (MQ=255)
                                                                                           gttcCTTCCGG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACt      <  2:2920869/100‑1 (MQ=255)
                                                                                             tccttccgG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgc    >  2:2434282/1‑100 (MQ=255)
                                                                                              ccttccgG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTg     >  1:1061505/1‑98 (MQ=255)
                                                                                               cttccgG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  2:1343167/1‑100 (MQ=255)
                                                                                               cttccgG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  2:1527276/1‑100 (MQ=255)
                                                                                               cttccgG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  2:178484/1‑100 (MQ=255)
                                                                                               cttccgG‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGATCGGAAGCGTGAGCGATTTCCACTgcg   >  2:1967089/1‑99 (MQ=255)
                                                                                                    |                                                                                                 
GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCC‑G‑G‑G‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649861‑650054

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