Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,153,184 A→G K317E (AAA→GAA)  metC → cystathionine beta‑lyase/L‑cysteine desulfhydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,153,1840AG100.0% 322.7 / NA 86K317E (AAA→GAA) metCcystathionine beta‑lyase/L‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/53);  total (33/53)

TGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGAAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATCGCCGCCATTCGCCCA  >  NC_000913/3153047‑3153315
                                                                                                                                         |                                                                                                                                   
tGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGgaa                                                                                                                                   <  2:794094/140‑1 (MQ=255)
  aaaGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGgaaaa                                                                                                                                 <  1:359485/140‑1 (MQ=255)
   aaGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGgaaaaa                                                                                                                                <  1:555265/140‑1 (MQ=255)
   aaGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGgaaaaa                                                                                                                                <  2:648460/140‑1 (MQ=255)
     gTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACt                                                                                                                              <  1:525764/140‑1 (MQ=255)
     gTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACt                                                                                                                              <  2:385724/140‑1 (MQ=255)
       ggCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCa                                                                                                                            <  1:1395282/140‑1 (MQ=255)
        gCTGAATGGCTGGCAGAACATCCGCAAGTTGTGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCaa                                                                                                                           >  1:737465/1‑140 (MQ=255)
        gCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCaa                                                                                                                           >  2:18362/1‑140 (MQ=255)
        gCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCaa                                                                                                                           >  2:338105/1‑140 (MQ=255)
        gCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCaa                                                                                                                           >  1:1122748/1‑140 (MQ=255)
        gCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCaa                                                                                                                           >  1:867129/1‑140 (MQ=255)
        gCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCaa                                                                                                                           >  1:629108/1‑140 (MQ=255)
        gCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCaa                                                                                                                           >  2:972416/1‑140 (MQ=255)
                gctggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAag                                                                                                                   <  1:1060829/140‑1 (MQ=255)
                  tggcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAagag                                                                                                                 >  1:767073/1‑140 (MQ=255)
                  gggcAGAACGTCCGCAAGTTGCTCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTGCCTTTGTGCTTAAGGAAAAACTCAATAATGAagag                                                                                                                 <  1:1223635/139‑1 (MQ=255)
                    gcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCt                                                                                                               <  1:131247/140‑1 (MQ=255)
                    gcAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCt                                                                                                               <  1:575657/140‑1 (MQ=255)
                       gAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGc                                                                                                            <  1:926380/140‑1 (MQ=255)
                        aaCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCg                                                                                                           <  2:1341079/140‑1 (MQ=255)
                         aCATCCGCAAGTTGTGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGa                                                                                                          <  2:737465/140‑1 (MQ=255)
                         aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGa                                                                                                          <  1:626994/140‑1 (MQ=255)
                         aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGa                                                                                                          <  1:994490/140‑1 (MQ=255)
                         aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGa                                                                                                          <  1:225871/140‑1 (MQ=255)
                         aCATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGa                                                                                                          <  2:778389/140‑1 (MQ=255)
                          cATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGaa                                                                                                         >  2:1073640/1‑140 (MQ=255)
                          cATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGaa                                                                                                         <  2:1197173/140‑1 (MQ=255)
                          cATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGaa                                                                                                         <  1:182411/140‑1 (MQ=255)
                          cATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGaa                                                                                                         <  1:1436440/140‑1 (MQ=255)
                          cATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGaa                                                                                                         <  1:518306/140‑1 (MQ=255)
                          cATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGaa                                                                                                         <  1:713748/140‑1 (MQ=255)
                          cATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGCCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGaa                                                                                                         >  2:815680/1‑140 (MQ=255)
                            tCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACt                                                                                                       <  2:106191/140‑1 (MQ=255)
                            tCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACt                                                                                                       <  2:659192/140‑1 (MQ=255)
                              cGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTAt                                                                                                     >  1:311164/1‑140 (MQ=255)
                              cGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTAt                                                                                                     >  2:378218/1‑140 (MQ=255)
                                 aaGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTg                                                                                                  <  2:1246282/140‑1 (MQ=255)
                                  aGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTgg                                                                                                 >  1:1051203/1‑140 (MQ=255)
                                   gTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGa                                                                                                >  1:375277/1‑140 (MQ=255)
                                    ttGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGAt                                                                                               <  1:381634/140‑1 (MQ=255)
                                       cgcgAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAAc                                                                                            >  1:336993/1‑140 (MQ=255)
                                        gcgAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACt                                                                                           <  2:101382/140‑1 (MQ=255)
                                          gAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTc                                                                                         <  1:912965/140‑1 (MQ=255)
                                           aGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCa                                                                                        <  1:1195584/140‑1 (MQ=255)
                                           aGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCa                                                                                        <  2:1469534/140‑1 (MQ=255)
                                             ttAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGt                                                                                      <  1:1078643/140‑1 (MQ=255)
                                             ttAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGt                                                                                      <  2:284251/140‑1 (MQ=255)
                                             ttAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGt                                                                                      <  1:1472918/140‑1 (MQ=255)
                                                accaccCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTa                                                                                   <  1:1003593/140‑1 (MQ=255)
                                                 ccaccCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTAt                                                                                  <  2:1468247/140‑1 (MQ=255)
                                                     ccTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAg                                                                              <  1:1073640/140‑1 (MQ=255)
                                                     ccTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAg                                                                              <  2:767073/140‑1 (MQ=255)
                                                     ccTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAg                                                                              <  1:344315/140‑1 (MQ=255)
                                                      cTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGc                                                                             >  1:1376638/1‑140 (MQ=255)
                                                      cTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGc                                                                             >  1:12045/1‑140 (MQ=255)
                                                      cTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGc                                                                             <  1:1061264/140‑1 (MQ=255)
                                                       tGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCa                                                                            <  1:338105/140‑1 (MQ=255)
                                                         ctctGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATg                                                                          <  1:1503237/140‑1 (MQ=255)
                                                              ccTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTa                                                                     >  1:59391/1‑140 (MQ=255)
                                                               cTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTAc                                                                    >  2:1097395/1‑140 (MQ=255)
                                                                   cAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGt                                                                <  1:18362/140‑1 (MQ=255)
                                                                       aaaGGTCACGAATTCTGGAAACGAGACTTTACGGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTgggg                                                            <  2:132805/140‑1 (MQ=255)
                                                                       aaaGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTTAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTgggg                                                            <  2:222741/140‑1 (MQ=255)
                                                                       aaaGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTgggg                                                            <  1:937591/140‑1 (MQ=255)
                                                                       aaaGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTgggg                                                            <  1:922317/140‑1 (MQ=255)
                                                                        aaGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGc                                                           <  2:442937/140‑1 (MQ=255)
                                                                        aaGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGc                                                           <  1:65317/140‑1 (MQ=255)
                                                                        aaGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGc                                                           >  1:1384972/1‑140 (MQ=255)
                                                                        aaGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGc                                                           >  2:707003/1‑140 (MQ=255)
                                                                           gTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCggg                                                        >  1:621005/1‑140 (MQ=255)
                                                                              aCGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTAt                                                     >  2:189372/1‑140 (MQ=255)
                                                                                gAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGa                                                   >  2:1225850/1‑140 (MQ=255)
                                                                                 aaTTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGaa                                                  <  2:1051203/140‑1 (MQ=255)
                                                                                    tCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCg                                               <  1:388833/140‑1 (MQ=255)
                                                                                    tCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCg                                               <  1:582594/140‑1 (MQ=255)
                                                                                             gagaCTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTg                                      <  1:1097395/140‑1 (MQ=255)
                                                                                                 cTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCaa                                  <  2:375277/140‑1 (MQ=255)
                                                                                                  tttACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCaaa                                 <  1:218545/140‑1 (MQ=255)
                                                                                                       aGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAAc                            <  1:891282/140‑1 (MQ=255)
                                                                                                       aGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAAc                            <  1:909670/140‑1 (MQ=255)
                                                                                                            gcagcGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGTTTCTGGCAAATCAACCAGaa                       >  2:1545991/1‑140 (MQ=255)
                                                                                                            gcagcGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGaa                       >  2:712513/1‑140 (MQ=255)
                                                                                                             cagcGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAAc                      >  2:1505841/1‑140 (MQ=255)
                                                                                                              agcGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACa                     >  1:843592/1‑140 (MQ=255)
                                                                                                                cGGGCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACata                   >  1:124415/1‑140 (MQ=255)
                                                                                                                   gCTATTTTCCTTTGTGCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATcg                >  1:773649/1‑140 (MQ=255)
                                                                                                                               tgtgCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATCGCCGCCATTCGcc    >  1:826392/1‑140 (MQ=255)
                                                                                                                               tgtgCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATCGCCGCCATTCGcc    >  2:888290/1‑140 (MQ=255)
                                                                                                                                 tgCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGATCTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGGGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATCGCCGCCATTCGCCca  <  1:211984/140‑1 (MQ=255)
                                                                                                                                 tgCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATCGCCGCCATTCGCCca  >  2:562688/1‑140 (MQ=255)
                                                                                                                                 tgCTTAAGGAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATCGCCGCCATTCGCCca  >  2:1284514/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                   
TGAAAGTGGCTGAATGGCTGGCAGAACATCCGCAAGTTGCGCGAGTTAACCACCCTGCTCTGCCTGGCAGTAAAGGTCACGAATTCTGGAAACGAGACTTTACAGGCAGCAGCGGGCTATTTTCCTTTGTGCTTAAGAAAAAACTCAATAATGAAGAGCTGGCGAACTATCTGGATAACTTCAGTTTATTCAGCATGGCCTACTCGTGGGGCGGGTATGAATCGTTGATCCTGGCAAATCAACCAGAACATATCGCCGCCATTCGCCCA  >  NC_000913/3153047‑3153315

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