Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,999,231 T→A 15.8% 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,2310TA15.8% 252.0 / 37.6 101T84S (ACT→TCT) dcyDD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (41/44);  new base A (8/8);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 101 total aligned reads displayed.

GCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATT  >  NC_000913/1999133‑1999329
                                                                                                  |                                                                                                  
gCGATTGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTc                                                                                                   <  1:2519902/100‑1 (MQ=255)
    ttGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGTAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGcc                                                                                               <  1:2152667/100‑1 (MQ=255)
     tGCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCg                                                                                              <  1:451969/100‑1 (MQ=255)
      gCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCgn                                                                                             <  1:1312834/100‑2 (MQ=255)
      gCCGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGc                                                                                             <  1:1251744/100‑1 (MQ=255)
       ccGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGca                                                                                            <  2:2748161/100‑1 (MQ=255)
        cGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcac                                                                                           <  1:2654048/100‑1 (MQ=255)
        cGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcac                                                                                           <  2:136653/100‑1 (MQ=255)
        cGTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcag                                                                                           <  2:1266021/100‑2 (MQ=255)
         gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcaca                                                                                          <  2:652676/100‑1 (MQ=255)
         gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGcaca                                                                                          >  2:485179/1‑100 (MQ=255)
         gTTGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagt                                                                                          <  1:95883/100‑3 (MQ=255)
          ttGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACAt                                                                                         <  1:958027/100‑1 (MQ=255)
           tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGGTTTTCCAGCAGCGCCTCGTAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtcg                                                                                        <  1:785268/100‑5 (MQ=255)
           tGGTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtcg                                                                                        <  1:2304299/100‑5 (MQ=255)
            ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATgg                                                                                       <  1:2363239/100‑1 (MQ=255)
            ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATgg                                                                                       >  2:1090482/1‑100 (MQ=255)
            ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATg                                                                                        >  2:2594242/1‑99 (MQ=255)
            ggTTAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtcgc                                                                                       <  2:1834660/100‑6 (MQ=255)
              ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                     <  1:2142273/100‑1 (MQ=255)
              ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                     <  1:1522869/100‑1 (MQ=255)
              ttAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGtt                                                                                     <  1:1245752/100‑1 (MQ=255)
               tAAATAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtcgctgc                                                                                    <  1:10547/100‑9 (MQ=255)
                aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAg                                                                                   <  1:562665/100‑1 (MQ=255)
                aaaTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGACTGCCGcagtcgctgcc                                                                                   <  1:2438961/100‑10 (MQ=39)
                  aTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGAc                                                                                 >  1:622141/1‑100 (MQ=255)
                  aTAGTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGAc                                                                                 <  2:2450850/100‑1 (MQ=255)
                   tAGTTTTCTGCGGTCGGGCCAACAGGCTTGTCCAGCAGAGCCACGCGGTACTGACCGAGTTTGGCAGAAACTTCGGCAGTCTGCCGCACAGAGTTAGAc                                                                                 >  2:1896371/1‑99 (MQ=255)
                     gTTTTCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTgg                                                                              >  2:2595350/1‑100 (MQ=255)
                      ttttCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGa                                                                             <  1:341794/100‑1 (MQ=255)
                         tCTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCg                                                                          <  2:1292608/100‑1 (MQ=255)
                          cTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGc                                                                         >  1:1187088/1‑100 (MQ=255)
                          cTGCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGc                                                                         >  2:2233342/1‑100 (MQ=255)
                            gCGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGccc                                                                       >  2:2644322/1‑100 (MQ=255)
                             cGGTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGcccc                                                                      <  2:539701/100‑1 (MQ=255)
                              ggTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCg                                                                     >  2:251797/1‑100 (MQ=255)
                              ggTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCg                                                                     >  2:1366570/1‑100 (MQ=255)
                              ggTTGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCg                                                                     >  2:641102/1‑100 (MQ=255)
                                ttGTGCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGc                                                                   <  2:2599276/100‑1 (MQ=255)
                                  gtgCCAATAGGTTTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAg                                                                 <  1:1757556/100‑1 (MQ=255)
                                  gtgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAg                                                                 <  1:1208677/100‑1 (MQ=255)
                                  gtgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAg                                                                 >  2:37022/1‑100 (MQ=255)
                                   tgCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGt                                                                <  2:952942/100‑1 (MQ=255)
                                    gCCAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGATTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTa                                                               >  1:346170/1‑100 (MQ=255)
                                      cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAAt                                                             >  1:2621524/1‑100 (MQ=255)
                                      cAATAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAAt                                                             >  2:526439/1‑100 (MQ=255)
                                        aTAGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCa                                                           >  2:783231/1‑100 (MQ=255)
                                          aGGATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGc                                                         <  2:2714837/100‑1 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        <  1:424718/100‑1 (MQ=255)
                                           ggATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        <  2:1538613/100‑1 (MQ=255)
                                            gATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCg                                                        <  2:2432144/99‑1 (MQ=255)
                                                  ccAGCAGCGCCACGCAGTGCAGACCTAGTTTGGGAGCGAGTGCGGCAGTCTGCCGCACATGGTTAGAGTGGATCGCCCCGGCAGTAATCAGCGTATCgn                                                  <  2:1818606/99‑2 (MQ=255)
                                                   cagcagCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCa                                                <  2:224969/100‑1 (MQ=255)
                                                   cagcagCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTAAGACTGGATCGCCCTGGCAGTAATCAGCGTATCGGCa                                                >  2:1851559/1‑100 (MQ=255)
                                                     gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                              <  1:59477/100‑1 (MQ=255)
                                                     gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                              >  1:2605525/1‑100 (MQ=255)
                                                     gcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCAcc                                              <  2:2088174/100‑1 (MQ=255)
                                                      cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCt                                             >  2:218370/1‑100 (MQ=255)
                                                         cgcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCa                                          <  1:2502174/100‑1 (MQ=255)
                                                          gcCACGCAGTGCAGACCGAGTTTGGCCGCGCCTGCGGCAGTCTGCCGCACAAGGTTAGACTGGATCACCCCCGCAGTAACCAGCGTATCGGCACCCTCa                                          >  2:409796/1‑99 (MQ=255)
                                                          gcCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCAc                                         <  2:2663597/100‑1 (MQ=255)
                                                            cACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGc                                       >  2:2363181/1‑100 (MQ=255)
                                                              cGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCag                                     <  2:1955427/100‑1 (MQ=255)
                                                               gCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                    >  1:1570230/1‑100 (MQ=255)
                                                               gCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCaga                                    <  1:1877014/100‑1 (MQ=255)
                                                                cAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCagag                                   <  1:1410221/100‑1 (MQ=255)
                                                                  gTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGCTCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCa                                 >  2:855103/1‑100 (MQ=255)
                                                                   tGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCAt                                <  1:2595368/100‑1 (MQ=255)
                                                                    gCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                               >  1:2427941/1‑100 (MQ=255)
                                                                    gCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATc                               >  2:2733861/1‑100 (MQ=255)
                                                                      aGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTg                             >  1:127206/1‑100 (MQ=255)
                                                                       gACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGc                            >  1:2263819/1‑100 (MQ=255)
                                                                           gAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCcgcg                        <  1:2562012/100‑1 (MQ=255)
                                                                             gTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCgag                      >  2:2025648/1‑100 (MQ=255)
                                                                               ttGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGaa                    <  1:488876/100‑1 (MQ=255)
                                                                                tGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGaaa                   >  2:2687795/1‑100 (MQ=255)
                                                                                 ctaaccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAt                  >  2:102086/10‑100 (MQ=39)
                                                                                  gCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGCCTGGATCGCCGCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAAtt                 >  1:278148/1‑100 (MQ=255)
                                                                                   cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTc                >  2:2098788/1‑100 (MQ=255)
                                                                                   cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTc                <  1:2257958/100‑1 (MQ=255)
                                                                                   cAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTc                >  1:476264/1‑100 (MQ=255)
                                                                                    aGCGACTGCGGCAGTCTGCCACACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTcc               >  1:2586483/1‑100 (MQ=255)
                                                                                     gCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa              <  2:2218085/100‑1 (MQ=255)
                                                                                     ccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa              >  2:2123620/6‑100 (MQ=255)
                                                                                     ccatgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCa              >  2:39132/6‑100 (MQ=255)
                                                                                      catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             >  1:529032/5‑100 (MQ=255)
                                                                                      catgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             >  2:531654/5‑100 (MQ=255)
                                                                                       gACTGCGGCAGTCTGCCGCACATGGTTAGACTCGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAg             <  2:1268223/99‑1 (MQ=255)
                                                                                       atgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGc            >  2:2100256/4‑100 (MQ=255)
                                                                                       atgtGCGGCAGACTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGc            >  1:1029161/4‑100 (MQ=255)
                                                                                       atgtGCGGCAGACTACCGCACATGGTTAGACTGGATAGCCTCGGCAGTAATAAGCGTATCGGCACCTTAACGCCGATCATCGGCCGCGAGAAATTCCAgg            >  2:391928/4‑99 (MQ=255)
                                                                                        aCTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCt           <  2:379964/100‑1 (MQ=255)
                                                                                           gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc        >  1:256400/1‑100 (MQ=255)
                                                                                           gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc        >  1:1122846/1‑100 (MQ=255)
                                                                                           gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc        >  1:1300787/1‑100 (MQ=255)
                                                                                           gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc        <  2:246445/100‑1 (MQ=255)
                                                                                           gCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTAc        >  2:2140940/1‑100 (MQ=255)
                                                                                            cGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACg       >  2:2520180/1‑100 (MQ=255)
                                                                                             gccAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGt      >  1:1863902/3‑100 (MQ=255)
                                                                                                 gTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAAtt  >  1:1826871/1‑100 (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: