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.

GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649861‑650054
                                                                                                    |                                                                                                 
gttgttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCATGCGAACAGCAGTTCAGttccttcc                                                                                                    >  1:2688648/1‑100 (MQ=255)
   gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggg                                                                                              <  1:1101039/100‑1 (MQ=255)
    ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggg                                                                                            <  2:1839685/100‑1 (MQ=255)
    ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggg                                                                                            <  1:2013577/100‑1 (MQ=255)
     tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggggg                                                                                           >  1:2740726/1‑99 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggagcg                                                                                       >  2:1303034/1‑95 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggagcg                                                                                       >  1:817054/1‑95 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggagcg                                                                                       >  1:3034679/1‑95 (MQ=255)
         gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggagcg                                                                                       >  1:1438394/1‑95 (MQ=255)
          tCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGt                                                                                       >  1:2675700/1‑100 (MQ=255)
           cGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGtt                                                                                      >  2:894070/1‑100 (MQ=255)
             gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCt                                                                                    >  1:1018653/1‑100 (MQ=255)
             gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCt                                                                                    <  2:1092521/100‑1 (MQ=255)
             gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCt                                                                                    >  2:922345/1‑100 (MQ=255)
              gatctaCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTg                                                                                   <  1:1475834/95‑1 (MQ=255)
              aaccggcGGTATCGAACGACGGAATACCCGCTGCGCGCGCGGCCTGCAAAATGGAACAGGCTGCGAACAGGAGTTCAGGTCCTGCCGGGGGAGCGTTCTg                                                                                   <  1:3031385/94‑1 (MQ=255)
              aaTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTg                                                                                   <  1:2211789/100‑1 (MQ=255)
               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGt                                                                                  >  2:2196004/1‑100 (MQ=255)
               aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGt                                                                                  <  1:2805950/100‑1 (MQ=255)
               aTAGACGGTATCGAACGCATGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGt                                                                                  >  2:3063176/1‑100 (MQ=255)
                tAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTa                                                                                 <  2:642875/100‑1 (MQ=255)
                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                <  2:3047061/100‑1 (MQ=255)
                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                <  2:3198467/100‑1 (MQ=255)
                 aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                <  1:1613590/100‑1 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                <  2:591215/99‑1 (MQ=255)
                  gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                               >  1:910509/1‑100 (MQ=255)
                   aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTGCAGTTCCTTCCGGGGGAGCGTTCTGTACGc                                                                              <  1:2296277/100‑1 (MQ=255)
                    cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                             <  1:1808397/100‑1 (MQ=255)
                      gTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                             <  1:2619940/98‑1 (MQ=255)
                       tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGt                                                                          <  2:231022/100‑1 (MQ=255)
                         tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTg                                                                        <  2:2385995/100‑1 (MQ=255)
                         tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTg                                                                        <  2:131307/100‑1 (MQ=255)
                          cGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgc                                                                       >  1:1491632/1‑100 (MQ=255)
                           gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcg                                                                      >  2:2528437/1‑100 (MQ=255)
                            aaCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcgc                                                                     <  2:1200376/100‑1 (MQ=255)
                             aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                    <  1:19054/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACTGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcac                                                                   >  1:1906582/1‑100 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGGGCGTTCTGTACGCAGGTTGCGcac                                                                   <  1:97784/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcac                                                                   >  1:2688937/1‑100 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcac                                                                   <  1:1349275/100‑1 (MQ=255)
                              cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATCGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcac                                                                   >  2:2442338/1‑100 (MQ=255)
                               gCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcaca                                                                  >  2:2964772/1‑100 (MQ=255)
                                ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  2:2909123/1‑100 (MQ=255)
                                 cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATa                                                                >  1:1103359/1‑100 (MQ=255)
                                 cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATa                                                                >  1:753594/1‑100 (MQ=255)
                                  tGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAg                                                               >  1:266168/1‑100 (MQ=255)
                                  tGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAg                                                               <  2:2642790/100‑1 (MQ=255)
                                   gAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGt                                                              <  1:612239/100‑1 (MQ=255)
                                       aCCCGCATAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          >  1:1831995/1‑100 (MQ=255)
                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          <  2:2688979/100‑1 (MQ=255)
                                       aCCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACACAGtcttc                                                          <  2:1825684/100‑1 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:2605451/1‑100 (MQ=255)
                                        cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:1331986/1‑100 (MQ=255)
                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        <  2:2801315/100‑1 (MQ=255)
                                         ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        <  2:1420651/100‑1 (MQ=255)
                                          cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  1:479301/1‑100 (MQ=255)
                                           gCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCa                                                      <  1:922357/100‑1 (MQ=255)
                                             agagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCAcc                                                    <  2:3222584/100‑1 (MQ=255)
                                              gagCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCg                                                   >  1:344062/1‑100 (MQ=255)
                                                gcgcgcGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAg                                                 <  1:2877902/100‑1 (MQ=255)
                                                       gCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATa                                          >  1:336070/1‑100 (MQ=255)
                                                         cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATTGTCTTCTGCACCGAGGGCGATAcc                                        >  1:348536/1‑100 (MQ=255)
                                                          tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                       >  1:1733878/1‑100 (MQ=255)
                                                           gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  1:435884/1‑100 (MQ=255)
                                                            cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGAGGGAGCGTTCTGTACGAAGGTTCCGCACATAGGCTTCTGCACCGAGGGCGATACCGAt                                     >  1:2063414/1‑100 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  2:1491654/100‑1 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  1:1158535/100‑1 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    >  1:37776/1‑100 (MQ=255)
                                                             aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  1:1273807/100‑1 (MQ=255)
                                                               aaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  >  2:1786250/1‑100 (MQ=255)
                                                                 tGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                >  1:1438156/1‑100 (MQ=255)
                                                                 tGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGTGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                <  1:2952075/100‑1 (MQ=255)
                                                                   gAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACg                               >  2:926026/1‑99 (MQ=255)
                                                                   gAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGc                              >  2:1732658/1‑100 (MQ=255)
                                                                     aCAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            <  1:1974952/100‑1 (MQ=255)
                                                                     aCAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            >  1:600304/1‑100 (MQ=255)
                                                                      cAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCg                           <  2:377127/100‑1 (MQ=255)
                                                                      cAGCGTGCGAACAACAGTTCAGGTCCTTCCGGGGGAGCGTTCTGTACGCAGCTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  <  1:2076917/93‑1 (MQ=255)
                                                                        gCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                         >  2:3175471/1‑100 (MQ=255)
                                                                          gTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAg                       >  2:110707/1‑100 (MQ=255)
                                                                           tgtgAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                      >  2:1485181/4‑100 (MQ=255)
                                                                           tGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                      >  2:2697444/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  1:2203969/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  1:337551/1‑100 (MQ=255)
                                                                            gCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     <  1:715946/100‑1 (MQ=255)
                                                                              gAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTg                   <  2:512433/100‑1 (MQ=255)
                                                                                aCAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAg                 <  2:2060975/100‑1 (MQ=255)
                                                                                aCAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAg                 <  2:1945477/100‑1 (MQ=255)
                                                                                  agcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  2:2068228/100‑1 (MQ=255)
                                                                                      gttcagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAttt           >  2:1181885/1‑100 (MQ=255)
                                                                                        tcagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTcc         <  2:3406770/100‑1 (MQ=255)
                                                                                         cagttcCTTCCGGGCGAGCGCTCTGTACGCAGGTTGCGCACATAGTCTTCTGCCCCGAGCGCGATCCCGATCACCCGCTCCGAAGCGTGAGCGATTTCCa        >  1:3279144/1‑100 (MQ=255)
                                                                                          agttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGCGCGATCCCGATCAAACGCCCGGAAGCGTCAGCGATTTCCAc       >  1:561777/1‑100 (MQ=255)
                                                                                           gttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACt      <  2:2920869/100‑1 (MQ=255)
                                                                                             tccttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgc    >  2:2434282/1‑100 (MQ=255)
                                                                                              ccttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTg     >  1:1061505/1‑98 (MQ=255)
                                                                                               cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  2:1343167/1‑100 (MQ=255)
                                                                                               cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  2:1527276/1‑100 (MQ=255)
                                                                                               cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  2:178484/1‑100 (MQ=255)
                                                                                               cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGATCGGAAGCGTGAGCGATTTCCACTgcg   >  2:1967089/1‑99 (MQ=255)
                                                                                                    |                                                                                                 
GTTGTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  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: