Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 383,259 C→T R226Q (CGG→CAG)  yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,2590CT100.0% 275.1 / NA 83R226Q (CGG→CAG) yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/44);  total (39/44)

ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAA  >  NC_000913/383123‑383382
                                                                                                                                        |                                                                                                                           
aCCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCGCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGa                                                                                                                           >  2:41589/1‑139 (MQ=255)
       cGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGt                                                                                                                    >  1:1370281/1‑139 (MQ=255)
         gCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAg                                                                                                                  >  2:797398/1‑139 (MQ=255)
                 aaaGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCGCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                          <  1:41589/139‑1 (MQ=255)
                      cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     <  2:811795/139‑1 (MQ=255)
                      cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     <  2:1324014/139‑1 (MQ=255)
                      cTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     <  2:1370281/139‑1 (MQ=255)
                                ggCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                           <  1:797398/139‑1 (MQ=255)
                                ggCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                           >  1:626101/1‑139 (MQ=255)
                                              gcgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                             <  1:608682/139‑1 (MQ=255)
                                               cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCg                                                                                                               <  1:1005933/104‑1 (MQ=255)
                                               cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCg                                                                                                               >  2:1005933/1‑104 (MQ=255)
                                               cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     >  2:130174/1‑114 (MQ=255)
                                               cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     <  1:130174/114‑1 (MQ=255)
                                               cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                            >  2:1192787/1‑139 (MQ=255)
                                               cgTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTc                                                                            >  2:737974/1‑139 (MQ=255)
                                                 tATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGa                                                                          <  2:149795/139‑1 (MQ=255)
                                                   tCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                        <  2:674615/139‑1 (MQ=255)
                                                        caATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGt                                                                   <  2:1042471/139‑1 (MQ=255)
                                                                     ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     >  2:617452/1‑92 (MQ=255)
                                                                     ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     <  1:617452/92‑1 (MQ=255)
                                                                     ccAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGaaa                                                      >  2:1261818/1‑139 (MQ=255)
                                                                              ccACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                             <  1:1422541/139‑1 (MQ=255)
                                                                                aCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                           <  2:1050105/139‑1 (MQ=255)
                                                                                aCTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACg                                           <  1:1261818/139‑1 (MQ=255)
                                                                                      cACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                           <  2:411399/85‑1 (MQ=255)
                                                                                      cACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                           >  1:411399/1‑85 (MQ=255)
                                                                                            ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                             >  1:310963/1‑93 (MQ=255)
                                                                                            ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                             <  2:310963/93‑1 (MQ=255)
                                                                                            ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                   >  1:728696/1‑135 (MQ=255)
                                                                                            ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                   <  2:728696/135‑1 (MQ=255)
                                                                                            ggCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGt                               >  2:412008/1‑139 (MQ=255)
                                                                                                     aTTCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACAt                      <  2:226808/139‑1 (MQ=255)
                                                                                                      ttCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATc                     >  2:1213153/1‑139 (MQ=255)
                                                                                                       tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                           >  1:754640/1‑68 (MQ=255)
                                                                                                       tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGt                                                                                           <  2:754640/68‑1 (MQ=255)
                                                                                                       tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                        <  2:592974/119‑1 (MQ=255)
                                                                                                       tCCTGTGGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                        >  1:592974/1‑119 (MQ=255)
                                                                                                           gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                        <  2:219354/83‑1 (MQ=255)
                                                                                                           gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACt                                                                        >  1:219354/1‑83 (MQ=255)
                                                                                                           gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                        >  2:978318/1‑115 (MQ=255)
                                                                                                           gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCaa                                        <  1:978318/115‑1 (MQ=255)
                                                                                                           gtgGGTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCa                >  1:9408/1‑139 (MQ=255)
                                                                                                             gggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAgg              <  1:1409644/139‑1 (MQ=255)
                                                                                                              ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                <  1:1206588/56‑1 (MQ=255)
                                                                                                              ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                >  2:1206588/1‑56 (MQ=255)
                                                                                                              ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                >  2:654578/1‑56 (MQ=255)
                                                                                                              ggTATACGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                <  1:654578/56‑1 (MQ=255)
                                                                                                                 ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     <  1:418897/48‑1 (MQ=255)
                                                                                                                 ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     >  2:418897/1‑48 (MQ=255)
                                                                                                                 ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGc                                          >  1:331995/1‑107 (MQ=255)
                                                                                                                 ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGc                                          <  2:331995/107‑1 (MQ=255)
                                                                                                                 ataCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCAcc          >  1:734295/1‑139 (MQ=255)
                                                                                                                  taCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                             <  2:1340627/71‑1 (MQ=255)
                                                                                                                  taCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCtt                                                                             >  1:1340627/1‑71 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAAccttgat                                               <  1:1117187/100‑6 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAAccttgat                                               >  2:1117187/1‑95 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                   <  1:256638/112‑1 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACAt                                   >  2:256638/1‑112 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCa            >  1:746641/1‑135 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCa            <  2:746641/135‑1 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc        <  1:1213153/139‑1 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc        <  2:1473170/139‑1 (MQ=255)
                                                                                                                   aCGATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGc        <  2:425345/139‑1 (MQ=255)
                                                                                                                     gATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTg      <  1:938541/139‑1 (MQ=255)
                                                                                                                     gATGCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTg      <  2:626101/139‑1 (MQ=255)
                                                                                                                        gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                             <  2:241172/97‑1 (MQ=255)
                                                                                                                        gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTa                                             >  1:241172/1‑97 (MQ=255)
                                                                                                                        gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGata   >  2:1001855/1‑139 (MQ=255)
                                                                                                                        gCAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGata   >  1:458180/1‑139 (MQ=255)
                                                                                                                         ccgCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                          >  1:882068/3‑35 (MQ=255)
                                                                                                                         cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGc                                                                                                          <  2:882068/35‑1 (MQ=255)
                                                                                                                         cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGt                                                                   <  1:32980/74‑1 (MQ=255)
                                                                                                                         cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGt                                                                   >  2:32980/1‑74 (MQ=255)
                                                                                                                         cAGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGataa  >  2:253432/1‑139 (MQ=255)
                                                                                                                          aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     <  1:329149/39‑1 (MQ=255)
                                                                                                                          aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGccacca                                                                                                     >  2:329149/1‑39 (MQ=255)
                                                                                                                          aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                >  1:14265/1‑44 (MQ=255)
                                                                                                                          aGCGATTCGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTc                                                                                                <  2:14265/44‑1 (MQ=255)
                                                                                                                                 cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATa                                      >  1:1050137/1‑95 (MQ=255)
                                                                                                                                 cGCCTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATa                                      <  2:1050137/95‑1 (MQ=255)
                                                                                                                                   ccTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCga                                                          <  2:788365/73‑1 (MQ=255)
                                                                                                                                   ccTGCTGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCga                                                          >  1:788365/1‑73 (MQ=255)
                                                                                                                                        |                                                                                                                           
ACCGAAGCGGCTAAATAAAAGTCTTTTATGCAGGCGCATACAGACCGCGTATCCCACAATCCAACGCCGCCAGCGACGCCACTCCTCACGCGGGCTGTTGCATTCCTGTGGGTATACGATGCAGCGATTCGCCTGCCGAATACGGTAGCCGTTTGCCACCAATGTCCAGGTGAGATCAAGGTCTTCGACTTTAGTACGATCCGAGAAACCGAACTTACGCAATACATCAGTACGGAACATCCCGCAGGCACCGCTGATAA  >  NC_000913/383123‑383382

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