Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,999,231 T→A 18.6% T84S (ACT→TCT)  dcyD ← D‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,999,2310TA18.6% 228.4 / 39.9 101T84S (ACT→TCT) dcyDD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (40/40);  new base A (10/8);  total (52/49)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 101 total aligned reads displayed.

GCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATT  >  NC_000913/1999133‑1999329
                                                                                                  |                                                                                                  
gCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTc                                                                                                   <  2:1666909/100‑1 (MQ=255)
   aTTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGc                                                                                                <  2:1451696/100‑1 (MQ=255)
   aTTGCCGTTGGTTAAATAGTTTTCTGCGGGTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGc                                                                                                <  1:1104267/100‑1 (MQ=255)
    ttGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGcc                                                                                               <  1:1918910/100‑1 (MQ=255)
      gCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGc                                                                                             >  2:136417/1‑100 (MQ=255)
       ccGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGca                                                                                            <  2:1670596/100‑1 (MQ=255)
         gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagt                                                                                          <  1:854037/100‑3 (MQ=255)
         gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagt                                                                                          <  2:324105/100‑3 (MQ=255)
          ttGGTTTAATAGTTTTATGCGGTTGTGCCAATCGGATTTTTCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagt                                                                                          <  1:2305739/99‑3 (MQ=255)
          ttGGTTAAATAGTTTTCTGCGGTTGTGCCTATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtc                                                                                         <  1:2142357/100‑4 (MQ=255)
          ttGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACAt                                                                                         <  1:2584966/100‑1 (MQ=255)
          ttGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtc                                                                                         <  1:587560/100‑4 (MQ=255)
          ttGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtc                                                                                         <  1:702428/100‑4 (MQ=255)
          ttGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtc                                                                                         <  1:816427/100‑4 (MQ=255)
           tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATg                                                                                        >  1:2095041/1‑100 (MQ=255)
            ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACAt                                                                                         >  1:1875244/1‑98 (MQ=255)
            ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtc                                                                                         <  1:437922/98‑4 (MQ=255)
              ttAAATAGTTTTCTGGGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                     <  2:2554856/100‑1 (MQ=255)
                aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAg                                                                                   >  1:862411/1‑100 (MQ=255)
                aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAg                                                                                   >  1:2162657/1‑100 (MQ=255)
                 aaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGa                                                                                  <  2:709235/100‑1 (MQ=255)
                  aTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGCCTGCGGCAGTCTGCCGCACATGGTTAGAc                                                                                 >  2:2181798/1‑100 (MQ=255)
                  aTAGGTTTCTGCGGTTGTGCCAATAGGGTTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGAc                                                                                 >  1:1511237/1‑100 (MQ=255)
                      ttttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGa                                                                             <  2:37929/100‑1 (MQ=255)
                       tttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                            >  2:746595/1‑100 (MQ=255)
                       tttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                            <  2:862414/100‑1 (MQ=255)
                       tttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGAt                                                                            <  2:2162671/100‑1 (MQ=255)
                        ttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATc                                                                           <  1:1370167/100‑1 (MQ=255)
                        tgcggCGGTTGTGCCAATAGGGTTTTCCAGCCGCGCGACGCAGTGCAGACCGAGTGTGGCAGCGACTGCGGCAGTCTGCCGTACATGGTTAGACTGGATc                                                                           <  1:1156445/96‑1 (MQ=255)
                             cGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcccc                                                                      <  1:2482878/100‑1 (MQ=255)
                              ggTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGCGTTTGACAGCGACTGCACCAGCCTCCCGCACATAGTGAGACTGGCTc                                                                           >  1:2572018/1‑94 (MQ=255)
                                ttGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                   <  2:67798/100‑1 (MQ=255)
                                   tgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGt                                                                >  2:1616713/1‑100 (MQ=255)
                                    gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTa                                                               >  2:1888141/1‑100 (MQ=255)
                                    gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAATTTGGCAGCGACTGCGGCACTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTa                                                               >  1:1415259/1‑100 (MQ=255)
                                      cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAAt                                                             >  2:121248/1‑100 (MQ=255)
                                      cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAAt                                                             >  2:1767897/1‑100 (MQ=255)
                                      cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAAt                                                             >  2:1767835/1‑100 (MQ=255)
                                       aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATc                                                            >  1:689623/1‑100 (MQ=255)
                                       aaTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTACGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATc                                                            >  1:182467/1‑100 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        <  1:314640/100‑1 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        <  1:282133/100‑1 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        >  2:2279172/1‑100 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        <  1:929793/100‑1 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        >  2:1462864/1‑100 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        >  1:970669/1‑100 (MQ=255)
                                            gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGt                                                       >  2:177241/1‑100 (MQ=255)
                                             aTTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTa                                                      <  1:884465/100‑1 (MQ=255)
                                             aTTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTa                                                      >  2:115120/1‑100 (MQ=255)
                                              ttttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTAt                                                     >  2:911253/1‑100 (MQ=255)
                                              ttttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTAt                                                     >  1:579897/1‑100 (MQ=255)
                                               tttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATc                                                    <  1:201514/100‑1 (MQ=255)
                                               tttCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATc                                                    <  1:25032/100‑1 (MQ=255)
                                                 tCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCgg                                                  >  2:1069090/1‑100 (MQ=255)
                                                     gcagcGCCACGCAGTGCAGGCCGAGTTGGGCAGCGGCTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                              <  1:80833/100‑1 (MQ=255)
                                                     gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                              <  2:764828/100‑1 (MQ=255)
                                                     gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                              <  1:530386/100‑1 (MQ=255)
                                                     gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                              <  1:803123/100‑1 (MQ=255)
                                                      cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCt                                             >  2:1333468/1‑100 (MQ=255)
                                                      cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCt                                             >  1:1617183/1‑100 (MQ=255)
                                                       agcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCtt                                            <  2:754648/100‑1 (MQ=255)
                                                         cgcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCa                                          <  2:2026279/100‑1 (MQ=255)
                                                          gcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                         >  1:413407/1‑100 (MQ=255)
                                                            cACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGc                                       <  1:868187/100‑1 (MQ=255)
                                                             aCGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCa                                      >  2:2047573/1‑100 (MQ=255)
                                                              cGCAGTGCAGACCGAGCTTGGCAGCGACTGCGGCCGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCag                                     >  2:1795927/1‑100 (MQ=255)
                                                                cAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCagag                                   <  1:2544398/100‑1 (MQ=255)
                                                                   tGCAGACCGAGTTTGGCAGCGACTGCGGCAGCCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCAt                                >  1:1150181/1‑100 (MQ=255)
                                                                    gCAGACCGAGTTTGGCAGCGACTGCCGCAGTCTGCCGCACATTGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                               <  1:2065770/100‑1 (MQ=255)
                                                                     cAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCt                              >  2:919245/1‑100 (MQ=255)
                                                                      aGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTg                             >  2:1193672/1‑100 (MQ=255)
                                                                      aGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTg                             >  2:1801864/1‑100 (MQ=255)
                                                                      aGACCGAGTTTGGCAGCGACTGCGGCAGCCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTg                             <  2:1150185/100‑1 (MQ=255)
                                                                       gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGc                            >  2:1618489/1‑100 (MQ=255)
                                                                        aCCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGcc                           >  2:59320/1‑100 (MQ=255)
                                                                         ccGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcg                          >  2:972907/1‑100 (MQ=255)
                                                                           gAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcgcg                        <  2:990278/100‑1 (MQ=255)
                                                                           gAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcgcg                        >  2:1520744/1‑100 (MQ=255)
                                                                               ttGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGaa                    >  1:922448/1‑100 (MQ=255)
                                                                                 ggCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAt                  >  2:345138/1‑100 (MQ=255)
                                                                                 ctaaccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAt                  >  2:472155/10‑100 (MQ=39)
                                                                                  gCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAtt                 >  1:938535/1‑100 (MQ=255)
                                                                                   aaccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTc                >  1:1259674/8‑100 (MQ=255)
                                                                                    aGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc               <  2:808408/100‑1 (MQ=255)
                                                                                      catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             >  2:1951799/5‑100 (MQ=255)
                                                                                      catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             >  2:2445457/5‑100 (MQ=255)
                                                                                      catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             >  2:1728893/5‑100 (MQ=255)
                                                                                      catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             >  1:247936/5‑100 (MQ=255)
                                                                                      cGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             <  2:615862/100‑1 (MQ=255)
                                                                                       atgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGc            >  2:2432933/4‑100 (MQ=255)
                                                                                        tgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGGCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCt           >  1:1382782/3‑100 (MQ=255)
                                                                                        tgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCt           >  1:1849705/3‑100 (MQ=255)
                                                                                         cTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCtt          <  2:970672/100‑1 (MQ=255)
                                                                                         cTGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCtt          >  2:1700050/1‑100 (MQ=255)
                                                                                          tGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTa         >  2:887653/1‑100 (MQ=255)
                                                                                           gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc        <  1:2105639/100‑1 (MQ=255)
                                                                                           gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc        <  2:1946016/100‑1 (MQ=255)
                                                                                            cGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACg       <  2:1270327/100‑1 (MQ=255)
                                                                                             ggCAGTCTGCCGCCCATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGt      <  2:820775/100‑1 (MQ=255)
                                                                                                 gTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTGCAGCTTACGTAAtt  <  1:18179/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
GCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATT  >  NC_000913/1999133‑1999329

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: