Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 255,659 A→G C20R (TGT→CGT)  pepD ← aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913255,6590AG100.0% 268.7 / NA 73C20R (TGT→CGT) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/31);  total (42/31)

CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGA  >  NC_000913/255530‑255793
                                                                                                                                 |                                                                                                                                      
cGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                               >  1:823098/1‑139 (MQ=255)
cGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                               >  1:105534/1‑139 (MQ=255)
     gCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCaaaaa                                                                                                                          <  2:210547/139‑1 (MQ=255)
       aGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAata                                                                                                                        >  1:343452/1‑139 (MQ=255)
                   cAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGc                                                                                                            >  2:431588/1‑139 (MQ=255)
                      gATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                         <  1:431588/139‑1 (MQ=255)
                               tacctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAt                                                                                                <  1:775319/139‑1 (MQ=255)
                                acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                                                                                               >  2:283623/1‑139 (MQ=255)
                                acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                                                                                               >  1:549759/1‑139 (MQ=255)
                                  ctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                     <  2:356479/131‑1 (MQ=255)
                                  ctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                     >  1:356479/1‑131 (MQ=255)
                                  ctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTga                                                                                             >  2:907433/1‑139 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                          >  2:648386/1‑139 (MQ=255)
                                     aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACa                                                                                         >  2:507301/1‑139 (MQ=255)
                                       cgcgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAg                                                                                        >  1:553496/1‑139 (MQ=255)
                                       cgcgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAg                                                                                        >  1:219975/1‑139 (MQ=255)
                                         cgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGtt                                                                                      <  2:92366/139‑1 (MQ=255)
                                           ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                           <  2:834644/132‑1 (MQ=255)
                                           ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                           >  1:834644/1‑132 (MQ=255)
                                           ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                                    >  2:297346/1‑139 (MQ=255)
                                               aCATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGaca                                                                                <  2:553496/139‑1 (MQ=255)
                                                         cTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGtct                                                                      <  1:3590/139‑1 (MQ=255)
                                                             ctctTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCtt                                                                  >  2:864426/1‑139 (MQ=255)
                                                                   tGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAatat                                                                                                                       <  2:1005096/80‑1 (MQ=255)
                                                                   tGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAatat                                                                                                                       >  1:1005096/1‑80 (MQ=255)
                                                                    gCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATc                                                           <  1:283623/139‑1 (MQ=255)
                                                                    gCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATc                                                           <  2:823098/139‑1 (MQ=255)
                                                                     cccAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCa                                                          >  2:874818/1‑139 (MQ=255)
                                                                         accaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                      <  2:549759/139‑1 (MQ=255)
                                                                                   tATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                     >  1:418886/1‑82 (MQ=255)
                                                                                   tATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                     <  2:418886/82‑1 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                         <  2:274575/76‑1 (MQ=255)
                                                                                     ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                         >  1:274575/1‑76 (MQ=255)
                                                                                      tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                      >  1:730275/1‑126 (MQ=255)
                                                                                      tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                      <  2:730275/126‑1 (MQ=255)
                                                                                       cAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                        <  1:297346/139‑1 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       <  1:874818/139‑1 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       <  1:106071/139‑1 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       >  2:657898/1‑139 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       >  2:93064/1‑139 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       >  1:399494/1‑139 (MQ=255)
                                                                                        aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       >  2:1040173/1‑139 (MQ=255)
                                                                                         gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAg                                              <  1:586159/131‑1 (MQ=255)
                                                                                         gCGAGTTGCTCTTCATGAAAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAg                                              >  2:586159/1‑131 (MQ=255)
                                                                                                   ctTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGATATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGt                            <  1:648386/139‑1 (MQ=255)
                                                                                                         gATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAg                      <  1:16834/139‑1 (MQ=255)
                                                                                                           tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg                                                                         <  1:190021/86‑1 (MQ=255)
                                                                                                           tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAg                                                                         >  2:190021/1‑86 (MQ=255)
                                                                                                           tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTa                                            >  1:487992/1‑115 (MQ=255)
                                                                                                           tAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTa                                            <  2:487992/115‑1 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                 >  1:556618/1‑139 (MQ=255)
                                                                                                              gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                 >  2:256723/1‑139 (MQ=255)
                                                                                                              gACGGGAGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGGGGAGATAATTGAGCCAGTTCAGACACGATAAGTCCCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                 >  2:945198/1‑139 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt               >  1:36533/1‑139 (MQ=255)
                                                                                                                cGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgt               >  1:775051/1‑139 (MQ=255)
                                                                                                                 gggTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgtt              <  2:105685/139‑1 (MQ=255)
                                                                                                                 gggTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgtt              >  1:105685/1‑139 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                                    <  2:814131/67‑1 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                                    >  1:814131/1‑67 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt            >  1:1050765/1‑139 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt            >  2:45227/1‑139 (MQ=255)
                                                                                                                   gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGAAGGTCACAGGGTTAGCAGAAAAtgttgt            >  2:460060/1‑139 (MQ=255)
                                                                                                                        gAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca       <  1:657898/139‑1 (MQ=255)
                                                                                                                        gAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca       <  1:495490/139‑1 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                           >  1:948542/1‑50 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                           <  2:948542/50‑1 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGATAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca       <  1:562468/134‑1 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       <  2:552989/102‑1 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                       >  1:552989/1‑102 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca       >  2:562468/1‑134 (MQ=255)
                                                                                                                             gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGa  >  1:189691/1‑139 (MQ=255)
                                                                                                                             gAACGGATTTTGACAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca       >  2:157537/1‑134 (MQ=255)
                                                                                                                             gAACGGATTTTGACAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca       <  1:157537/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                      
CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGA  >  NC_000913/255530‑255793

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