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.4% 1036.1 / ‑5.4 329E1215G (GAG→GGG) metHB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (173/154);  total (174/155)
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 337 total aligned reads displayed.

CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGT  >  NC_000913/4224960‑4225139
                                                                                          |                                                                                         
cGGACAGCAAGTATTACGTTGAGGCACAACTCCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgc                                                                                         <  1:2467142/93‑1 (MQ=255)
 ggACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGgcg                                                                                        >  1:3488751/1‑93 (MQ=255)
   aCAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                      >  1:4346677/1‑93 (MQ=255)
   aCAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgct                                                                                      >  1:7705057/1‑93 (MQ=255)
    cAGCAAGTCCTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctg                                                                                     <  1:36758/93‑1 (MQ=255)
    cAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctg                                                                                     >  1:2133825/1‑93 (MQ=255)
    cAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctg                                                                                     >  2:2081959/1‑93 (MQ=255)
      gCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggc                                                                                   >  1:836142/1‑93 (MQ=255)
       cAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggct                                                                                  >  1:7479128/1‑93 (MQ=255)
        aaGTACTACGCTGTAGCACAAATTCGGCGCGATCAGGTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCgctggctg                                                                                 <  2:3870303/93‑1 (MQ=255)
         aGTTCTACGGTGTAGAACAAATTCGGCGGAAGCAGGTAGAAAATTAGCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGctggctgg                                                                                <  1:535007/93‑1 (MQ=255)
           tactacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCa                                                                              >  1:7283650/1‑93 (MQ=255)
           tactacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCa                                                                              >  2:5678353/1‑93 (MQ=255)
            gctacGCTGAAGCAAAAATCCAGCGCGACGGGGTTGAAGATTCCGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                             <  2:317954/92‑1 (MQ=255)
            actacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCAc                                                                             >  2:1426111/1‑93 (MQ=255)
              tacGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCg                                                                           >  1:7731164/1‑93 (MQ=255)
               acGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGa                                                                          >  2:4762158/1‑93 (MQ=255)
                  cTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                       >  1:6942229/1‑93 (MQ=255)
                  cTGTAGCACAAATCCAGCGCGATCAGGTGAAAGATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATc                                                                       <  1:100449/93‑1 (MQ=255)
                   tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                      <  2:1626811/93‑1 (MQ=255)
                   tGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCt                                                                      <  2:2342523/93‑1 (MQ=255)
                    gTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTg                                                                     >  1:7887706/1‑93 (MQ=255)
                     tAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgg                                                                    <  2:3352083/93‑1 (MQ=255)
                       gCAAAAATTCAGCGCGATCAGGTTGAAGTTTACGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTgggg                                                                  <  1:6338046/93‑1 (MQ=255)
                        caaaaaTCCAGCGCGATCAGGTTGAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGt                                                                 <  2:6589410/90‑1 (MQ=255)
                        caaaaaTCCAGCGCGATCAGGTGAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGt                                                                 <  2:4507375/90‑1 (MQ=255)
                         acaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                                <  1:2731128/93‑1 (MQ=255)
                         acaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                                >  1:6510009/1‑93 (MQ=255)
                         acaAATTCAGCGCGATCAGGTTAAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                                <  1:41504/93‑1 (MQ=255)
                         aaaaaTTCAGCGCGATCAGTTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                                <  1:2757289/91‑1 (MQ=255)
                         aaaaaTTCAGCGCGATCAGGTTGAAGATTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTa                                                                <  2:5945339/91‑1 (MQ=255)
                          aaaaTTCGGCGCAATCAGTTTGAAGATGATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTAt                                                               <  2:4148072/92‑1 (MQ=255)
                            aaTTCAGCGCGGTCAGGTTGAAAATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                             <  1:411725/93‑1 (MQ=255)
                            aaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGa                                                             >  1:3044916/1‑93 (MQ=255)
                              ttCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcg                                                           >  1:5776513/1‑93 (MQ=255)
                               tCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                          >  2:1399070/1‑93 (MQ=255)
                               tCAGCGCGATCAGGTGGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGAcgc                                                          <  1:2341877/93‑1 (MQ=255)
                                  gcgcGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGa                                                       >  2:2663703/1‑93 (MQ=255)
                                  gcgcGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGa                                                       >  2:7109152/1‑93 (MQ=255)
                                     ggATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTg                                                    <  2:6852829/92‑1 (MQ=255)
                                     cgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGCATCTGGGGTATGACGCGGAGTg                                                    >  2:1193467/1‑93 (MQ=255)
                                       aTCGGGTTGAAGTTTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                  <  1:1575391/93‑1 (MQ=255)
                                       aTCAGTTTGAAGATTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGAt                                                  <  1:7427691/93‑1 (MQ=255)
                                       aTCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGAt                                                  >  2:4251695/1‑93 (MQ=255)
                                         gaGGTTGAAGTTTAGGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTc                                                <  1:7764639/92‑1 (MQ=255)
                                         cAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTc                                                >  1:4548714/1‑93 (MQ=255)
                                          aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTca                                               >  1:7145163/1‑93 (MQ=255)
                                            gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTccca                                             >  2:1678825/1‑93 (MQ=255)
                                            gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCGGGCACCGAATCTGGGGCATGACGCGGACTGATTcaca                                             >  1:876811/1‑93 (MQ=255)
                                              tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACaaa                                           <  2:3588713/93‑1 (MQ=255)
                                               gAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAAt                                          >  2:4093933/1‑93 (MQ=255)
                                                aaGATTCTCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATc                                         <  1:3939184/93‑1 (MQ=255)
                                                  gATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTg                                       >  2:274989/1‑93 (MQ=255)
                                                  gATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTg                                       <  2:6866606/93‑1 (MQ=255)
                                                  atttATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTg                                       <  1:5734008/91‑1 (MQ=255)
                                                   aTTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGt                                      >  2:5612769/1‑93 (MQ=255)
                                                   aTTATCCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGt                                      <  2:3097935/93‑1 (MQ=255)
                                                    ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTc                                     >  2:5378136/1‑93 (MQ=255)
                                                     tATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCa                                    >  1:5566789/1‑93 (MQ=255)
                                                     gaaacccGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATACACAAATCTGTCa                                    <  2:6987659/89‑1 (MQ=255)
                                                       tGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACt                                  >  2:1834561/1‑93 (MQ=255)
                                                         cccGCCGTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACttt                                <  1:1960768/93‑1 (MQ=255)
                                                           cgccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTc                              >  2:5848349/1‑93 (MQ=255)
                                                            gccgTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTcc                             >  1:1251364/1‑93 (MQ=255)
                                                               gTAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTa                          <  1:2534037/93‑1 (MQ=255)
                                                                tAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTac                         >  1:6843718/1‑93 (MQ=255)
                                                                tAAAGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACAGATTCACAAATCTGTCACTTTTCCTTac                         >  2:2416295/1‑93 (MQ=255)
                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                        >  1:1177570/1‑93 (MQ=255)
                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                        >  2:6070862/1‑93 (MQ=255)
                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                        <  1:3845077/93‑1 (MQ=255)
                                                                 aaaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTaca                        >  1:6249544/1‑93 (MQ=255)
                                                                  aaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaa                       <  2:5030362/93‑1 (MQ=255)
                                                                  aaGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaa                       <  2:5189729/93‑1 (MQ=255)
                                                                   aGGTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTacaac                      <  2:61311/93‑1 (MQ=255)
                                                                     gTATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTAcaacaa                    >  2:101867/1‑93 (MQ=255)
                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   >  1:7217004/1‑93 (MQ=255)
                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   <  1:5148794/93‑1 (MQ=255)
                                                                      tATGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaa                   >  2:527540/1‑93 (MQ=255)
                                                                       aTGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaac                  >  1:7849219/1‑93 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                 <  2:7705056/93‑1 (MQ=255)
                                                                        tGAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACaacaaaca                 <  2:4346676/93‑1 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg                <  2:2529710/93‑1 (MQ=255)
                                                                         gAGCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAg                <  2:6805046/93‑1 (MQ=255)
                                                                           gCGTTACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAggg              >  1:6189701/1‑93 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  2:3797622/1‑93 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  2:5063180/1‑93 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           <  2:7211466/93‑1 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  1:2236219/1‑93 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           >  1:3292336/1‑93 (MQ=255)
                                                                              ttACCGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGgcgc           <  2:225846/93‑1 (MQ=255)
                                                                                  cGAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa       <  1:2848780/93‑1 (MQ=255)
                                                                                   gAAGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt      <  1:1400414/93‑1 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    <  1:7392655/93‑1 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    <  1:4665659/93‑1 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    >  1:4840643/1‑93 (MQ=255)
                                                                                     aGTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa    >  2:4818418/1‑93 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACATCAAACAGGGCGCTCAATGAg   <  2:7497273/93‑1 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg   >  2:4023016/1‑93 (MQ=255)
                                                                                      gTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAg   <  1:4954451/93‑1 (MQ=255)
                                                                                       ttGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgc  <  2:7118433/93‑2 (MQ=255)
                                                                                          |                                                                                         
CGGACAGCAAGTACTACGCTGTAGCACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGT  >  NC_000913/4224960‑4225139

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