Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,889,648 A→G M151T (ATG→ACG)  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,6480AG97.0% 338.1 / ‑0.7 101M151T (ATG→ACG) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (54/44);  total (57/44)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 102 total aligned reads displayed.

AGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAAC  >  NC_000913/2889551‑2889743
                                                                                                 |                                                                                               
aGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTg                                                                                               <  2:1433232/100‑1 (MQ=255)
 ggCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCgtgt                                                                                              >  1:2677429/1‑100 (MQ=255)
   cAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCg                                                                                            >  2:691165/1‑100 (MQ=255)
   cAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGGTCGCGGTGCAGAGTACGTTACGGTTCGTGTCg                                                                                            >  2:1045935/1‑100 (MQ=255)
    aaCAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGt                                                                                           <  2:1577284/100‑1 (MQ=255)
    aaCAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGt                                                                                           >  2:612434/1‑100 (MQ=255)
      cAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGGGCAGCTGCGTCTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTa                                                                                         <  1:2658481/100‑1 (MQ=255)
        gaggCTGAGAAGTCTGTTGCGCCCACTCGTACGCTTCCGCGTGCATCTGCGACTGGTAATGGTTCCAGGTGCACAGTACGTTACGGTTCGTg                                                                                               <  2:2386467/89‑1 (MQ=255)
        gATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGtt                                                                                          <  2:2262413/97‑1 (MQ=255)
        gATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGc                                                                                       <  1:470709/100‑1 (MQ=255)
          tGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTg                                                                                     >  2:1014608/1‑100 (MQ=255)
           gCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGt                                                                                    >  2:540571/1‑100 (MQ=255)
            cTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTc                                                                                   >  2:2694518/1‑100 (MQ=255)
            cTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTc                                                                                   >  1:529010/1‑100 (MQ=255)
             tCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCg                                                                                  >  2:1227168/1‑100 (MQ=255)
              cAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGc                                                                                 >  2:2563306/1‑100 (MQ=255)
              cAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGc                                                                                 >  1:1508317/1‑100 (MQ=255)
              cAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGc                                                                                 >  2:1547765/1‑100 (MQ=255)
               aGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGcc                                                                                <  2:508520/100‑1 (MQ=255)
                gAAATCTTCTTCGCCCACTGGTACGCTTCCGCGTGCAGCTGCGGCTCGAAAGGTGTCGAGGTGCAGAGTACAGTACGGTTCGTGGGGTTAGCTGTCGc                                                                                 >  2:457941/1‑98 (MQ=255)
                 aaaTCTTCTTCGCCCACTCGTACGCTTCCGCGTTCAGCTGCGACTCGAAAGGCCTCGAGACGCAGAGTACGTTGCGGTTGGTGTCGTTAGCTATc                                                                                   >  2:418556/1‑95 (MQ=255)
                  aaTCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAgc                                                                             >  2:448839/1‑100 (MQ=255)
                  aaTCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAgc                                                                             <  2:1871730/100‑1 (MQ=255)
                   aTCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAgcg                                                                            >  1:365041/1‑100 (MQ=255)
                    tcttctTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAgcgc                                                                           <  2:136293/100‑1 (MQ=255)
                     cttcttCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAgcgc                                                                           >  2:1802267/1‑99 (MQ=255)
                     cttcttCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCa                                                                          >  2:1614207/1‑100 (MQ=255)
                       tcttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATc                                                                        >  1:2908040/1‑100 (MQ=255)
                       tcttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATc                                                                        >  1:2666887/1‑100 (MQ=255)
                        cttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCa                                                                       >  1:100988/1‑100 (MQ=255)
                         ttcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCaa                                                                      <  1:2548699/100‑1 (MQ=255)
                          tcGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAg                                                                     >  2:224041/1‑100 (MQ=255)
                             cccACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAgacc                                                                  >  2:2897395/1‑100 (MQ=255)
                              ccACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAgacc                                                                  <  2:588156/99‑1 (MQ=255)
                               cACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAgaccga                                                                >  2:1318935/1‑100 (MQ=255)
                                aCTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAgaccgac                                                               <  2:2018063/100‑1 (MQ=255)
                                 cTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAgaccgacc                                                              <  2:299886/100‑1 (MQ=255)
                                 cTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAgaccgacc                                                              >  1:630606/1‑100 (MQ=255)
                                    gTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGAGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAg                                                           <  2:1629074/100‑1 (MQ=255)
                                        gCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgca                                                       >  2:819014/1‑100 (MQ=255)
                                        gCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgc                                                        <  2:521991/99‑1 (MQ=255)
                                         cTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgcag                                                      >  2:2087453/1‑100 (MQ=255)
                                         cTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgcag                                                      <  1:833867/100‑1 (MQ=255)
                                           tCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgcagca                                                    >  1:1519594/1‑100 (MQ=255)
                                           tCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTGCGTGTCGTCAGCTGTCGCCAGCGCATCAAGACCGACCGAGTgcagca                                                    <  2:1881287/100‑1 (MQ=255)
                                                 tGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGt                                              <  1:691171/100‑1 (MQ=255)
                                                 tGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGt                                              <  2:2695064/100‑1 (MQ=255)
                                                   cAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGtcc                                            >  2:1395017/1‑98 (MQ=255)
                                                   cAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGc                                            >  1:457577/1‑100 (MQ=255)
                                                    aGCTGCGANTCGTAAGGGTTCGAGGTNCAGANTACGTTACGGTTCGTGTGGTTAGCTGTAGCCAGCGCATCAAGACCGACCGAGTGCAGCATCNGGTGc                                            >  1:2803322/1‑99 (MQ=255)
                                                      cTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                         >  2:2837137/1‑100 (MQ=255)
                                                        gCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCgg                                       <  2:2927483/100‑1 (MQ=255)
                                                        gCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCgg                                       >  2:191215/1‑100 (MQ=255)
                                                         tgACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGt                                      >  1:1462334/2‑100 (MQ=255)
                                                          gACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGtt                                     >  1:711900/1‑100 (MQ=255)
                                                          gACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGtt                                     <  1:2507260/100‑1 (MQ=255)
                                                            cTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTg                                   >  2:2210520/1‑100 (MQ=255)
                                                              cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc                                 >  2:2678997/1‑100 (MQ=255)
                                                              cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc                                 >  1:292728/1‑100 (MQ=255)
                                                              cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc                                 <  2:2766487/100‑1 (MQ=255)
                                                               gTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                >  2:2164496/1‑100 (MQ=255)
                                                                 aaGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGtt                              >  1:2130087/1‑100 (MQ=255)
                                                                  aGGGTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTc                             <  1:2694542/100‑1 (MQ=255)
                                                                   gggTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCTCATCAAGACCGAACGAGTGCAGCACCTGGTGCACCGGTTCGACGt                               >  2:2571546/1‑97 (MQ=255)
                                                                   gggTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCt                            >  1:731318/1‑100 (MQ=255)
                                                                   gggTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGATTGACGTTCt                            >  2:1984729/1‑100 (MQ=255)
                                                                    ggTTCGAGGTGCGGAGTACGCTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGCCCGAGTGCAGCATCTGGTGCCCCGGTTTGACGt                               >  2:1858644/1‑96 (MQ=255)
                                                                    ggTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                           >  1:2362298/1‑100 (MQ=255)
                                                                    ggTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                           >  2:2184365/1‑100 (MQ=255)
                                                                    ggTTCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                           >  1:431778/1‑100 (MQ=255)
                                                                       tCGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCttttt                        >  2:2578083/1‑100 (MQ=255)
                                                                        cGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGGAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                       >  2:356759/1‑100 (MQ=255)
                                                                        cGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                       <  2:23368/100‑1 (MQ=255)
                                                                        cGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                       >  2:2390204/1‑100 (MQ=255)
                                                                        cGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                       <  2:1127551/100‑1 (MQ=255)
                                                                        cGAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTc                       <  2:485433/100‑1 (MQ=255)
                                                                         gAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCa                      <  1:87801/100‑1 (MQ=255)
                                                                         gAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCa                      <  1:1839165/100‑1 (MQ=255)
                                                                         gAGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCa                      <  2:910014/100‑1 (MQ=255)
                                                                          aGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAg                     >  2:2499672/1‑100 (MQ=255)
                                                                          aGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAg                     <  2:870691/100‑1 (MQ=255)
                                                                          aGGTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAg                     <  1:2451080/100‑1 (MQ=255)
                                                                            gTGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGaa                   <  2:1923459/100‑1 (MQ=255)
                                                                             tGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAAt                  <  2:954578/100‑1 (MQ=255)
                                                                             tGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAAt                  >  2:2231338/1‑100 (MQ=255)
                                                                             tGCAGAGTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAAt                  <  1:774665/100‑1 (MQ=255)
                                                                                  agTACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGt             <  1:540578/100‑1 (MQ=255)
                                                                                    tACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTgg           <  2:2017880/100‑1 (MQ=255)
                                                                                    tACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTgg           <  1:2732514/100‑1 (MQ=255)
                                                                                    tACGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTgg           <  1:224044/100‑1 (MQ=255)
                                                                                     aCGTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGa          >  2:2887756/1‑100 (MQ=255)
                                                                                      cGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGaa         >  2:1009334/1‑100 (MQ=255)
                                                                                       gTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAAc        >  1:806769/1‑100 (MQ=255)
                                                                                       gTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAAc        <  1:2887782/100‑1 (MQ=255)
                                                                                       gTTACGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAAc        >  1:1434775/1‑100 (MQ=255)
                                                                                           cGGTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa    <  2:946960/100‑1 (MQ=255)
                                                                                            ggTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaaa   <  2:2087917/100‑1 (MQ=255)
                                                                                            ggTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaaa   >  2:1016470/1‑100 (MQ=255)
                                                                                             gTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaatc  >  2:1841400/1‑98 (MQ=255)
                                                                                             gTTCGTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAAc  <  1:2404508/100‑1 (MQ=255)
                                                                                                 |                                                                                               
AGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAAC  >  NC_000913/2889551‑2889743

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: ‑