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.7% 1178.3 / ‑6.5 358E1215G (GAG→GGG) metHB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (203/154);  minor base C (1/0);  total (204/154)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 375 total aligned reads displayed.

CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGT  >  NC_000913/4224960‑4225139
                                                                                          |                                                                                         
  gACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgcgc                                                                                       >  2:1720911/1‑93 (MQ=255)
   acaacaAGTACTCCGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                      <  2:7411969/93‑1 (MQ=255)
    caacAAGTACTACGCTGAAGCAAAAATTCAGCGCGATCGGCTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctg                                                                                     <  2:9743905/90‑1 (MQ=255)
    cAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGTgctg                                                                                     >  1:1530304/1‑93 (MQ=255)
    cAGAAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGTTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctg                                                                                     <  2:3714635/93‑1 (MQ=255)
     aGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctgg                                                                                    >  1:2423154/1‑93 (MQ=255)
     aGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctgg                                                                                    <  2:8698060/93‑1 (MQ=255)
      gaaaGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                   <  2:8884125/91‑1 (MQ=255)
       caggTACTACGCTGAAGCACAAATGCAGCGCGATCGGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggct                                                                                  <  1:6930343/90‑1 (MQ=255)
       cAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggct                                                                                  >  1:656282/1‑93 (MQ=255)
         aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                                <  2:7505405/93‑1 (MQ=255)
         aGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGCGCGctggctgg                                                                                >  1:5904029/1‑93 (MQ=255)
           tactacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGCATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCa                                                                              >  2:8295806/1‑93 (MQ=255)
            attacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                             <  1:8724914/91‑1 (MQ=255)
            actccGCTGTAGCACAAATTCAGCGCGACCAGGTTGAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                             <  1:4767456/93‑1 (MQ=255)
             ctacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAcc                                                                            <  1:5590317/93‑1 (MQ=255)
             ctacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATACCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAcc                                                                            >  2:4397203/1‑93 (MQ=255)
                  cTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                       >  2:2557841/1‑93 (MQ=255)
                  cTGTAGCAAAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                       <  1:3887905/93‑1 (MQ=255)
                    gTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                     <  1:4608604/93‑1 (MQ=255)
                    gTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                     >  1:3498969/1‑93 (MQ=255)
                    gTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                     >  2:9468851/1‑93 (MQ=255)
                    gTAGCACAAATTCAGCGCGATCAGGTGGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                     <  2:7594987/93‑1 (MQ=255)
                      aGCACAAATTCAGCGCGATCAGGTGGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTggg                                                                   <  1:3664315/93‑1 (MQ=255)
                       gCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                  >  2:1543563/1‑93 (MQ=255)
                       gCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                  >  2:5392415/1‑93 (MQ=255)
                       gCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGACCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                  >  2:5262707/1‑93 (MQ=255)
                        cacaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGt                                                                 >  1:169871/1‑93 (MQ=255)
                          caAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTAt                                                               >  1:5207349/1‑93 (MQ=255)
                          aaaaTTCGGGGCGATCGGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTAt                                                               <  1:8424000/92‑1 (MQ=255)
                           aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATg                                                              <  2:2935913/93‑1 (MQ=255)
                           aaaTCCGGCGCGCTCAGGTTAAATTTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATg                                                              <  1:9658968/93‑1 (MQ=255)
                            atttCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                             <  1:2190531/91‑1 (MQ=255)
                            aaTTCAGCGCGATCAGTTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                             <  1:5444860/93‑1 (MQ=255)
                            aaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGa                                                             >  1:4415261/1‑93 (MQ=255)
                             aTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAc                                                            <  1:7957997/93‑1 (MQ=255)
                              ttCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcg                                                           >  2:7214826/1‑93 (MQ=255)
                               tCAGCGCGATCAGGTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                          <  2:6703137/93‑1 (MQ=255)
                                  gcgcGAGCAGTTTGAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGa                                                       <  1:4827868/93‑1 (MQ=255)
                                      gATCAGTTAAAAGATTACGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                   <  2:3983632/93‑1 (MQ=255)
                                      gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGa                                                   >  2:1327583/1‑93 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                  >  1:8010667/1‑93 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                  >  2:6463792/1‑93 (MQ=255)
                                        tCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAtt                                                 >  1:9761826/1‑93 (MQ=255)
                                         cggtttGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTc                                                <  2:5581741/89‑1 (MQ=255)
                                          agtttGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                               <  1:1416144/90‑1 (MQ=255)
                                          aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                               >  2:2836624/1‑93 (MQ=255)
                                           ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcac                                              >  2:4985099/1‑93 (MQ=255)
                                            gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcaca                                             >  1:8215510/1‑93 (MQ=255)
                                            gTTGAAGATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTcaca                                             <  1:6689426/93‑1 (MQ=255)
                                              tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaaa                                           <  2:8014506/93‑1 (MQ=255)
                                               gAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGTCTGATTCACAAAt                                          >  2:5098501/1‑93 (MQ=255)
                                               aaaGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAAt                                          <  1:7319187/92‑1 (MQ=255)
                                                aaGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATc                                         <  1:7696346/93‑1 (MQ=255)
                                                aaGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATc                                         <  2:9885665/93‑1 (MQ=255)
                                                 aGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCCTAAATCt                                        >  2:4206242/1‑93 (MQ=255)
                                                 aGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCt                                        <  1:8352939/93‑1 (MQ=255)
                                                  aaTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTg                                       <  2:8930140/92‑1 (MQ=255)
                                                     tATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCa                                    >  2:5886139/1‑93 (MQ=255)
                                                      tggCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCAc                                   <  1:4299857/91‑1 (MQ=255)
                                                        gCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCCCtt                                 >  1:4963248/1‑93 (MQ=255)
                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                               <  1:9982152/93‑1 (MQ=255)
                                                          ccgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACtttt                               >  2:3364391/1‑93 (MQ=255)
                                                            gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             >  1:4103329/1‑93 (MQ=255)
                                                            gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             >  1:3150774/1‑93 (MQ=255)
                                                             ccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCt                            <  1:7819847/93‑1 (MQ=255)
                                                             ccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCt                            >  2:6043446/1‑93 (MQ=255)
                                                              cgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCtt                           <  1:5730070/93‑1 (MQ=255)
                                                               gTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                          >  1:7051629/1‑93 (MQ=255)
                                                                tAAAGGTATGCGCGTTACCGACGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTac                         >  1:2007173/1‑93 (MQ=255)
                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                        >  1:8887548/1‑93 (MQ=255)
                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                        <  1:619719/93‑1 (MQ=255)
                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaaaaaaaa                   >  2:1067668/1‑93 (MQ=255)
                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   >  2:5415148/1‑93 (MQ=255)
                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   <  2:2076975/93‑1 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                 >  1:9622301/1‑93 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                 >  2:9293672/1‑93 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg                >  2:5768276/1‑93 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg               <  1:2509173/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg               <  2:6890359/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg               <  2:6923245/93‑1 (MQ=255)
                                                                          aGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAgg               <  1:6820326/93‑1 (MQ=255)
                                                                            cGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgc             <  2:8144693/93‑1 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  2:8431833/1‑93 (MQ=255)
                                                                                aCCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc         >  2:9532154/1‑93 (MQ=255)
                                                                                 ccGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCa        >  1:1183277/1‑93 (MQ=255)
                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa       <  2:7095935/93‑1 (MQ=255)
                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa       >  1:8564220/1‑93 (MQ=255)
                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa       >  1:1021560/1‑93 (MQ=255)
                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa       <  2:2618823/93‑1 (MQ=255)
                                                                                   gAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt      >  2:3512721/1‑93 (MQ=255)
                                                                                    aaGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg     <  1:9232753/93‑1 (MQ=255)
                                                                                    aaGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATg     >  1:2937224/1‑93 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    >  2:4525428/1‑93 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg   >  1:6521297/1‑93 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg   <  1:2646949/93‑1 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCGTACAACAAACAGGGCGGTAATTGAg   <  2:3611515/93‑1 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGAGTATGACGCAGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAACGAg   <  2:4647911/93‑1 (MQ=255)
                                                                                       ttGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgg  >  1:8908324/1‑92 (MQ=255)
                                                                                       ttGGGCGCTGGCTGGCACCGAATCGGGGGTATCACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgc  <  2:9032148/93‑2 (MQ=255)
                                                                                          |                                                                                         
CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGT  >  NC_000913/4224960‑4225139

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