Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 16.8% intergenic (‑181/+127) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA16.8% 602.6 / 95.0 209intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (82/92);  new base A (33/2);  total (115/94)
Fisher's exact test for biased strand distribution p-value = 4.95e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

Only 100 of 223 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                         
 ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                         >  1:3173474/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        >  2:3150108/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        <  2:2845430/100‑1 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                       <  1:2962055/100‑1 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                      >  1:545022/1‑100 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                     >  1:237660/1‑100 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    >  1:686652/1‑100 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGAGTTa                                                                                                    >  2:3425971/1‑100 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                >  2:2478934/1‑100 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                >  2:1669326/1‑100 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                            <  1:3198107/100‑1 (MQ=255)
                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                          >  1:1776977/1‑100 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                      >  1:921884/1‑100 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:1490677/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:3184218/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                     <  1:4138406/100‑2 (MQ=37)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTCTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:2559503/100‑1 (MQ=255)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                    <  1:3362380/100‑1 (MQ=255)
                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                   <  2:3786337/100‑1 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGGTAACGCCTTATCCGGCCTACaa                                                                                >  1:105899/1‑100 (MQ=255)
                            cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                              <  1:3568787/100‑1 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                            <  1:3646792/100‑1 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                           >  1:3664227/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                           <  1:3253302/100‑1 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTCTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                            >  2:3492175/1‑99 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                         <  1:2944035/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                         >  1:1699585/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                         >  1:2927163/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                       <  1:1367573/100‑1 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                       <  1:3900806/100‑1 (MQ=255)
                                      agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                    <  1:3796728/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                   >  1:3424797/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                   >  1:1595894/1‑100 (MQ=255)
                                       agCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                   >  2:1754591/2‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                               <  2:3087705/100‑1 (MQ=255)
                                            aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                              <  1:3314644/100‑1 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                             >  1:2310467/1‑100 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                           <  2:3829710/100‑1 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAGGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                          >  2:1571776/1‑100 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                        >  1:2673898/1‑100 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                        <  2:1950911/98‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                         <  2:2903065/97‑1 (MQ=255)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑Ga                                                   >  1:2452121/1‑100 (MQ=255)
                                                       gTATACAGATTGTTGATAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAt                                                 <  2:2437113/100‑1 (MQ=255)
                                                       gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAt                                                 <  2:1310932/100‑1 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCAGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAt‑t                                               >  2:2497849/1‑100 (MQ=255)
                                                          taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATTCGGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAt‑t‑‑‑tt                                          >  1:2268689/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCg                                        <  2:3722251/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCg                                        <  2:2927146/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCg                                        <  2:2268675/100‑1 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTa                                      <  2:4151717/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGcc                                  >  1:4052532/1‑100 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGcc                                  <  1:2124849/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATGCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGcc                                  <  2:3923117/100‑1 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTg                                >  2:1890898/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTg                                >  2:110380/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTg                                <  1:485548/100‑1 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTg                                >  2:2378451/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTg                                >  2:857927/1‑100 (MQ=255)
                                                                    cgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTg                                >  2:281776/2‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTGCAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGa                               <  2:2304540/100‑1 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGa                               >  1:2729033/1‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGa                               <  1:75442/100‑1 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGa                               >  1:2519310/1‑100 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATa                             <  2:2646751/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATa                             <  1:4047399/100‑1 (MQ=255)
                                                                          aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGc                          <  2:3223122/100‑1 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCg                         >  2:533968/1‑100 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCg                         >  1:3743788/1‑100 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCgg                        >  2:75980/1‑100 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGa                       <  1:1440388/100‑1 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGa                       >  1:871732/1‑100 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGa                       >  1:233793/1‑99 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                    >  1:2047747/1‑100 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTAAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAgcg                    >  1:2585880/1‑100 (MQ=255)
                                                                                  ttGTTTATGCCGGATGCGCCGTTAACGCCTTATCCGGCCTACAAACTCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCa                  <  1:1157764/100‑1 (MQ=255)
                                                                                    ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGC‑A‑TGAG‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCa                  >  1:1294540/7‑98 (MQ=17)
                                                                                     tgtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGC‑A‑GAGA‑T‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATc                >  1:3992390/6‑99 (MQ=11)
                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCa               >  2:4062215/5‑99 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑A‑AGAT‑T‑‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCAg              >  2:1066045/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTACAAGAACGTGCAAATTCAACAAATTGTGATGCAT‑C‑‑‑CT‑GTAGGCCTGATAAGCGTAGCGCATCa               >  1:1023285/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGC‑A‑‑ATA‑T‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:3794943/4‑100 (MQ=18)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg             >  1:2918897/1‑100 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGC‑‑‑‑AAT‑TCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg             >  1:3687641/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGC‑A‑GGGA‑C‑‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGc            >  1:2551347/3‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGC‑ATGGGC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:3332077/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGT‑G‑GGCT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1516650/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑A‑GGAG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:3830939/2‑100 (MQ=18)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:1754180/1‑100 (MQ=255)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  2:933234/1‑100 (MQ=255)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGC‑A‑GAAA‑T‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  1:1641078/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTC‑TGCAAATTCAACAAATTGC‑‑‑GAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  1:743393/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGA‑T‑CGAT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:3517283/1‑99 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGC‑A‑GGAA‑A‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCag          >  1:2784636/1‑99 (MQ=18)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:1699575/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:919456/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:3608705/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   >  2:2098187/1‑100 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  1:3401998/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:3730951/1‑100 (MQ=255)
                                                                                                  |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑