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.8% 1465.5 / ‑5.6 396P510L (CCG→CTG) fdhFselenopolypeptide subunit of formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (215/180);  total (215/181)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 404 total aligned reads displayed.

GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506
                                                                                          |                                                                                         
gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         >  1:4543255/1‑93 (MQ=255)
gTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCagg                                                                                         <  2:3383565/93‑1 (MQ=255)
 tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        <  1:6589838/93‑1 (MQ=255)
 tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        >  1:6279435/1‑93 (MQ=255)
 tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        >  1:6164943/1‑93 (MQ=255)
 tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        >  1:7371544/1‑93 (MQ=255)
 gccccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        <  2:5451986/92‑1 (MQ=255)
  ccccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGca                                                                                       >  2:6885741/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTGTAGCGACAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcacaaa                                                                                   <  1:8983123/93‑3 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  2:3598974/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  1:69191/1‑93 (MQ=255)
      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   <  1:1941114/93‑1 (MQ=255)
        aTCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATg                                                                                 >  2:7344870/1‑93 (MQ=255)
        aTCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATg                                                                                 >  2:9579780/1‑93 (MQ=255)
         tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATgt                                                                                >  2:2068671/1‑92 (MQ=255)
         tCGGCATCTGAAGTATCGAGGCAAGGGCCCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                >  1:6104184/1‑93 (MQ=255)
           ggCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACg                                                                              >  2:2476410/1‑93 (MQ=255)
            gCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGc                                                                             >  2:326679/1‑93 (MQ=255)
             cATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCa                                                                            >  2:4648135/1‑93 (MQ=255)
             cATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCa                                                                            >  1:1747814/1‑93 (MQ=255)
              aTCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCaa                                                                           >  2:3958522/1‑93 (MQ=255)
               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          >  1:9971211/1‑93 (MQ=255)
               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          <  1:2348962/93‑1 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        >  2:9216059/1‑93 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        <  2:7014264/93‑1 (MQ=255)
                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        <  2:1536399/93‑1 (MQ=255)
                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       <  2:4484005/93‑1 (MQ=255)
                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       <  1:7680426/93‑1 (MQ=255)
                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       <  1:1538258/93‑1 (MQ=255)
                    aGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATc                                                                     >  1:5028102/1‑93 (MQ=255)
                    aGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATc                                                                     <  2:7995954/93‑1 (MQ=255)
                      tATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   <  1:4815810/93‑1 (MQ=255)
                      tATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   <  1:421045/93‑1 (MQ=255)
                      gaTCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   <  1:406365/92‑1 (MQ=255)
                        tCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAg                                                                 >  2:4998950/1‑93 (MQ=255)
                         cgcgGCAAGGCCACTGAATGAAGCTCAGTGCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGa                                                                <  2:230493/93‑1 (MQ=255)
                          gcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAt                                                               <  2:4888951/93‑1 (MQ=255)
                          gcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAt                                                               >  2:7531583/1‑93 (MQ=255)
                           cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtc                                                              <  2:9390159/93‑1 (MQ=255)
                            ggCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtct                                                             >  1:9573865/1‑93 (MQ=255)
                             gCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtctc                                                            >  1:3457609/1‑93 (MQ=255)
                             gCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtctc                                                            >  1:1163095/1‑93 (MQ=255)
                                 ggCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                        >  2:6396606/1‑93 (MQ=255)
                                 ggCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                        >  2:8176525/1‑93 (MQ=255)
                                      cTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgtt                                                   >  2:9852119/1‑93 (MQ=255)
                                       tgaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttg                                                  <  2:8389208/93‑1 (MQ=255)
                                       tgaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttg                                                  >  2:2917540/1‑93 (MQ=255)
                                       tgaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttg                                                  >  2:1365239/1‑93 (MQ=255)
                                       tgaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttg                                                  <  1:8458281/93‑1 (MQ=255)
                                        gaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgt                                                 <  1:7157275/93‑1 (MQ=255)
                                        gaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgt                                                 <  2:1174150/93‑1 (MQ=255)
                                          atgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttg                                               >  2:6869694/1‑93 (MQ=255)
                                           tgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttgt                                              >  2:3026566/1‑93 (MQ=255)
                                           tgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttgt                                              >  1:5806348/1‑93 (MQ=255)
                                           tgaaGCCCAGTTCGCCCATTTCCTCGTAAGTCGCACCGTAGAAATCCAGGCCCAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttgt                                              <  1:78595/93‑1 (MQ=255)
                                            gaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTa                                             <  2:1892040/93‑1 (MQ=255)
                                            gaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTa                                             <  1:5581291/93‑1 (MQ=255)
                                               gCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                          >  2:5662307/1‑93 (MQ=255)
                                               gCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                          >  2:1413421/1‑93 (MQ=255)
                                               gCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                          >  1:5651016/1‑93 (MQ=255)
                                                cccAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGc                                         <  2:6141846/93‑1 (MQ=255)
                                                   aGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                      <  2:6746185/93‑1 (MQ=255)
                                                    gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  2:7159235/1‑93 (MQ=255)
                                                    gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  1:441923/1‑93 (MQ=255)
                                                    gTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                     >  1:9088253/1‑93 (MQ=255)
                                                     ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    >  1:248602/1‑93 (MQ=255)
                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   >  1:5183935/1‑93 (MQ=255)
                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  2:9083900/93‑1 (MQ=255)
                                                        gCCCATTTTCTCGTAAGTCGCGCCGTAGAAATCCAGGCACAGATGACGCAGCTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATa                                 <  1:9508934/93‑1 (MQ=255)
                                                            aTTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAccc                             >  2:743322/1‑93 (MQ=255)
                                                            aTTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAccc                             >  1:8078459/1‑93 (MQ=255)
                                                             ttttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCa                            >  2:10088/1‑93 (MQ=255)
                                                             ttttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCa                            <  1:3289249/93‑1 (MQ=255)
                                                             ttttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCa                            >  1:1050351/1‑93 (MQ=255)
                                                               ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATa                          <  2:5356913/93‑1 (MQ=255)
                                                               ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATa                          <  1:870949/93‑1 (MQ=255)
                                                                tctcGTAAGTCTCCCCGTAGAAATCCAGGCACAGATGACGCAACTCATCACAGATCTCCTGCGTGTTGTTGTAGTGCATCAGATAGCCCATAc                         <  2:6239000/93‑1 (MQ=255)
                                                                tctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATAc                         <  1:8260739/93‑1 (MQ=255)
                                                                 ctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACg                        <  2:4814753/93‑1 (MQ=255)
                                                                  tcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgg                       <  1:5265797/93‑1 (MQ=255)
                                                                  tcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgg                       <  2:6680049/93‑1 (MQ=255)
                                                                   cGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACggg                      <  1:3695979/93‑1 (MQ=255)
                                                                    gTAAGTCGCACCGTAGAAATCCACGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     <  2:5496184/93‑1 (MQ=255)
                                                                     tAAGTCGCACCGTAGCACTCCAGGCACAGATGACGCACCTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTg                    <  2:5924215/93‑1 (MQ=255)
                                                                      aaGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTgg                   <  1:3910022/93‑1 (MQ=255)
                                                                       aGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGc                  <  1:7610648/93‑1 (MQ=255)
                                                                       aGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGc                  <  2:7864689/93‑1 (MQ=255)
                                                                       aGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGc                  <  1:885475/93‑1 (MQ=255)
                                                                         tCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGa                <  1:2071660/93‑1 (MQ=255)
                                                                         tCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGa                <  2:3681241/93‑1 (MQ=255)
                                                                             cccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            <  1:9387694/92‑1 (MQ=255)
                                                                             aCCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTc            <  2:793074/93‑1 (MQ=255)
                                                                              ccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCa           <  2:8862847/93‑1 (MQ=255)
                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         >  1:2624781/1‑93 (MQ=255)
                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         <  1:982262/93‑1 (MQ=255)
                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         <  2:6535589/93‑1 (MQ=255)
                                                                                    aaaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatg     <  1:2903767/93‑1 (MQ=255)
                                                                                     aaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatga    >  2:8543315/1‑93 (MQ=255)
                                                                                     aaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatga    <  1:1276735/93‑1 (MQ=255)
                                                                                       tCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAtt  <  2:4575193/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506

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