Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,889,702 A→G 21.0% V133A (GTG→GCG)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,7020AG21.0% 358.9 / 91.1 195V133A (GTG→GCG) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base A (81/72);  new base G (25/16);  total (107/88)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 197 total aligned reads displayed.

AGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGT  >  NC_000913/2889603‑2889801
                                                                                                   |                                                                                                   
aGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCCTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCa                                                                                                     >  1:2510920/1‑100 (MQ=255)
    tcGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCgg                                                                                                 >  2:1246981/2‑100 (MQ=255)
         tCGTAAGGGTTCGTGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGa                                                                                            <  2:314301/100‑1 (MQ=255)
         tCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGa                                                                                            <  1:1914820/100‑1 (MQ=255)
         tCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGa                                                                                            <  1:1819013/100‑1 (MQ=255)
          cGTAATGGTTCGAGTTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAGGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGa                                                                                            <  1:3438203/99‑1 (MQ=255)
          cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGAc                                                                                           >  2:1555916/1‑100 (MQ=255)
          cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc                                                                                           >  2:2651755/1‑100 (MQ=255)
          cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAAACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc                                                                                           >  1:84119/1‑100 (MQ=255)
          cGTAAGGGTTCGAGGTGAAGAGTACGTTACGGATCATGTCATTAGCTGTCGCCAGCGCATCAAGACAGACCGAGTGCAGCATCTGGAGCACCGGTGTGAc                                                                                           >  2:2452228/1‑100 (MQ=255)
           gTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                          >  2:2853788/1‑100 (MQ=255)
           gTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                          <  1:615054/100‑1 (MQ=255)
            tAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                          >  1:672492/1‑99 (MQ=255)
             aaGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGtt                                                                                        >  1:2676249/1‑100 (MQ=255)
             aaGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCCCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACggt                                                                                        >  2:1869345/1‑98 (MQ=255)
              aGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTc                                                                                       <  2:3681374/100‑1 (MQ=255)
                ggTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                                                                                     <  1:3653488/100‑1 (MQ=255)
                  ttCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCtttt                                                                                   >  2:2880512/1‑100 (MQ=255)
                    cGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTc                                                                                 <  1:3388517/100‑1 (MQ=255)
                    cGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                                                                                 <  1:356470/100‑1 (MQ=255)
                      aGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAg                                                                               <  1:3263993/100‑1 (MQ=255)
                       ggTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGa                                                                              >  2:1063973/1‑100 (MQ=255)
                        gTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTccg                                                                               <  1:2659040/98‑3 (MQ=255)
                          gCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATg                                                                           >  2:2187494/1‑100 (MQ=255)
                             gagTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCg                                                                        >  2:1457656/1‑100 (MQ=255)
                                 aCGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGa                                                                    <  2:80417/100‑1 (MQ=255)
                                 aCGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGa                                                                    >  1:598287/1‑100 (MQ=255)
                                  cGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGaa                                                                   <  1:1976841/100‑1 (MQ=255)
                                  cGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGaa                                                                   <  2:2844465/100‑1 (MQ=255)
                                  cGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGaa                                                                   <  1:2141510/100‑1 (MQ=255)
                                  cGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGaa                                                                   >  1:2976934/1‑100 (MQ=255)
                                    ttACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACt                                                                 <  2:1395609/100‑1 (MQ=255)
                                       cGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa                                                              <  2:2733532/100‑1 (MQ=255)
                                       cGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa                                                              >  1:1080609/1‑100 (MQ=255)
                                        ggTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaaa                                                             >  2:1162357/1‑100 (MQ=255)
                                          ttCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACg                                                           <  1:3715190/100‑1 (MQ=255)
                                            cATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAAAATGCCGTGAAACTGAAACGTc                                                         >  1:922242/1‑100 (MQ=255)
                                              tGTCGTTAGCTGTCGCCAGCGCATCAAGACCGATCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTg                                                       <  1:472168/100‑1 (MQ=255)
                                              tGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTg                                                       >  1:3703110/1‑100 (MQ=255)
                                                tCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGc                                                     <  2:7031/100‑1 (MQ=255)
                                                tCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGc                                                     <  1:791642/100‑1 (MQ=255)
                                                  gTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCgg                                                   <  1:2987216/100‑1 (MQ=255)
                                                    tAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGtt                                                 >  2:2620226/1‑100 (MQ=255)
                                                    tAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGtt                                                 >  2:2152098/1‑100 (MQ=255)
                                                    ctgCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGtt                                                 <  1:3078276/98‑1 (MQ=255)
                                                      gCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTgg                                               >  2:1938307/1‑100 (MQ=255)
                                                      gCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCATCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTgg                                               <  1:1997940/100‑1 (MQ=255)
                                                      gCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAGACGTCTGGCGGTTgg                                               >  1:3869221/1‑100 (MQ=255)
                                                      gCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGttag                                               >  2:3078248/1‑98 (MQ=255)
                                                        tGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGTACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTc                                             <  1:323218/100‑1 (MQ=255)
                                                         tctgCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCa                                            <  2:485396/97‑1 (MQ=255)
                                                         gTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCa                                            >  1:2032609/1‑100 (MQ=255)
                                                           cGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCCCCGGTTTGCCGTTCTTTTTCAGAATGCGGTGGAACTGAAACGTCTGGCGGTTGGTCAgg                                          >  1:2685033/1‑100 (MQ=255)
                                                           cGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATTCCGTGAAACTGAAACGTCTGGCGGTTGGTCAg                                           >  1:3324681/1‑99 (MQ=255)
                                                            gCCAGCGCATCAAGACCGACTGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTGTTTCAGAATGCCGTGGAAGTGAAACGTCTGGCGGTTGGTCAAGc                                         <  2:1029482/100‑1 (MQ=255)
                                                            gCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCTGTTTGACGTTCTTTTTCAGAATGCCTTGGAACTGAAACGTCTGGCGGTTGGTCAGGc                                         >  1:1386323/1‑100 (MQ=255)
                                                            gCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGATTTGACGTTCTTTTTCAGAATGACGTGGAACTGAAACGTCTGGCGGTTGGTCAGGc                                         >  1:239222/1‑100 (MQ=255)
                                                            gCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGc                                         <  2:25007/100‑1 (MQ=255)
                                                             cttccgatctaaGTCCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCg                                        <  1:2273070/90‑1 (MQ=255)
                                                             ccAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACTGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCg                                        <  1:3171919/100‑1 (MQ=255)
                                                               aGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCTGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCTGTTGGTCAgg                                          <  2:3655602/96‑1 (MQ=255)
                                                                 cgcATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtg                                    <  2:615455/100‑1 (MQ=255)
                                                                 cgcATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtg                                    <  2:3136034/100‑1 (MQ=255)
                                                                   cATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTTGAACTAAAACGTCTGGCGGTTGGTCAGGCGAAtg                                    <  2:236895/98‑1 (MQ=255)
                                                                    aTCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgctg                                 >  1:2661290/1‑100 (MQ=255)
                                                                     tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTTGTCAGGCGAAtgctgc                                <  2:1771391/100‑1 (MQ=255)
                                                                     tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgctgc                                <  1:1947722/100‑1 (MQ=255)
                                                                     tCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgctgc                                <  1:3760343/100‑1 (MQ=255)
                                                                     tCAAGACCGACCGAGTGCAACATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgctgc                                >  2:3764479/1‑100 (MQ=255)
                                                                      cAAGACCGACCGAGTGCAGCATCTGTTTCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGCACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgctgac                               >  2:45544/1‑98 (MQ=255)
                                                                      cAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGcc                               <  2:3818431/100‑1 (MQ=255)
                                                                      cAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGcc                               >  2:3598859/1‑100 (MQ=255)
                                                                      cAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAAGCGAATGCTGcc                               >  2:3369776/1‑100 (MQ=255)
                                                                            cgaccgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGAt                         <  2:9003/100‑1 (MQ=255)
                                                                             gaccgaGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCATGTAACTGAAACGTCTGGCGGTTGGTCAGGCGAATACTTCCATAGAt                         >  1:3127494/1‑99 (MQ=255)
                                                                             gaccgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGTTTGGTCAGGCGAATGCTGACATAGAt                         <  2:514182/99‑1 (MQ=255)
                                                                              accgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATgg                       <  2:3188472/100‑1 (MQ=255)
                                                                                cgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATaagatcgg                     >  1:103544/1‑93 (MQ=255)
                                                                                 gaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGgtgt                    >  2:291272/1‑100 (MQ=255)
                                                                                   gTGCAGCATCTGGTGCACCGGTTTGACTTTCTTTTTCAGAATGCCTTTGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGttt                  <  1:1143257/100‑1 (MQ=255)
                                                                                   gTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTGAGGCGAATGCTGCCATAGATGGTGttt                  <  1:2853814/100‑1 (MQ=255)
                                                                                   gTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGttt                  <  1:1787042/100‑1 (MQ=255)
                                                                                   gTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGttt                  <  2:2247387/100‑1 (MQ=255)
                                                                                    tGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGCCGAATGCTGCCATAGATGGCGtttt                 >  2:1320135/1‑100 (MQ=255)
                                                                                     tctgcaTCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTc                <  1:1602678/97‑1 (MQ=255)
                                                                                     gaagcaTCTGGTGAACCGGTTTGACGTTCTTATTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTc                >  1:2640915/3‑100 (MQ=255)
                                                                                        gcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAATTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc             <  1:551805/100‑1 (MQ=255)
                                                                                        gcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTTAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc             >  2:198638/1‑100 (MQ=255)
                                                                                        gcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc             >  1:1203324/1‑100 (MQ=255)
                                                                                          aTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCgg           >  2:1992250/1‑100 (MQ=255)
                                                                                          aTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGAGTTTTCACCgg           >  1:1910110/1‑100 (MQ=255)
                                                                                           tCTGGTGCGCCGGTTTGACGTTCTTTTTCAGATTTCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAACg            >  2:2629315/1‑98 (MQ=255)
                                                                                            cTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCa         <  1:2697224/100‑1 (MQ=255)
                                                                                              ggTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTACCATAGATGGTGTTTTCACCGGCaa        >  1:1041337/1‑99 (MQ=255)
                                                                                                tGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAAtt     >  2:2774592/1‑100 (MQ=255)
                                                                                                tGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAAtt     >  2:3296388/1‑100 (MQ=255)
                                                                                                 gCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCACAttt    >  2:3647257/1‑100 (MQ=255)
                                                                                                 gCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAAttt    >  2:143536/1‑100 (MQ=255)
                                                                                                 gCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGAGGTTGGTCAGGCGAATGCTGCAATAGATGGTGTTTTCACCGGCAAAttt    >  2:3913867/1‑100 (MQ=255)
                                                                                                  cACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTg   >  2:2237314/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGT  >  NC_000913/2889603‑2889801

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