Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,126,464 A→G V71V (GTA→GTG metB → cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,126,4640AG100.0% 811.6 / NA 248V71V (GTA→GTGmetBcystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (124/124);  total (124/124)

Only 100 of 265 total aligned reads displayed.

GCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTA  >  NC_000913/4126374‑4126553
                                                                                          |                                                                                         
 cgcgcgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtgctt                                                                                        >  2:3271722/1‑93 (MQ=255)
  gcgcgcATGTTTTCTCGGCTCGCGGCAACCCAACGCGCGATGTGGTTCAGTGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtgctta                                                                                       <  1:4442816/93‑1 (MQ=255)
  gcgcgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtgctta                                                                                       <  2:6001136/93‑1 (MQ=255)
   cgcgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtgcttac                                                                                      <  2:2517134/93‑1 (MQ=255)
   cgcgcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtgcttac                                                                                      >  1:6059019/1‑93 (MQ=255)
      gcATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTaa                                                                                   <  2:475360/93‑1 (MQ=255)
       cATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAAt                                                                                  >  1:1112601/1‑93 (MQ=255)
       cATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAAt                                                                                  >  1:3731046/1‑93 (MQ=255)
         tGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATAc                                                                                >  1:1287281/1‑93 (MQ=255)
         tGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATAc                                                                                <  2:4692677/93‑1 (MQ=255)
         tGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATAc                                                                                >  1:4756025/1‑93 (MQ=255)
           aTTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCg                                                                              <  1:4254146/93‑1 (MQ=255)
           aTTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCg                                                                              <  1:3327134/93‑1 (MQ=255)
             tACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGc                                                                            <  2:5781603/93‑1 (MQ=255)
             tACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGc                                                                            <  1:5032149/93‑1 (MQ=255)
               cTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCAt                                                                          <  2:920779/93‑1 (MQ=255)
                tcgcgtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATg                                                                         >  1:990105/1‑93 (MQ=255)
                  gcgtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTc                                                                       >  2:4597029/1‑93 (MQ=255)
                   cgtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTcc                                                                      <  1:2920429/93‑1 (MQ=255)
                   cgtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTcc                                                                      <  2:4758208/93‑1 (MQ=255)
                   cgtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTcc                                                                      <  1:2628293/93‑1 (MQ=255)
                    gtcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCcg                                                                     >  1:837275/1‑93 (MQ=255)
                     tcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCcgc                                                                    >  2:3428184/1‑93 (MQ=255)
                     tcgcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCcgc                                                                    >  1:1929198/1‑93 (MQ=255)
                     tcgcgGCAACACAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCcgg                                                                    >  1:1083332/1‑92 (MQ=255)
                       gcgGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGa                                                                  >  2:5828524/1‑93 (MQ=255)
                         ggCGACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGAtt                                                                <  2:195585/93‑1 (MQ=255)
                          gCACCCCAGCGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTc                                                               <  2:1816366/93‑1 (MQ=255)
                          gCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTc                                                               <  2:341161/93‑1 (MQ=255)
                          gCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTc                                                               >  1:1709604/1‑93 (MQ=255)
                          gCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTc                                                               >  2:1447387/1‑93 (MQ=255)
                          gCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTc                                                               <  2:4043903/93‑1 (MQ=255)
                            aaCCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGGAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCAc                                                             >  1:2440094/1‑93 (MQ=255)
                            aaCCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCAc                                                             <  1:3299054/93‑1 (MQ=255)
                            aaCCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCAc                                                             >  1:3178830/1‑93 (MQ=255)
                             aCCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCAcc                                                            <  2:3870244/93‑1 (MQ=255)
                              cccAACGCGCGACGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCt                                                           >  1:5234279/1‑93 (MQ=255)
                               ccAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTg                                                          <  1:72564/93‑1 (MQ=255)
                                cAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTgg                                                         >  1:461218/1‑93 (MQ=255)
                                 aaCGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGt                                                        >  2:521419/1‑93 (MQ=255)
                                 aaCGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGt                                                        <  1:5138829/93‑1 (MQ=255)
                                 aaCGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGt                                                        >  1:4476922/1‑93 (MQ=255)
                                   cgcgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTaa                                                      <  1:1385351/93‑1 (MQ=255)
                                    gcgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAAc                                                     <  2:5670106/93‑1 (MQ=255)
                                    gcgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAAc                                                     >  1:2011285/1‑93 (MQ=255)
                                     cgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACg                                                    <  1:2096660/93‑1 (MQ=255)
                                     cgcgATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACg                                                    <  2:2164168/93‑1 (MQ=255)
                                        gATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGAcc                                                 >  2:4185447/1‑93 (MQ=255)
                                        gATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGAcc                                                 >  2:6238469/1‑93 (MQ=255)
                                         aTGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCg                                                <  2:3381657/93‑1 (MQ=255)
                                            tgGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCt                                             >  1:2503700/1‑93 (MQ=255)
                                               ttCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCtttt                                          <  2:2973082/93‑1 (MQ=255)
                                                tCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCttttt                                         <  2:4875089/93‑1 (MQ=255)
                                                  aGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGa                                       <  1:1494080/93‑1 (MQ=255)
                                                   gCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGaa                                      >  1:3594052/1‑93 (MQ=255)
                                                   gCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGaa                                      <  1:3043082/93‑1 (MQ=255)
                                                   gCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGaa                                      <  1:3392843/93‑1 (MQ=255)
                                                    cGTGCGCTGGCAGAACTGGAAGGTGGTTCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGCCCGTCTTTTTGaaa                                     <  2:3232662/93‑1 (MQ=255)
                                                       gcgcTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCt                                  <  2:4382315/93‑1 (MQ=255)
                                                        cgcTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTg                                 <  2:2633539/93‑1 (MQ=255)
                                                         gcTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTgg                                <  1:4090251/93‑1 (MQ=255)
                                                           tGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCg                              >  2:2572853/1‑93 (MQ=255)
                                                           tGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCg                              <  2:1275133/93‑1 (MQ=255)
                                                           tGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCg                              >  2:2843736/1‑93 (MQ=255)
                                                           tGGCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCg                              >  2:5106916/1‑93 (MQ=255)
                                                            ggCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGa                             <  1:5329748/93‑1 (MQ=255)
                                                            ggCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGa                             <  2:2823918/93‑1 (MQ=255)
                                                            ggCAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGa                             >  2:748059/1‑93 (MQ=255)
                                                              cAGAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATc                           >  1:4706729/1‑93 (MQ=255)
                                                                gAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATctg                         >  2:5510634/1‑93 (MQ=255)
                                                                gAACTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATctg                         <  1:4800882/93‑1 (MQ=255)
                                                                 aaCTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATctgc                        <  1:5431506/93‑1 (MQ=255)
                                                                  aCTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATctgct                       <  1:2797371/93‑1 (MQ=255)
                                                                   cTGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATctgctg                      >  2:3769579/1‑93 (MQ=255)
                                                                    tGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTgg                     >  1:3986454/1‑93 (MQ=255)
                                                                    tGGAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTgg                     >  2:6128366/1‑93 (MQ=255)
                                                                     gtaaGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGCCTTTTTGAAACCTGGCGATCTGCTGGt                    >  2:1238430/3‑93 (MQ=255)
                                                                     ggAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACTTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGt                    <  2:2312646/93‑1 (MQ=255)
                                                                     ggAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGt                    >  2:2454076/1‑93 (MQ=255)
                                                                     ggAAGGGGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCCCCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGt                    <  2:1648622/93‑1 (MQ=255)
                                                                      gAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGtt                   >  1:4101238/1‑93 (MQ=255)
                                                                      gAAGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGtt                   >  2:3892543/1‑93 (MQ=255)
                                                                        aGGTGGTGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTgc                 <  2:534137/93‑1 (MQ=255)
                                                                         ggtggtGCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTgcg                >  2:213062/1‑93 (MQ=255)
                                                                          gtggtgCTGGTGCAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTgcgc               >  1:2578139/1‑93 (MQ=255)
                                                                           tggtgctggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgcc              <  1:3274438/93‑1 (MQ=255)
                                                                           tggtgctggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgcc              <  2:2437911/93‑1 (MQ=255)
                                                                           tggtgctggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgcc              <  1:5938490/93‑1 (MQ=255)
                                                                            ggtgctggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccg             >  1:4738734/1‑93 (MQ=255)
                                                                             gtgctggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccgc            <  1:2761235/93‑1 (MQ=255)
                                                                             gtgctggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccgc            <  2:6177681/93‑1 (MQ=255)
                                                                                 tggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGa        >  1:881533/1‑93 (MQ=255)
                                                                                  ggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGAc       <  1:220108/93‑1 (MQ=255)
                                                                                  ggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGAc       <  2:5011310/93‑1 (MQ=255)
                                                                                  ggtgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGAc       >  2:1178681/1‑93 (MQ=255)
                                                                                    tgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTg     >  1:1795153/1‑93 (MQ=255)
                                                                                    tgcAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTg     >  1:6278452/1‑93 (MQ=255)
                                                                                      cAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCt   >  2:3181010/1‑93 (MQ=255)
                                                                                      cAGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCt   >  1:4141971/1‑93 (MQ=255)
                                                                                       aGTGCTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTa  <  2:726356/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTA  >  NC_000913/4126374‑4126553

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: