Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,926,782 G→A *221* (TAG→TAA
D686D (GAC→GAT
exoX →
ptrB ←
exodeoxyribonuclease 10, DNA exonuclease X
protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,926,7820GA100.0% 241.2 / NA 74*221* (TAG→TAA
D686D (GAC→GAT
exoX
ptrB
exodeoxyribonuclease 10, DNA exonuclease X
protease II
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/41);  total (33/41)

TGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAA  >  NC_000913/1926647‑1926917
                                                                                                                                       |                                                                                                                                       
tGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAAAGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATcc                                                                                                                                      <  2:1045132/139‑1 (MQ=255)
 ggCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCg                                                                                                                                     <  1:1519531/139‑1 (MQ=255)
     aaTACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAgg                                                                                                                                 >  2:972049/1‑139 (MQ=255)
      aTACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGc                                                                                                                                <  1:1567599/139‑1 (MQ=255)
           gTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                                           <  1:1207506/139‑1 (MQ=255)
             ggCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCgg                                                                                                                         <  1:612725/139‑1 (MQ=255)
                   gCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATg                                                                                                                   <  1:972049/139‑1 (MQ=255)
                      gTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTc                                                                                                                >  1:147478/1‑139 (MQ=255)
                       tttCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTcc                                                                                                               <  1:1241088/139‑1 (MQ=255)
                         tCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCt                                                                                                             >  2:1618022/1‑139 (MQ=255)
                            gACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTgcg                                                                                                          <  2:194302/139‑1 (MQ=255)
                                   ccGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgcgc                                                                                                   >  2:1437726/1‑139 (MQ=255)
                                      aaCGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                                <  1:599711/139‑1 (MQ=255)
                                      aaCGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                                <  1:482057/139‑1 (MQ=255)
                                           gATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaa                                                                                           <  2:461767/139‑1 (MQ=255)
                                           gATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaa                                                                                           <  2:1634614/139‑1 (MQ=255)
                                                 ggCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                                <  1:798223/128‑1 (MQ=255)
                                                 ggCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                                >  2:798223/1‑128 (MQ=255)
                                                      tCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGt                                                                                                                              <  2:280633/93‑1 (MQ=255)
                                                      tCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGt                                                                                                                              >  1:280633/1‑93 (MQ=255)
                                                       cTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCa                                                                               >  2:379934/1‑139 (MQ=255)
                                                          cgcTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCg                                                                            <  2:147478/139‑1 (MQ=255)
                                                            cTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCt                                                                          >  2:800691/1‑139 (MQ=255)
                                                                  aTTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                    <  1:379934/139‑1 (MQ=255)
                                                                       ataaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTa                                                               <  2:563781/139‑1 (MQ=255)
                                                                         aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                    >  1:1453060/1‑132 (MQ=255)
                                                                         aaCCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                    <  2:1453060/132‑1 (MQ=255)
                                                                           ccTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGAt                                                           >  1:99653/1‑139 (MQ=255)
                                                                             tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCg                                                                                                                                     <  2:994636/63‑1 (MQ=255)
                                                                             tGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCg                                                                                                                                     >  1:994636/1‑63 (MQ=255)
                                                                              ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCata                                                                                    >  1:1395217/1‑111 (MQ=255)
                                                                              ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCata                                                                                    <  2:1395217/111‑1 (MQ=255)
                                                                              ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                    <  2:728796/127‑1 (MQ=255)
                                                                              ggACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCt                                                                    >  1:728796/1‑127 (MQ=255)
                                                                               gACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgcgc                                                                                                   <  1:1004176/95‑1 (MQ=255)
                                                                               gACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAgcgc                                                                                                   >  2:1004176/1‑95 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATg                                                                                                                   >  1:574082/1‑74 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATg                                                                                                                   <  2:574082/74‑1 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                                >  2:289653/1‑93 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGAc                                                                                                <  1:289653/93‑1 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTc                                                                                 >  1:984096/1‑108 (MQ=255)
                                                                                    cATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTc                                                                                 <  2:984096/108‑1 (MQ=255)
                                                                                     aTGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTaa                                                       <  2:66763/133‑1 (MQ=255)
                                                                                     aTGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTaa                                                       >  1:66763/1‑133 (MQ=255)
                                                                                      tGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGAcc                                                >  2:1035632/1‑139 (MQ=255)
                                                                                         gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAg                                                                                        <  1:306630/96‑1 (MQ=255)
                                                                                         gCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAg                                                                                        >  2:306630/1‑96 (MQ=255)
                                                                                             ggAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTa                                         <  1:1618022/139‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTa                                                                                                                      >  1:815113/1‑61 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTa                                                                                                                      <  2:815113/61‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaa                                                                                           <  1:865973/88‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaa                                                                                           >  2:865973/1‑88 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGt                                                                >  2:368590/1‑115 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGt                                                                <  1:368590/115‑1 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                                        >  2:945243/1‑139 (MQ=255)
                                                                                              gAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTAc                                        >  1:918737/1‑139 (MQ=255)
                                                                                                         aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                             <  2:216410/139‑1 (MQ=255)
                                                                                                         aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                             >  1:216410/1‑139 (MQ=255)
                                                                                                         aaCACTGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCt                             <  2:875575/139‑1 (MQ=255)
                                                                                                              tGAAACATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTc                        <  2:628371/139‑1 (MQ=255)
                                                                                                                aaaCATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                                           >  2:1619906/1‑38 (MQ=255)
                                                                                                                aaaCATTATCTGGAAAATACTTAATCCGCAGGCGTAGc                                                                                                                           <  1:1619906/38‑1 (MQ=255)
                                                                                                                aaaCATTATCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCa                      >  2:1198619/1‑139 (MQ=255)
                                                                                                                        tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaaa                                                                                          <  1:1116283/63‑1 (MQ=255)
                                                                                                                        tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGaaa                                                                                          >  2:1116283/1‑63 (MQ=255)
                                                                                                                        tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAg                                                     <  2:1432341/100‑1 (MQ=255)
                                                                                                                        tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAg                                                     >  1:1432341/1‑100 (MQ=255)
                                                                                                                        tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa              >  2:1627813/1‑139 (MQ=255)
                                                                                                                        tCTGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTa              <  1:800691/139‑1 (MQ=255)
                                                                                                                          tGGAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACa            >  2:998036/1‑139 (MQ=255)
                                                                                                                            gAAAATACTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACaga          <  1:1198619/139‑1 (MQ=255)
                                                                                                                                  aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                   <  2:812166/124‑1 (MQ=255)
                                                                                                                                  aCTTAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGt                   >  1:812166/1‑124 (MQ=255)
                                                                                                                                    ttAATCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATaaa  <  2:687250/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
TGGCAAATACCGTGGCAAAGCGGTTTCCGACGTTGCCGAACGCGATCCGGGCTATCTGCGCTGGTTATTTAATAACCTGGACAGCATGAGCCCGGAGCTGCGTTTAACACTGAAACATTATCTGGAAAATACTTAGTCCGCAGGCGTAGCGGGTAATGTTCCCTGCGCCAGCGCGACCAGAAAAGCATATTCCATCGCTACGCCTTCGTACGATTTAAAGCGACCAGATTTACCGCCATGGCCTGAGTCCATGTCGGTACAGAGCAATAAA  >  NC_000913/1926647‑1926917

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