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

GTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  NC_000913/649864‑650054
                                                                                                  |                                                                                                
gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGGGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggg                                                                                              <  2:174291/100‑1 (MQ=255)
gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggg                                                                                              <  1:2355907/100‑1 (MQ=255)
gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggg                                                                                              <  1:516094/100‑1 (MQ=255)
 ttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggg                                                                                            <  1:2436689/100‑1 (MQ=255)
  tAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggg                                                                                            >  1:489644/1‑99 (MQ=255)
   aGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggggga                                                                                          <  2:947511/100‑3 (MQ=255)
    gCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggag                                                                                         >  2:829534/1‑97 (MQ=255)
      gacGGAATAGACGATATCGAACGCCTGAATAACCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggagcg                                                                                       >  1:1730911/3‑95 (MQ=255)
      gTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggagcg                                                                                       >  1:1658704/1‑95 (MQ=255)
       tCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGt                                                                                       >  1:730768/1‑100 (MQ=255)
        cGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGtt                                                                                      >  1:1777562/1‑100 (MQ=255)
        cGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGtt                                                                                      <  1:1367043/100‑1 (MQ=255)
         ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTc                                                                                     <  1:1494883/100‑1 (MQ=255)
         ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTc                                                                                     <  1:512050/100‑1 (MQ=255)
         ggAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTc                                                                                     <  1:760580/100‑1 (MQ=255)
          gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCt                                                                                    <  1:2479997/100‑1 (MQ=255)
          gAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCt                                                                                    <  1:2505605/100‑1 (MQ=255)
           aaTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCTGTTCAGTTCCTTCCGGGGGAGCGTTCTg                                                                                   <  1:2366850/100‑1 (MQ=255)
            aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCTCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGt                                                                                  <  1:1889287/100‑1 (MQ=255)
            aTAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGt                                                                                  <  1:2148470/100‑1 (MQ=255)
              aGACGGTATCGAACGTCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                <  1:2517562/100‑1 (MQ=255)
              aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                >  2:102700/1‑100 (MQ=255)
              aGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTAc                                                                                <  2:818454/100‑1 (MQ=255)
               gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                               >  1:49130/1‑100 (MQ=255)
                 cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                             >  2:1971738/1‑100 (MQ=255)
                 cGGTATCCAACACCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                             >  2:1341259/1‑100 (MQ=255)
                  ggTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAg                                                                            <  2:780983/100‑1 (MQ=255)
                    tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGt                                                                          <  1:103472/100‑1 (MQ=255)
                    tATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGt                                                                          <  1:103449/100‑1 (MQ=255)
                      tCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTg                                                                        >  2:13845/1‑100 (MQ=255)
                       cGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgc                                                                       <  1:1664954/100‑1 (MQ=255)
                        gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcg                                                                      >  1:489221/1‑100 (MQ=255)
                        gAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgcg                                                                      >  1:1254954/1‑100 (MQ=255)
                          aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                    <  1:231514/100‑1 (MQ=255)
                           cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcac                                                                   <  2:2502547/100‑1 (MQ=255)
                            gCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTCGCGcaca                                                                  <  1:1341268/100‑1 (MQ=255)
                             ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 >  1:1303964/1‑100 (MQ=255)
                             ccTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACAt                                                                 <  1:1214728/100‑1 (MQ=255)
                                 aaTACCCGCAGAGCGCGCGGCCAGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtc                                                             >  1:1076243/1‑100 (MQ=255)
                                   tACCCGCAAAGCGCGCAGCCCGCAAAGTGGAACAGCGTGCGAACAGCAGTTCAGTTCCGTCGGGGAGCACGTACTGTACGCAGGTTGCGCACATat                                                               >  2:262679/1‑95 (MQ=255)
                                     cccGCAGAGGGCGCGACCTGCAAAATCAAACAGCGTGCGAACAGCAGTTCAGCTCCTTCCGACGGAGCGTTCTGTACGCAGGGCGCGCACAtcggct                                                            >  2:738069/1‑92 (MQ=255)
                                     cccGCAGAGCGCGCTGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttc                                                          <  1:2034775/99‑1 (MQ=255)
                                     cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCCGTTCCTTCCGGGGGAGCGTTCTGTGCGCAGGTTGCGCACCTAGtctcc                                                          >  2:1097785/1‑97 (MQ=255)
                                     cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         <  1:242378/100‑1 (MQ=255)
                                     cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtcttct                                                         >  1:96704/1‑100 (MQ=255)
                                      ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        <  2:1933391/100‑1 (MQ=255)
                                      ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        >  1:473993/1‑100 (MQ=255)
                                      ccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTg                                                        >  1:1854386/1‑100 (MQ=255)
                                       cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  1:1267237/1‑100 (MQ=255)
                                       cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGc                                                       >  1:2249303/1‑100 (MQ=255)
                                             gcgcgcGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCGGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAg                                                 >  2:508203/1‑100 (MQ=255)
                                                cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGTAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  1:2428657/100‑1 (MQ=255)
                                                cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  1:1868377/100‑1 (MQ=255)
                                                cgcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                              <  1:1381703/100‑1 (MQ=255)
                                                 gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             <  1:13846/100‑1 (MQ=255)
                                                 gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                             <  2:46324/100‑1 (MQ=255)
                                                  cgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCTGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGa                                            >  2:2500742/1‑100 (MQ=255)
                                                   ggCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGAt                                           >  1:110164/1‑100 (MQ=255)
                                                    gCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATa                                          >  2:305726/1‑100 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGTGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        >  1:1866789/1‑100 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTGCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        >  2:1368355/1‑100 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:2103583/100‑1 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:528537/100‑1 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                        <  2:1952443/100‑1 (MQ=255)
                                                       tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                       >  1:1166513/1‑100 (MQ=255)
                                                        gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  1:1515451/1‑100 (MQ=255)
                                                        gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                      >  2:2346625/1‑100 (MQ=255)
                                                         cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGAt                                     <  1:1335287/100‑1 (MQ=255)
                                                          aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    >  1:390524/1‑100 (MQ=255)
                                                          aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                    <  1:2346640/100‑1 (MQ=255)
                                                           aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  1:1352103/1‑100 (MQ=255)
                                                           aaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCa                                   >  1:1691754/1‑100 (MQ=255)
                                                            aaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  <  2:2042529/100‑1 (MQ=255)
                                                            aaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCaa                                  <  1:1695095/100‑1 (MQ=255)
                                                             aTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAAc                                <  1:1161541/100‑1 (MQ=255)
                                                                  aCAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            <  1:2528038/100‑1 (MQ=255)
                                                                  aCAGCGTGCGAACAGCAGGTCAGTGCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCAGATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTc                            <  1:714924/100‑1 (MQ=255)
                                                                   cAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCg                           <  2:946451/100‑1 (MQ=255)
                                                                    aGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCgg                          <  2:2556262/100‑1 (MQ=255)
                                                                     gCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                         >  1:1411928/1‑100 (MQ=255)
                                                                      cGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGaa                        <  2:2009376/100‑1 (MQ=255)
                                                                      cGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGaa                        <  2:290409/100‑1 (MQ=255)
                                                                       gTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGCTACCGATCAAACGCTCGGAAg                       >  2:575293/1‑100 (MQ=255)
                                                                       gTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAg                       <  2:1738032/100‑1 (MQ=255)
                                                                         gCGAACAGTAGTTCAGTTCCTTCCGGGGGGGCTTTCTTCACGCAGGTTGCCCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCg                           <  2:2062599/94‑1 (MQ=255)
                                                                         gCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     <  1:1696562/100‑1 (MQ=255)
                                                                          cGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCg                     >  2:79467/1‑99 (MQ=255)
                                                                           gAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTg                   <  1:1482940/100‑1 (MQ=255)
                                                                              cagcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATCCCGATCACACGCTCGGCAGCGTAAGc                >  1:2136492/1‑100 (MQ=255)
                                                                              cagcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGc                <  2:865412/100‑1 (MQ=255)
                                                                              cagcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGc                <  1:2220127/100‑1 (MQ=255)
                                                                              cagcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGc                >  1:239237/1‑100 (MQ=255)
                                                                               agcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  1:540379/100‑1 (MQ=255)
                                                                                gcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg               <  2:1166220/99‑1 (MQ=255)
                                                                                gcagTTCAGTTCCTTCCGGGCGAGCGTTCTGTACGCAGGTTGCGCACCTAGTCTTCTGCACCGGGGGCGATACCGATCAAACGCTCGGCAGCGTGAGCGa              >  2:591450/1‑100 (MQ=255)
                                                                                 cagttcagttCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAt             >  2:1885241/1‑100 (MQ=255)
                                                                                 cagttcagttCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAt             <  2:2542081/100‑1 (MQ=255)
                                                                                      cagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa        >  2:316664/1‑100 (MQ=255)
                                                                                            cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc  >  1:2424322/1‑100 (MQ=255)
                                                                                                  |                                                                                                
GTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGC  >  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: