Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 649,964:1 +G 100% 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.G94.4% 299.4 / 9.0 89coding (531/909 nt)citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base G (35/49);  total (39/50)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑G‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649864‑650054
                                                                                                  |                                                                                                
gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGGGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg                                                                                              <  2:174291/100‑1 (MQ=255)
gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg                                                                                              <  1:2355907/100‑1 (MQ=255)
gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg                                                                                              <  1:516094/100‑1 (MQ=255)
 ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg‑g                                                                                            <  1:2436689/100‑1 (MQ=255)
  tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg‑g                                                                                            >  1:489644/1‑99 (MQ=255)
   aGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg‑gga                                                                                          <  2:947511/100‑3 (MQ=255)
    gCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg‑ggag                                                                                         >  2:829534/1‑97 (MQ=255)
      gacGGAATAGACGATATCGAACGCCTGAATAACCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg‑ggagcg                                                                                       >  1:1730911/3‑95 (MQ=255)
      gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑g‑gg‑ggagcg                                                                                       >  1:1658704/1‑95 (MQ=255)
       tCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGt                                                                                       >  1:730768/1‑100 (MQ=255)
        cGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGtt                                                                                      >  1:1777562/1‑100 (MQ=255)
        cGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGtt                                                                                      <  1:1367043/100‑1 (MQ=255)
         ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTc                                                                                     <  1:1494883/100‑1 (MQ=255)
         ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTc                                                                                     <  1:512050/100‑1 (MQ=255)
         ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTc                                                                                     <  1:760580/100‑1 (MQ=255)
          gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCt                                                                                    <  1:2479997/100‑1 (MQ=255)
          gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCt                                                                                    <  1:2505605/100‑1 (MQ=255)
           aaTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCTGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTg                                                                                   <  1:2366850/100‑1 (MQ=255)
            aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCTCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGt                                                                                  <  1:1889287/100‑1 (MQ=255)
            aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGt                                                                                  <  1:2148470/100‑1 (MQ=255)
              aGACGGTATCGAACGTCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTAc                                                                                <  1:2517562/100‑1 (MQ=255)
              aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTAc                                                                                >  2:102700/1‑100 (MQ=255)
              aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTAc                                                                                <  2:818454/100‑1 (MQ=255)
               gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACg                                                                               >  1:49130/1‑100 (MQ=255)
                 cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCa                                                                             >  2:1971738/1‑100 (MQ=255)
                 cGGTATCCAACACCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCa                                                                             >  2:1341259/1‑100 (MQ=255)
                  ggTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAg                                                                            <  2:780983/100‑1 (MQ=255)
                    tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGt                                                                          <  1:103472/100‑1 (MQ=255)
                    tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGt                                                                          <  1:103449/100‑1 (MQ=255)
                      tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTg                                                                        >  2:13845/1‑100 (MQ=255)
                       cGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTgc                                                                       <  1:1664954/100‑1 (MQ=255)
                        gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTgcg                                                                      >  1:489221/1‑100 (MQ=255)
                        gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTgcg                                                                      >  1:1254954/1‑100 (MQ=255)
                          aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGca                                                                    <  1:231514/100‑1 (MQ=255)
                           cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGcac                                                                   <  2:2502547/100‑1 (MQ=255)
                            gCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTCGCGcaca                                                                  <  1:1341268/100‑1 (MQ=255)
                             ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  1:1303964/1‑100 (MQ=255)
                             ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 <  1:1214728/100‑1 (MQ=255)
                                 aaTACCCGCAGAGCGCGCGGCCAGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtc                                                             >  1:1076243/1‑100 (MQ=255)
                                   tACCCGCAAAGCGCGCAGCCCGCAAAGTGGAACAGCGTGCGAACAGCAGTTCAGTTCCGTC‑G‑G‑GGAGCACGTACTGTACGCAGGTTGCGCACATat                                                               >  2:262679/1‑95 (MQ=255)
                                     cccGCAGAGGGCGCGACCTGCAAAATCAAACAGCGTGCGAACAGCAGTTCAGCTCCTTC‑C‑GACG‑GAGCGTTCTGTACGCAGGGCGCGCACAtcggct                                                            >  2:738069/1‑92 (MQ=255)
                                     cccGCAGAGCGCGCTGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          <  1:2034775/99‑1 (MQ=255)
                                     cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCCGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTGCGCAGGTTGCGCACCTAGtctcc                                                          >  2:1097785/1‑97 (MQ=255)
                                     cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         <  1:242378/100‑1 (MQ=255)
                                     cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:96704/1‑100 (MQ=255)
                                      ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        <  2:1933391/100‑1 (MQ=255)
                                      ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        >  1:473993/1‑100 (MQ=255)
                                      ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        >  1:1854386/1‑100 (MQ=255)
                                       cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  1:1267237/1‑100 (MQ=255)
                                       cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  1:2249303/1‑100 (MQ=255)
                                             gcgcgcGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCGG‑G‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAg                                                 >  2:508203/1‑100 (MQ=255)
                                                cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGTAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  1:2428657/100‑1 (MQ=255)
                                                cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  1:1868377/100‑1 (MQ=255)
                                                cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  1:1381703/100‑1 (MQ=255)
                                                 gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             <  1:13846/100‑1 (MQ=255)
                                                 gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             <  2:46324/100‑1 (MQ=255)
                                                  cgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCTGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGa                                            >  2:2500742/1‑100 (MQ=255)
                                                   ggCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGAt                                           >  1:110164/1‑100 (MQ=255)
                                                    gCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATa                                          >  2:305726/1‑100 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGTGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        >  1:1866789/1‑100 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTGCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        >  2:1368355/1‑100 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:2103583/100‑1 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:528537/100‑1 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:1952443/100‑1 (MQ=255)
                                                       tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                       >  1:1166513/1‑100 (MQ=255)
                                                        gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  1:1515451/1‑100 (MQ=255)
                                                        gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  2:2346625/1‑100 (MQ=255)
                                                         cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGAt                                     <  1:1335287/100‑1 (MQ=255)
                                                          aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    >  1:390524/1‑100 (MQ=255)
                                                          aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  1:2346640/100‑1 (MQ=255)
                                                           aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  1:1352103/1‑100 (MQ=255)
                                                           aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  1:1691754/1‑100 (MQ=255)
                                                            aaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  <  2:2042529/100‑1 (MQ=255)
                                                            aaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  <  1:1695095/100‑1 (MQ=255)
                                                             aTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑G‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                <  1:1161541/100‑1 (MQ=255)
                                                                  aCAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            <  1:2528038/100‑1 (MQ=255)
                                                                  aCAGCGTGCGAACAGCAGGTCAGTGCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCAGATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            <  1:714924/100‑1 (MQ=255)
                                                                   cAGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCg                           <  2:946451/100‑1 (MQ=255)
                                                                    aGCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCgg                          <  2:2556262/100‑1 (MQ=255)
                                                                     gCGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                         >  1:1411928/1‑100 (MQ=255)
                                                                      cGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGaa                        <  2:2009376/100‑1 (MQ=255)
                                                                      cGTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGaa                        <  2:290409/100‑1 (MQ=255)
                                                                       gTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGCTACCGATCAAACGCTCGGAAg                       >  2:575293/1‑100 (MQ=255)
                                                                       gTGCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAg                       <  2:1738032/100‑1 (MQ=255)
                                                                         gCGAACAGTAGTTCAGTTCCTTC‑CGG‑GG‑GGGCTTTCTTCACGCAGGTTGCCCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCg                           <  2:2062599/94‑1 (MQ=255)
                                                                         gCGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     <  1:1696562/100‑1 (MQ=255)
                                                                          cGAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  2:79467/1‑99 (MQ=255)
                                                                           gAACAGCAGTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTg                   <  1:1482940/100‑1 (MQ=255)
                                                                              cagcagTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATCCCGATCACACGCTCGGCAGCGTAAGc                >  1:2136492/1‑100 (MQ=255)
                                                                              cagcagTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGc                <  2:865412/100‑1 (MQ=255)
                                                                              cagcagTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGc                <  1:2220127/100‑1 (MQ=255)
                                                                              cagcagTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGc                >  1:239237/1‑100 (MQ=255)
                                                                               agcagTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  1:540379/100‑1 (MQ=255)
                                                                                gcagTTCAGTTCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  2:1166220/99‑1 (MQ=255)
                                                                                gcagTTCAGTTCCTTC‑C‑G‑GGCGAGCGTTCTGTACGCAGGTTGCGCACCTAGTCTTCTGCACCGGGGGCGATACCGATCAAACGCTCGGCAGCGTGAGCGa              >  2:591450/1‑100 (MQ=255)
                                                                                 cagttcagttCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAt             >  2:1885241/1‑100 (MQ=255)
                                                                                 cagttcagttCCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAt             <  2:2542081/100‑1 (MQ=255)
                                                                                      cagttcCTTC‑CGG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa        >  2:316664/1‑100 (MQ=255)
                                                                                            cttc‑cgG‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  1:2424322/1‑100 (MQ=255)
                                                                                                  |                                                                                                
GTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTC‑C‑G‑GG‑GAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649864‑650054

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: ‑