Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,842,032 G→T R67R (CGC→CGA hycG ← hydrogenase 3 and formate hydrogenase complex, HycG subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,842,0320GT99.8% 2179.7 / NA 556R67R (CGC→CGAhycGhydrogenase 3 and formate hydrogenase complex, HycG subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (268/286);  total (269/287)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 595 total aligned reads displayed.

CACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGAACGTTTGA  >  NC_000913/2841934‑2842132
                                                                                                  |                                                                                                    
cacgcacgCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCg                                                                                                    <  1:10893259/101‑1 (MQ=255)
   gcacgcaGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                                                                                                 >  2:183235/1‑101 (MQ=255)
    cacgcaGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGtnnn                                                                                                >  2:8811414/1‑98 (MQ=255)
    cacgcaGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                                                >  1:7004337/1‑101 (MQ=255)
      cgcagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTcnnn                                                                                              >  2:12094861/1‑98 (MQ=255)
      cgcagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTcnnn                                                                                              >  2:12744879/1‑98 (MQ=255)
        nnncgcagCGGATCGCATTGCACGGGTGACCGCGCCGATAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCAt                                                                                            <  2:13605708/98‑1 (MQ=255)
        cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTgnnn                                                                                            >  2:14196807/1‑98 (MQ=255)
        cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCAt                                                                                            <  1:979401/101‑1 (MQ=255)
        cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCAt                                                                                            <  1:9462924/101‑1 (MQ=255)
        cagcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGGTCTGCAt                                                                                            >  2:3845213/1‑101 (MQ=255)
         nnngcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATc                                                                                           <  2:10268504/98‑1 (MQ=255)
         agcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGGTATGCATc                                                                                           >  2:3547293/1‑101 (MQ=255)
            gcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCaaa                                                                                        >  1:4857189/1‑101 (MQ=255)
             cagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAAc                                                                                       <  1:11403341/101‑1 (MQ=255)
             cagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAAc                                                                                       <  1:6373975/101‑1 (MQ=255)
              agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCaannn                                                                                      >  2:13020064/1‑98 (MQ=255)
              agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACa                                                                                      <  1:7880993/101‑1 (MQ=255)
              agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACa                                                                                      >  1:8189094/1‑101 (MQ=255)
               nnncaTCGCATTGCACCGGTGACAGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCATTCTGCATCAAACAg                                                                                     <  2:13154582/97‑1 (MQ=255)
               ggggATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAAAGACTTTAATGCCGAATCGTTCTGCATCAAACAg                                                                                     <  1:11263511/101‑1 (MQ=255)
                gggATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgc                                                                                    <  1:5172668/101‑1 (MQ=255)
                 ggATCGCATTCCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTATTGCCGAATCGTTCTGCATCAAACAgcg                                                                                   <  1:13094071/101‑1 (MQ=255)
                  gATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgnnn                                                                                  >  2:9237415/1‑98 (MQ=255)
                  gATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                                                  >  1:4303384/1‑101 (MQ=255)
                  gATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                                                  <  2:6806/101‑1 (MQ=255)
                   aTCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc                                                                                 >  1:3572639/1‑101 (MQ=255)
                    tCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggcg                                                                                >  1:6104869/1‑101 (MQ=255)
                        aTTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg                                                                            >  2:5510979/1‑101 (MQ=255)
                        aTTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg                                                                            >  1:1544859/1‑101 (MQ=255)
                         ttGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGc                                                                           <  2:4140373/101‑1 (MQ=255)
                          tGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCg                                                                          >  2:2359597/1‑101 (MQ=255)
                           gcccGGGTGACCGCGCCGATAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCCGCATCAAACAGCGGCGAAAGCGt                                                                         <  1:5402492/98‑1 (MQ=255)
                             aCGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCgnnn                                                                       >  2:10011454/1‑98 (MQ=255)
                                 gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAACCAGCGGCGAAAGCGTGGCGaa                                                                   <  2:1049531/101‑1 (MQ=255)
                                 gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaa                                                                   <  1:11800976/101‑1 (MQ=255)
                                   gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaaaa                                                                 >  2:2783820/1‑101 (MQ=255)
                                    aCCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAat                                                                >  2:4731503/1‑101 (MQ=255)
                                       gcgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttc                                                             <  2:5402324/101‑1 (MQ=255)
                                        cgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcg                                                            >  2:3478476/1‑101 (MQ=255)
                                         gcCNGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgn                                                           >  1:3104269/1‑100 (MQ=255)
                                           cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatt                                                         <  2:5659459/101‑1 (MQ=255)
                                             gTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                                                       <  1:3503169/101‑1 (MQ=255)
                                             gTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                                                       <  1:4593265/101‑1 (MQ=255)
                                             gTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                                                       <  1:8670193/101‑1 (MQ=255)
                                             gTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGAGGCGAAAGCGTGGCGAAAatttcgattta                                                       >  1:13413517/1‑100 (MQ=255)
                                               aaaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGc                                                     >  1:14140321/1‑101 (MQ=255)
                                                aaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAattt                                                              <  1:14509764/91‑1 (MQ=255)
                                                aaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                                                    <  2:202347/101‑1 (MQ=255)
                                                 aCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCaa                                                   <  1:6722168/101‑1 (MQ=255)
                                                 aCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCaa                                                   <  1:5521137/101‑1 (MQ=255)
                                                  nnntAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCagc                                                  <  2:9681199/98‑3 (MQ=255)
                                                  cAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAc                                                  <  1:7418966/101‑1 (MQ=255)
                                                   nnnaaaaTATCCGCATCACGCGGTGAAGGAACGACTTTAATGCAGAATCGTTCTGCATAAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAcc                                                 <  2:10898701/98‑1 (MQ=255)
                                                   aGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAcc                                                 <  2:1180680/101‑1 (MQ=255)
                                                   aGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAcc                                                 >  2:3989333/1‑101 (MQ=255)
                                                    gTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCaa                                                   <  1:9901256/98‑1 (MQ=255)
                                                     tAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGt                                               <  2:660493/101‑1 (MQ=255)
                                                      aaaaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCAGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCgct                                              >  1:12200428/1‑99 (MQ=255)
                                                       aaaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTg                                             >  2:6134624/1‑101 (MQ=255)
                                                        nnnatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGc                                            <  2:12465152/98‑1 (MQ=255)
                                                           atCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgc                                         >  1:3312925/1‑101 (MQ=255)
                                                           atCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgc                                         <  1:9032271/101‑1 (MQ=255)
                                                           atCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgc                                         >  1:11510346/1‑101 (MQ=255)
                                                            tCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgcc                                        <  2:3117100/101‑1 (MQ=255)
                                                               gCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgccgcc                                     <  1:846445/101‑1 (MQ=255)
                                                               gCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgccgcc                                     <  2:4587305/101‑1 (MQ=255)
                                                               gCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgccgcc                                     >  1:13209002/1‑101 (MQ=255)
                                                                cATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGccgccg                                    <  1:8862083/101‑1 (MQ=255)
                                                                 aTGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgccgcnnn                                   >  2:13961112/1‑98 (MQ=255)
                                                                 aTGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTGGCAgccgcnnn                                   >  2:9040928/1‑98 (MQ=255)
                                                                    nnncgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGt                                <  2:7473789/98‑1 (MQ=255)
                                                                     cgcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTc                               <  1:12499553/101‑1 (MQ=255)
                                                                      gcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcc                              >  1:12738473/1‑101 (MQ=255)
                                                                       cgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGtnnn                             >  2:7312823/1‑98 (MQ=255)
                                                                        ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCGGCAGTCCag                            >  2:6264750/1‑100 (MQ=255)
                                                                        ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcnnn                            >  2:7748114/1‑98 (MQ=255)
                                                                        gctGAAGGACCGAGTTTAAGGCCGAATCTTGCGGCATCAAACAGCGGCGGAAAGGTGGCGAAAATTTCGCTTTCGCAACCGTTGCAGCCGCCGCAGTCCAc                            <  2:1631751/99‑1 (MQ=255)
                                                                         gTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAcg                           <  1:4703892/101‑1 (MQ=255)
                                                                         gTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAcg                           <  1:266203/101‑1 (MQ=255)
                                                                           gAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAcgcg                         >  1:3892186/1‑101 (MQ=255)
                                                                              ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTa                      <  1:10325069/101‑1 (MQ=255)
                                                                                nnngACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTaaa                    <  2:10698846/98‑1 (MQ=255)
                                                                                 aCGACTTTAATGCCGAATCGTTTTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAAc                   <  1:13907989/101‑1 (MQ=255)
                                                                                 aCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAAc                   >  1:10056399/1‑101 (MQ=255)
                                                                                   gACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTaaannn                 >  2:7925174/1‑98 (MQ=255)
                                                                                    nnnttAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATa                <  2:11510346/98‑1 (MQ=255)
                                                                                    aCTTTACTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCCTGGCGAAAATTTCGATTTCGCACCCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATa                <  1:5875136/101‑1 (MQ=38)
                                                                                    aCTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATa                <  1:10011454/101‑1 (MQ=255)
                                                                                     nnntAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAg               <  2:13367693/98‑1 (MQ=255)
                                                                                     cTTTAATGCCGAATCGTTCTGCATCAAACAGCGACGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTTCAGCCGCCGCAGTCCACGCGGTAAACTTAg               >  1:6862794/1‑101 (MQ=255)
                                                                                         aaTGCCGCATAGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGTAAATTTCGTTTTCGCACCGTTTGCAGTCGCCGCAGTCCACGCGGTAAACATAGGCAg           <  2:4390418/101‑1 (MQ=255)
                                                                                          aTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGa          <  1:14315851/101‑1 (MQ=255)
                                                                                           nnncGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGaa         <  2:8389918/98‑1 (MQ=255)
                                                                                            nnngAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGAAc        <  2:14569635/98‑1 (MQ=255)
                                                                                            gCCGAATCGTTCTGCATCAACCAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGAAc        <  1:191199/101‑1 (MQ=255)
                                                                                            gCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGTTTTCGCAACCGTTTCAGCCGCCGCATTCCACGCGGTAATCTTAGGCatnnn        >  2:11402412/1‑97 (MQ=255)
                                                                                              cGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGTCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGAACGt      <  2:5183690/101‑1 (MQ=255)
                                                                                              cGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGAACGt      <  1:931677/101‑1 (MQ=255)
                                                                                               gAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGAAcnnn     >  2:14816083/1‑98 (MQ=255)
                                                                                                  |                                                                                                    
CACGCACGCAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGCAGAACGTTTGA  >  NC_000913/2841934‑2842132

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