Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,224,890 G→A G1162S (GGC→AGC)  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,224,8900GA84.2% 732.6 / 104.7 265G1162S (GGC→AGC) metHB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase
Reads supporting (aligned to +/- strand):  ref base G (16/26);  new base A (126/97);  total (142/123)
Fisher's exact test for biased strand distribution p-value = 4.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 270 total aligned reads displayed.

AGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTA  >  NC_000913/4224798‑4224982
                                                                                            |                                                                                            
  ggCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGc                                                                                            <  2:6073704/93‑1 (MQ=255)
  ggCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGc                                                                                            <  2:73939/93‑1 (MQ=255)
   gCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCa                                                                                           <  1:4021826/93‑1 (MQ=255)
   gCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCa                                                                                           <  2:3541738/93‑1 (MQ=255)
      tccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGa                                                                                        <  2:4604250/93‑1 (MQ=255)
      tccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGa                                                                                        >  2:5130321/1‑93 (MQ=255)
      tccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGa                                                                                        >  2:5418878/1‑93 (MQ=255)
      tccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGa                                                                                        >  1:3100215/1‑93 (MQ=255)
      tccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGa                                                                                        >  2:261783/1‑93 (MQ=255)
      tccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGa                                                                                        >  2:508971/1‑93 (MQ=255)
      tccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGa                                                                                        >  1:4287468/1‑93 (MQ=255)
       ccgtccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGaa                                                                                       <  1:1847065/93‑1 (MQ=255)
          tccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACt                                                                                    <  1:5129218/93‑1 (MQ=255)
          tccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACt                                                                                    >  1:718015/1‑93 (MQ=255)
          tccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCCGCTGGAAGTGGAAAAACACACTGGCATGAAACt                                                                                    >  1:6114331/1‑93 (MQ=255)
           ccgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTc                                                                                   <  1:4086217/93‑1 (MQ=255)
            cgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTca                                                                                  <  2:1177029/93‑1 (MQ=255)
            cgGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTca                                                                                  >  2:2455884/1‑93 (MQ=255)
             ggCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGCAGTGGAAAAACACACTAGCATGAAACTcac                                                                                 >  2:4913054/1‑93 (MQ=255)
             ggCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTcac                                                                                 <  1:2673409/93‑1 (MQ=255)
             ggCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTcac                                                                                 >  2:2506911/1‑93 (MQ=255)
              gCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTcaca                                                                                >  1:3123523/1‑93 (MQ=255)
                aCCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGa                                                                              <  2:6170075/93‑1 (MQ=255)
                 ccGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGaa                                                                             <  1:4466800/93‑1 (MQ=255)
                  cGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAAt                                                                            <  1:5720721/93‑1 (MQ=255)
                   gggCTATCCGGCCTGCCCGGAACATACGGAAAAATCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATc                                                                           >  1:1564181/1‑93 (MQ=255)
                   gggCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATc                                                                           <  1:1137486/93‑1 (MQ=255)
                   gggCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATc                                                                           <  1:351671/93‑1 (MQ=255)
                    ggCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCt                                                                          >  1:1821708/1‑93 (MQ=255)
                     gCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCtt                                                                         <  1:2582118/93‑1 (MQ=255)
                     gCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCtt                                                                         >  2:442813/1‑93 (MQ=255)
                      cTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCttt                                                                        <  1:1247789/93‑1 (MQ=255)
                      cTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCttt                                                                        <  1:5743223/93‑1 (MQ=255)
                       tATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTc                                                                       <  1:255008/93‑1 (MQ=255)
                        aTCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCg                                                                      >  1:2900015/1‑93 (MQ=255)
                          ccGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGcc                                                                    <  1:3855085/93‑1 (MQ=255)
                          ccGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGcc                                                                    <  1:6218958/93‑1 (MQ=255)
                            ggCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCAt                                                                  >  1:3302095/1‑93 (MQ=255)
                             gCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCAtg                                                                 >  2:3001701/1‑93 (MQ=255)
                               cTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCAtgtg                                                               >  1:616392/1‑93 (MQ=255)
                                tGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTgg                                                              >  2:25251/1‑93 (MQ=255)
                                tGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTgg                                                              >  2:3384104/1‑93 (MQ=255)
                                   ccGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGccc                                                           >  1:3775835/1‑93 (MQ=255)
                                   ccGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGccc                                                           >  1:2522303/1‑93 (MQ=255)
                                    cGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCg                                                          >  1:1301524/1‑93 (MQ=255)
                                    cGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCg                                                          >  2:109210/1‑93 (MQ=255)
                                     ggAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCgg                                                         >  2:3785813/1‑93 (MQ=255)
                                      gAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGt                                                        >  1:883523/1‑93 (MQ=255)
                                       aaCATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTg                                                       <  1:2041009/93‑1 (MQ=255)
                                       aaCATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTg                                                       >  1:3221065/1‑93 (MQ=255)
                                       aaCATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTg                                                       >  1:278595/1‑93 (MQ=255)
                                        aCATACGGAAAAAGCCACCATCTGAGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGc                                                      <  2:4995144/93‑1 (MQ=255)
                                         cATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCa                                                     <  1:4799401/93‑1 (MQ=255)
                                         cATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCa                                                     <  2:1426661/93‑1 (MQ=255)
                                           tACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATc                                                   <  1:4952635/93‑1 (MQ=255)
                                           tACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATc                                                   >  1:4283130/1‑93 (MQ=255)
                                           tACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATc                                                   <  1:785775/93‑1 (MQ=255)
                                           tACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATc                                                   >  2:5723648/1‑93 (MQ=255)
                                             cGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCgg                                                 >  2:6400985/1‑93 (MQ=255)
                                             cGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCgg                                                 >  2:4278506/1‑93 (MQ=255)
                                              ggAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGt                                                <  1:1889030/93‑1 (MQ=255)
                                               gAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGtt                                               <  2:1877231/93‑1 (MQ=255)
                                               gAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGtt                                               >  1:1390668/1‑93 (MQ=255)
                                               gAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGtt                                               <  1:5583058/93‑1 (MQ=255)
                                                 aaaaGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTc                                             >  1:2479510/1‑93 (MQ=255)
                                                    aGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCggg                                          >  2:5449754/1‑93 (MQ=255)
                                                     gCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGt                                         >  2:4766690/1‑93 (MQ=255)
                                                     gCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGACCGGTGCATCGGTTTCGGGt                                         >  1:5656324/1‑93 (MQ=255)
                                                     accaccaTCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGt                                         <  2:5593254/92‑1 (MQ=255)
                                                      ccaccaTCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGtt                                        >  2:867161/1‑93 (MQ=255)
                                                       caccaTCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTg                                       >  1:959/1‑93 (MQ=255)
                                                         ccaTCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGt                                     <  2:3099037/93‑1 (MQ=255)
                                                           aTCTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTAc                                   <  2:6447985/93‑1 (MQ=255)
                                                             cTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACtt                                 <  2:6514881/93‑1 (MQ=255)
                                                             cTGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACtt                                 >  1:4271660/1‑93 (MQ=255)
                                                              tGGGAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTc                                >  2:3493934/1‑93 (MQ=255)
                                                               gggAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCa                               >  2:1611202/1‑93 (MQ=255)
                                                                ggAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAg                              >  2:2160016/1‑93 (MQ=255)
                                                                ggAGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAg                              >  2:2096110/1‑93 (MQ=255)
                                                                  aGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGcc                            >  2:5976467/1‑93 (MQ=255)
                                                                  aGCTGCTGGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGcc                            <  1:602409/93‑1 (MQ=255)
                                                                    ctgctgGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCAc                          <  1:2844039/93‑1 (MQ=255)
                                                                    ctgctgGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCAc                          <  2:3329165/93‑1 (MQ=255)
                                                                        tgGAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCgg                      >  1:3622698/1‑93 (MQ=255)
                                                                         ggAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGa                     <  2:5207876/93‑1 (MQ=255)
                                                                          gAAGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGAc                    <  2:2858581/93‑1 (MQ=255)
                                                                           aaGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACa                   >  1:2369470/1‑93 (MQ=255)
                                                                            aGTGGAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAg                  <  1:6343241/93‑1 (MQ=255)
                                                                                gAAAAACACACTAGCATGAAACTCACAGAATGTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAg              >  1:5308358/1‑93 (MQ=255)
                                                                                gAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAg              >  1:5384822/1‑93 (MQ=255)
                                                                                gAAAAACACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAg              <  2:2268454/93‑1 (MQ=255)
                                                                                 aaaaaCACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGt             >  1:4645016/1‑93 (MQ=255)
                                                                                   aaaCACACTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGtac           >  2:3649175/1‑93 (MQ=255)
                                                                                     acacacTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGtacta         >  2:4022290/1‑93 (MQ=255)
                                                                                      cacacTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCCAGtactac        <  1:6493555/93‑1 (MQ=255)
                                                                                      cacacTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGtactac        >  1:108074/1‑93 (MQ=255)
                                                                                       acacTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACg       >  1:5027971/1‑93 (MQ=255)
                                                                                       acacTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACg       >  1:5474815/1‑93 (MQ=255)
                                                                                       acacTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACg       >  2:1954370/1‑93 (MQ=255)
                                                                                        cacTAGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGc      <  1:2266674/93‑1 (MQ=255)
                                                                                            |                                                                                            
AGGGCATCCGTCCGGCACCGGGCTATCCGGCCTGCCCGGAACATACGGAAAAAGCCACCATCTGGGAGCTGCTGGAAGTGGAAAAACACACTGGCATGAAACTCACAGAATCTTTCGCCATGTGGCCCGGTGCATCGGTTTCGGGTTGGTACTTCAGCCACCCGGACAGCAAGTACTACGCTGTA  >  NC_000913/4224798‑4224982

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: