Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,225,050 A→G E1215G (GAG→GGG)  metH → B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,0500AG99.5% 678.7 / ‑4.8 210E1215G (GAG→GGG) metHB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (123/86);  total (123/87)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 213 total aligned reads displayed.

CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4224960‑4225138
                                                                                          |                                                                                        
cGGACAGAAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTCCCGAAGTTGGgc                                                                                        <  1:3408585/93‑1 (MQ=255)
  gACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgcgc                                                                                      >  1:3088071/1‑93 (MQ=255)
      gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                  >  1:5697374/1‑93 (MQ=255)
      gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTACGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                  >  1:3265025/1‑93 (MQ=255)
       cAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggct                                                                                 >  1:3797475/1‑93 (MQ=255)
         agaactagGCTGTAGCAAAAATACAGCGGGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                               <  2:2542906/90‑1 (MQ=255)
         agaactacGGTGTAGCAAAAATTCAGCGCGATGAGGTTGAAGTTTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                               <  1:4286731/90‑1 (MQ=255)
         aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                               <  2:2750179/93‑1 (MQ=255)
          gTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCtggctggc                                                                              >  1:1546713/1‑93 (MQ=255)
          gTACTACCCTTTAGCACAAATTCAGCGCGACCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCtggctggc                                                                              <  1:2711048/93‑1 (MQ=255)
            actacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                            >  2:2117068/1‑93 (MQ=255)
            actacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                            >  2:4148491/1‑93 (MQ=255)
             ctgcGCTGTAGCACAAATTCAGCGCGATGAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAcc                                                                           <  1:3011706/90‑1 (MQ=255)
              tacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCg                                                                          >  1:4216101/1‑93 (MQ=255)
              tacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCg                                                                          >  2:3741291/1‑93 (MQ=255)
              tacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCg                                                                          >  2:5412690/1‑93 (MQ=255)
               ccGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGa                                                                         <  1:3582207/92‑1 (MQ=255)
               acGATGTCGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGa                                                                         <  1:412241/93‑1 (MQ=255)
                 gCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAAt                                                                       >  2:2784523/1‑93 (MQ=255)
                 gCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAAt                                                                       >  2:4484353/1‑93 (MQ=255)
                   tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                     >  2:4210855/1‑93 (MQ=255)
                   tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                     >  1:563407/1‑93 (MQ=255)
                   tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                     >  2:4420017/1‑93 (MQ=255)
                   tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                     >  1:5312616/1‑93 (MQ=255)
                   tGTAGCAAAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                     <  1:311184/93‑1 (MQ=255)
                    gTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                    >  1:3761018/1‑93 (MQ=255)
                     tAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgg                                                                   >  1:2229094/1‑93 (MQ=255)
                      agaacaAATTCAGCGCGATCAGTTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                  <  1:1855798/90‑1 (MQ=255)
                      accacaAATCCAGCGCGCTCAGGTTGAAAATTGTCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                  <  2:2447522/91‑1 (MQ=255)
                      aGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                  >  2:2204186/1‑93 (MQ=255)
                       gccaaaaTTCAGCAAGATAAGTTTGAAGTTTACGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                 <  1:4534119/89‑1 (MQ=255)
                       gCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                 >  1:725284/1‑93 (MQ=255)
                         aaaaaTTCGGTGCGAAAAGTTTAAAGATCTTGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                               <  2:3373338/91‑1 (MQ=255)
                         aaaaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                               <  2:6088307/91‑1 (MQ=255)
                           aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATg                                                             >  2:5134761/1‑93 (MQ=255)
                            atttCAGCGCGATCAGTTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                            <  2:4677479/91‑1 (MQ=255)
                            aaTTCANCTCGATCAGGTGAAATATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                            <  1:646829/93‑1 (MQ=255)
                            aaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                            >  1:5873667/1‑93 (MQ=255)
                             aTTCAGCGCGATAAGGTGGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAc                                                           <  2:5287286/93‑1 (MQ=255)
                              ttCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcg                                                          >  2:4946697/1‑93 (MQ=255)
                               tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                         >  1:5180688/1‑93 (MQ=255)
                                cAGCGCGAACAGGTTAAAGATTCTGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgcg                                                        <  1:6108106/93‑1 (MQ=255)
                                 aGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCgg                                                       >  1:109421/1‑93 (MQ=255)
                                  gcgcGATCAGATTGAAGATAATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGa                                                      <  1:5926157/93‑1 (MQ=255)
                                   cgcgATCAGTTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGAc                                                     <  1:89320/93‑1 (MQ=255)
                                   cgcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGAc                                                     <  1:2335505/93‑1 (MQ=255)
                                      gATGAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                  <  2:576556/93‑1 (MQ=255)
                                      gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                  >  2:3952989/1‑93 (MQ=255)
                                       acaaGGTTAAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 <  2:1163638/90‑1 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 >  2:1059288/1‑93 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                 >  2:4065126/1‑93 (MQ=255)
                                         cAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTc                                               >  2:2617200/1‑93 (MQ=255)
                                           ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcac                                             >  1:3937721/1‑93 (MQ=255)
                                           ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcac                                             >  2:4553188/1‑93 (MQ=255)
                                               gAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAAt                                         >  2:706925/1‑93 (MQ=255)
                                               gAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGCATGACGCGGACTGATTCACAAAt                                         >  2:822441/1‑93 (MQ=255)
                                                aaGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATc                                        <  1:1718705/93‑1 (MQ=255)
                                                   aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGt                                     >  2:3693709/1‑93 (MQ=255)
                                                   aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGt                                     >  2:3524736/1‑93 (MQ=255)
                                                    ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                    >  2:1433369/1‑93 (MQ=255)
                                                      aTGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCAc                                  <  2:2969257/93‑1 (MQ=255)
                                                       tGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACt                                 >  2:567679/1‑93 (MQ=255)
                                                       tGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACt                                 <  1:4805188/93‑1 (MQ=255)
                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                              >  1:2781442/1‑93 (MQ=255)
                                                            gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  2:2908166/93‑1 (MQ=255)
                                                            gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                            <  2:4982296/93‑1 (MQ=255)
                                                              cgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCtt                          >  1:3365657/1‑93 (MQ=255)
                                                               gTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                         <  2:5626905/93‑1 (MQ=255)
                                                                tAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTac                        <  2:157385/93‑1 (MQ=255)
                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                  >  1:22835/1‑93 (MQ=255)
                                                                       aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                 >  2:5838725/1‑93 (MQ=255)
                                                                       aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                 <  1:4792585/93‑1 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                >  1:2184059/1‑93 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                >  1:3020622/1‑93 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaaaaaaca                >  2:6016266/1‑93 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg               >  1:2156281/1‑93 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg               <  1:2838250/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg              >  1:3452107/1‑93 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg              <  2:4905050/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg              <  2:4050660/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg              >  2:5230354/1‑93 (MQ=255)
                                                                           gCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAggg             >  2:629106/1‑93 (MQ=255)
                                                                             gTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcg           >  1:5228404/1‑93 (MQ=255)
                                                                             gTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcg           >  2:5528054/1‑93 (MQ=255)
                                                                                aCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc        >  2:611389/1‑93 (MQ=255)
                                                                                 ccGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACACACAGGGCGCTCa       >  1:2491520/1‑93 (MQ=255)
                                                                                 ccGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGTGCGGTCa       >  2:3462284/1‑93 (MQ=255)
                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa      <  2:6176525/93‑1 (MQ=255)
                                                                                   gAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt     <  1:1027946/93‑1 (MQ=255)
                                                                                   gAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt     <  2:264027/93‑1 (MQ=255)
                                                                                   gAAGTTGGGCGCTGGCTGGAACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt     >  2:6162741/1‑93 (MQ=255)
                                                                                    aaGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg    <  2:592995/93‑1 (MQ=255)
                                                                                    aaGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg    <  2:1325246/93‑1 (MQ=255)
                                                                                    aaGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg    <  2:1846933/93‑1 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg  >  2:2153027/1‑92 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa   >  1:5284208/1‑93 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa   >  2:619255/1‑93 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  <  2:5056577/93‑1 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg  >  1:3085056/1‑93 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGCGCGCTCAATGAg  >  2:5013789/1‑93 (MQ=255)
                                                                                          |                                                                                        
CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAG  >  NC_000913/4224960‑4225138

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