Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,295,417 G→A P510L (CCG→CTG)  fdhF ← selenopolypeptide subunit of formate dehydrogenase H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,295,4170GA99.4% 1247.5 / ‑5.0 344P510L (CCG→CTG) fdhFselenopolypeptide subunit of formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (184/158);  total (185/159)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 353 total aligned reads displayed.

AAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295325‑4295506
                                                                                            |                                                                                         
  gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         >  1:6501725/1‑93 (MQ=255)
  gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         <  2:4910644/93‑1 (MQ=255)
   tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        >  1:6103974/1‑93 (MQ=255)
   tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        <  2:6083795/93‑1 (MQ=255)
    ccccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcc                                                                                       >  1:1068400/1‑92 (MQ=255)
    ccccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGca                                                                                       <  2:1943071/93‑1 (MQ=255)
    ccccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGca                                                                                       >  2:6301006/1‑93 (MQ=255)
      ccTGATCGGCATCTGAAGTATCGCGGCAGGGCCACTGAATGAAGCCCAGTTCGCCCATTATCTCGTAAGTCGCACCGTAGAAATCCaggacca                                                                                     >  1:7630413/1‑89 (MQ=255)
       cTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAg                                                                                    >  2:1715583/1‑93 (MQ=255)
        tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  2:1490199/1‑93 (MQ=255)
        tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  1:4563175/1‑93 (MQ=255)
        tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  2:6063084/1‑93 (MQ=255)
        tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  1:2235341/1‑93 (MQ=255)
          aTCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGAta                                                                                 >  1:1481188/1‑92 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACTGTAGAAATCCAGGCACAGATGa                                                                                <  2:3480885/93‑1 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  1:1251178/93‑1 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                >  2:7927512/1‑93 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  1:518003/93‑1 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  2:5582243/93‑1 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                >  1:7331370/1‑93 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  1:1985607/93‑1 (MQ=255)
           tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  1:789969/93‑1 (MQ=255)
             ggCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACg                                                                              >  1:4725523/1‑93 (MQ=255)
             ggCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACg                                                                              >  2:422892/1‑93 (MQ=255)
              gCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGc                                                                             >  1:1955836/1‑93 (MQ=255)
                aTCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCaa                                                                           >  2:511540/1‑93 (MQ=255)
                 tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          >  2:4740325/1‑93 (MQ=255)
                   tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        >  2:2865438/1‑93 (MQ=255)
                   tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCCCCCATTTTCTCGTAAGCCGCACCGTAGAAACCCAGCAACAGATGAGGTAACTc                                                                        >  2:2165559/1‑93 (MQ=255)
                       gTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATcc                                                                    <  2:1340633/93‑1 (MQ=255)
                        tATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   <  2:3582931/93‑1 (MQ=255)
                        tATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   >  2:4672200/1‑93 (MQ=255)
                              ggCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtct                                                             >  1:1367134/1‑93 (MQ=255)
                                 aaGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCt                                                          <  2:1825026/93‑1 (MQ=255)
                                  aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                         >  2:2951502/1‑93 (MQ=255)
                                  aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                         >  2:769656/1‑93 (MQ=255)
                                  aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                         <  1:6906645/93‑1 (MQ=255)
                                    gCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTggg                                                       <  1:2837225/93‑1 (MQ=255)
                                      cACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtg                                                     <  1:3927193/93‑1 (MQ=255)
                                      cACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtg                                                     >  2:24705/1‑93 (MQ=255)
                                       aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    >  2:7586953/1‑93 (MQ=255)
                                       aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    >  2:1029225/1‑93 (MQ=255)
                                        cTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGcgtt                                                   >  1:5712024/1‑93 (MQ=255)
                                             tgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTTCTGGGtgttgttgt                                              >  1:5933445/1‑93 (MQ=255)
                                             tgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCCCCTGGGtgttgttgt                                              <  1:7741085/93‑1 (MQ=255)
                                              gaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACACATCCCAGATCTCCTGGGTGTTGTTGTa                                             <  2:7231088/93‑1 (MQ=255)
                                               aaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAg                                            <  2:1437149/93‑1 (MQ=255)
                                                 gCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                          >  1:6723974/1‑93 (MQ=255)
                                                   ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCa                                        >  2:7765306/1‑93 (MQ=255)
                                                   ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCa                                        >  1:5154701/1‑93 (MQ=255)
                                                   ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCa                                        >  1:2726791/1‑93 (MQ=255)
                                                   ccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTACa                                        >  2:3119109/1‑93 (MQ=255)
                                                    cAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCAt                                       >  2:394748/1‑93 (MQ=255)
                                                     aGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                      <  2:1652566/93‑1 (MQ=255)
                                                      gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  1:4097509/1‑93 (MQ=255)
                                                      gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  2:6409338/1‑93 (MQ=255)
                                                      gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  2:700141/1‑93 (MQ=255)
                                                      gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  1:3075194/1‑93 (MQ=255)
                                                       ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    >  2:893602/1‑93 (MQ=255)
                                                       ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    >  1:6135056/1‑93 (MQ=255)
                                                        tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  1:3543520/93‑1 (MQ=255)
                                                          gCCCATTTTCTCTTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGTTGTTGTTGTTGTGCATCGGATa                                 >  2:5660844/1‑93 (MQ=255)
                                                          gCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATa                                 <  1:342377/93‑1 (MQ=255)
                                                            ccATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAc                               <  2:4563389/93‑1 (MQ=255)
                                                            ccATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAc                               >  1:5790002/1‑93 (MQ=255)
                                                               ttttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCa                            >  1:3670251/1‑93 (MQ=255)
                                                                 ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATa                          <  2:6769141/93‑1 (MQ=255)
                                                                 ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATa                          <  1:5314151/93‑1 (MQ=255)
                                                                  tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         <  2:7080333/93‑1 (MQ=255)
                                                                  tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         <  1:7492589/93‑1 (MQ=255)
                                                                  tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         >  1:6541912/1‑93 (MQ=255)
                                                                    tcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgg                       >  1:131617/1‑93 (MQ=255)
                                                                     cGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACggg                      <  1:930710/93‑1 (MQ=255)
                                                                     cGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACggg                      <  1:6272524/93‑1 (MQ=255)
                                                                     cGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACggg                      <  2:4498838/93‑1 (MQ=255)
                                                                      gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     <  2:532258/93‑1 (MQ=255)
                                                                      gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     <  1:2553304/93‑1 (MQ=255)
                                                                       tAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTg                    >  2:5888018/1‑93 (MQ=255)
                                                                          gTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCg                 >  1:4964334/1‑93 (MQ=255)
                                                                          gTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCg                 >  1:7034406/1‑93 (MQ=255)
                                                                          gTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCg                 <  2:6468497/93‑1 (MQ=255)
                                                                            cGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAt               >  1:5872870/1‑93 (MQ=255)
                                                                             gCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAtt              <  1:4299717/93‑1 (MQ=255)
                                                                               cccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            <  1:2496096/92‑1 (MQ=255)
                                                                               aCCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            <  2:3184212/93‑1 (MQ=255)
                                                                               aCCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            >  2:6225818/1‑93 (MQ=255)
                                                                                ccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCa           >  1:5486314/1‑93 (MQ=255)
                                                                                ccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCa           <  1:1628668/93‑1 (MQ=255)
                                                                                 cGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCAc          >  1:2285904/1‑93 (MQ=255)
                                                                                  gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCAcg         >  1:6007658/1‑92 (MQ=255)
                                                                                  gTAGAAATCCAGGCACAGATGACGCAACTCATCACAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         >  2:2376193/1‑93 (MQ=255)
                                                                                   tAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtg        >  1:3750593/1‑93 (MQ=255)
                                                                                    aGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga       <  1:798329/93‑1 (MQ=255)
                                                                                     gAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgat      >  1:5396713/1‑93 (MQ=255)
                                                                                      aaaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatg     <  1:4990290/93‑1 (MQ=255)
                                                                                       aaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATGCGGGTGGCGATTTCACtgatga    <  2:3277043/93‑1 (MQ=255)
                                                                                       aaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatga    <  2:119927/93‑1 (MQ=255)
                                                                                        aTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTgatgat   >  2:254889/1‑93 (MQ=255)
                                                                                         tCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAtt  >  1:3221978/1‑93 (MQ=255)
                                                                                         tCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAtt  >  2:4393403/1‑93 (MQ=255)
                                                                                            |                                                                                         
AAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295325‑4295506

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