Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 179,442 A→T 100% H69L (CAC→CTC)  dgt → deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913179,4420AT100.0% 256.4 / ‑6.7 82H69L (CAC→CTC) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (42/39);  minor base G (1/0);  total (43/39)

GGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGAT  >  NC_000913/179346‑179536
                                                                                                |                                                                                              
tgCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCtctc                                                                                             >  2:374701/2‑100 (MQ=255)
ggCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCtctc                                                                                             <  1:400943/100‑1 (MQ=255)
ggCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCtctc                                                                                             <  2:1322349/100‑1 (MQ=255)
 gCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCg                                                                                            >  2:1627472/1‑100 (MQ=255)
   gTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAACACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGAt                                                                                          <  1:293990/100‑1 (MQ=255)
    tATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATg                                                                                         >  1:1005704/1‑100 (MQ=255)
      tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGa                                                                                       <  1:1716572/100‑1 (MQ=255)
      tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGa                                                                                       <  1:658251/100‑1 (MQ=255)
       catcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGACCGCAATGCCGCCTTGCGCACGCGTCTTACCCGCTCGCTGGa                                                                                       >  2:593007/1‑99 (MQ=255)
        atcaCCTCTCCGGCAATTCGTCGTCTTCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCCTGCGCACGCGTCTTACCCTCTCGATGGAAg                                                                                     <  2:632107/100‑1 (MQ=255)
        atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAg                                                                                     >  1:113872/1‑100 (MQ=255)
         tcaACTCCCCGGCAATTCGTCGTCCGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAg                                                                                     <  2:2461501/99‑1 (MQ=255)
          caACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTc                                                                                   <  1:2543081/100‑1 (MQ=255)
          caACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTc                                                                                   >  1:823256/1‑100 (MQ=255)
               ctcCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCcagca                                                                              <  1:335450/100‑1 (MQ=255)
               ctcCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCcagca                                                                              <  1:349179/100‑1 (MQ=255)
               ctcCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTCGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCcagc                                                                               <  2:592205/99‑1 (MQ=255)
                 ccGGCGAGTCGTCGTCTGCAACAAACGGCCCAGGTTTTTCCACTGAAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCcag                                                                                <  2:1925474/96‑1 (MQ=255)
                  cGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGt                                                                           <  1:2263157/100‑1 (MQ=255)
                   ggCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTg                                                                          <  1:231026/100‑1 (MQ=255)
                   ggCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACACTCTCGATGGAAGTCCAGCAGGTg                                                                          >  2:1821043/1‑100 (MQ=255)
                     cAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTggg                                                                        <  1:2107027/100‑1 (MQ=255)
                        ttCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGgcg                                                                     >  1:1834944/1‑100 (MQ=255)
                        ttCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGgcg                                                                     >  2:2361438/1‑100 (MQ=255)
                          cgtcgtCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCt                                                                   <  2:1991489/100‑1 (MQ=255)
                           gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCt                                                                   >  1:310909/1‑99 (MQ=255)
                           gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTa                                                                  >  1:1223265/1‑100 (MQ=255)
                           gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTa                                                                  >  2:1113614/1‑100 (MQ=255)
                           gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCCCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTa                                                                  <  2:1987014/100‑1 (MQ=255)
                             cgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACa                                                                >  1:1620/1‑100 (MQ=255)
                             cgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACa                                                                <  1:2224683/100‑1 (MQ=255)
                              gtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACAt                                                               >  1:2002823/1‑100 (MQ=255)
                              gtcTGCAACAAAAGACCCAGGTGTTTTTTTTTTTTTGCAACGCCGCCGTGCGCACGCGTCTTACCCTCTCGCTGGAAGTCCACCAGGGGGGCCGCTAc                                                                 >  2:1490293/1‑98 (MQ=255)
                                cTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCg                                                             >  2:608256/1‑100 (MQ=255)
                                 tGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGc                                                            >  1:2507495/1‑100 (MQ=255)
                                 tGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGc                                                            >  1:1750291/1‑100 (MQ=255)
                                 tGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGc                                                            >  1:401567/1‑100 (MQ=255)
                                  gCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGcc                                                           <  1:167/100‑1 (MQ=255)
                                  gCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTAAATCGcc                                                           >  1:1764097/1‑100 (MQ=255)
                                   caacaaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCa                                                          >  1:2202886/1‑100 (MQ=255)
                                     acaaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCaaa                                                        >  2:1281389/1‑100 (MQ=255)
                                      caaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAg                                                       >  2:1036889/1‑100 (MQ=255)
                                      caaAAGACCCAGGTTTTTACACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTCACCCTCTCGCTGGAAGTCCAGCAGGAGGGGAGCTACATCGCAAAAg                                                       >  2:879981/1‑100 (MQ=255)
                                        aaaGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGaa                                                     <  1:2035293/100‑1 (MQ=255)
                                         aaGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGaaa                                                    <  1:2252836/100‑1 (MQ=255)
                                             cccAGGTTTTTCCACTGTAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGCGCGCTACATCGCCAAAGAAAtt                                                  >  2:817721/1‑98 (MQ=255)
                                             cccAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAAtttt                                                <  2:779727/100‑1 (MQ=255)
                                             cccAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAAtttt                                                >  2:819642/1‑100 (MQ=255)
                                             cccAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAAtt                                                  >  2:674506/1‑98 (MQ=255)
                                              ccAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTa                                               <  2:1285627/100‑1 (MQ=255)
                                               cAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTaa                                              <  2:2507480/100‑1 (MQ=255)
                                                caggttttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAg                                             >  1:203865/5‑100 (MQ=255)
                                                  gTTTTTCCACTGGATCGCATTGCCGCCGTGCTCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAg                                             <  2:1524678/98‑1 (MQ=255)
                                                      ttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACGTCGCCAAAGAAATTTTAATCCGTCt                                       >  1:2300938/1‑100 (MQ=255)
                                                      ttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTc                                        <  2:2205206/99‑1 (MQ=255)
                                                      ttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCt                                       >  2:2483005/1‑100 (MQ=255)
                                                       tCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGACGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTg                                      >  2:2288332/1‑100 (MQ=255)
                                                         cACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGaa                                    >  2:1702060/1‑100 (MQ=255)
                                                          aCTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGaaa                                   >  1:2239559/1‑100 (MQ=255)
                                                           cTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAag                                  <  2:1629589/100‑1 (MQ=255)
                                                             ggAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAagag                                >  1:1569152/1‑100 (MQ=255)
                                                               aGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCt                              <  1:2413652/100‑1 (MQ=255)
                                                                 cgcAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTa                            <  2:1949095/100‑1 (MQ=255)
                                                                   cAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTaaa                          >  1:168925/1‑100 (MQ=255)
                                                                   cAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTaaa                          <  1:507544/100‑1 (MQ=255)
                                                                    aaTGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAAt                         <  1:1702072/100‑1 (MQ=255)
                                                                       gccgccGTGCGCGCGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTAc                      >  1:1682639/1‑100 (MQ=255)
                                                                       gccgccGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTAc                      <  1:203947/100‑1 (MQ=255)
                                                                        ccgccgTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACt                     >  2:575050/1‑100 (MQ=255)
                                                                         cgccgTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACt                     <  2:1630249/99‑1 (MQ=255)
                                                                          gccgTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGCGGCGCGCTACATCGCCAACGAAATTTTAAGCCGGCTGAAAGAGCTTAAATTACTgg                   >  2:1352323/1‑100 (MQ=255)
                                                                            cgTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGa                  <  2:637777/99‑1 (MQ=255)
                                                                            cgTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGATACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACt                     <  2:1598315/96‑1 (MQ=255)
                                                                             gTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGaa                 <  2:1620/99‑1 (MQ=255)
                                                                              tGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGc               <  1:2412971/100‑1 (MQ=255)
                                                                               gcgcgtgcgTCTGACCCTCACGATGGAAGTCCAGCATGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGc               <  2:204119/93‑1 (MQ=255)
                                                                                cgcACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCAt             >  1:348820/1‑100 (MQ=255)
                                                                                    cgcgTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACgg         >  1:1077934/1‑100 (MQ=255)
                                                                                    cgcgTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACgg         <  1:1584675/100‑1 (MQ=255)
                                                                                      cgTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGcc       <  1:819649/100‑1 (MQ=255)
                                                                                         cTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTgg    >  1:83192/1‑100 (MQ=255)
                                                                                         cTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTgg    <  1:360777/100‑1 (MQ=255)
                                                                                          ttACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGa   >  2:2162489/1‑100 (MQ=255)
                                                                                           tACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGGTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGATTAAATTACTGGAAGCATACGGCCTGGAt  >  2:1535920/1‑100 (MQ=255)
                                                                                           tACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGAt  >  1:601787/1‑100 (MQ=255)
                                                                                                |                                                                                              
GGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGAT  >  NC_000913/179346‑179536

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