Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,545,565 G→A *228* (TAG→TAA)  gntX → DNA catabolic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,545,5650GA100.0% 240.0 / NA 74*228* (TAG→TAA) gntXDNA catabolic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/39);  total (35/39)

CCCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGA  >  NC_000913/3545431‑3545682
                                                                                                                                      |                                                                                                                      
cccGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGcc                                                                                                                    >  2:188906/1‑139 (MQ=255)
       aaGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                             <  1:10818/139‑1 (MQ=255)
       aaGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                             <  1:188906/139‑1 (MQ=255)
       aaGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCgatga                                                                                                             <  2:872186/139‑1 (MQ=255)
           tCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                         <  1:1156404/139‑1 (MQ=255)
             accATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        <  2:118872/137‑1 (MQ=255)
             accATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        >  1:118872/2‑138 (MQ=255)
              ccATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        >  2:897695/1‑137 (MQ=255)
              ccATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        <  2:881488/137‑1 (MQ=255)
              ccATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        >  1:881488/1‑137 (MQ=255)
              ccATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        <  1:897695/137‑1 (MQ=255)
                     gTGATTGTGGATGATTTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACTCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATaa                                                                                               <  1:879455/139‑1 (MQ=255)
                            tgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                         >  2:394958/1‑122 (MQ=255)
                            tgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                         <  1:394958/122‑1 (MQ=255)
                            tgGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTa                                                                                       >  2:913666/1‑139 (MQ=255)
                                tgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        >  1:1236450/1‑119 (MQ=255)
                                tgatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGc                                                                                                        <  2:1236450/119‑1 (MQ=255)
                                 gatgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATa                                                                                  <  1:83733/139‑1 (MQ=255)
                                  atgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTNTCGAACNTTNTNAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAg                                                                                 >  1:1315728/1‑139 (MQ=255)
                                   tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTaaa                                                                                     <  2:748762/134‑1 (MQ=255)
                                   tgTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTaaa                                                                                     >  1:748762/1‑134 (MQ=255)
                                    gTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTc                                                                               <  1:913666/139‑1 (MQ=255)
                                    gTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTc                                                                               >  1:25533/1‑139 (MQ=255)
                                    gTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTc                                                                               >  2:898514/1‑139 (MQ=255)
                                      cGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCaa                                                                             <  1:1008178/139‑1 (MQ=255)
                                          accaccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                         <  1:10817/108‑1 (MQ=255)
                                          accaccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                         >  2:10817/1‑108 (MQ=255)
                                              ccGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAgg                                                                     >  2:1205250/1‑139 (MQ=255)
                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    <  2:25533/139‑1 (MQ=255)
                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    <  2:640820/139‑1 (MQ=255)
                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    <  1:1205250/139‑1 (MQ=255)
                                               cGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    <  1:1268685/139‑1 (MQ=255)
                                                 gAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGACGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCa                                                                  <  1:587269/139‑1 (MQ=255)
                                                         gTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATg                                                          >  1:318662/1‑139 (MQ=255)
                                                           cGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGAt                                                        <  1:898514/139‑1 (MQ=255)
                                                            gCAGAGATTTCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGa                                                         >  1:977451/1‑137 (MQ=255)
                                                            gCAGAGATTACGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGa                                                         <  2:977451/137‑1 (MQ=255)
                                                                      cgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACt                                                             >  1:136111/1‑123 (MQ=255)
                                                                      cgcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACt                                                             <  2:136111/123‑1 (MQ=255)
                                                                       gcAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGa                                            <  1:161417/139‑1 (MQ=255)
                                                                                     gCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAgcgc                              >  2:646271/1‑139 (MQ=255)
                                                                                        aTGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    <  2:966798/98‑1 (MQ=255)
                                                                                        aTGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    >  1:966798/1‑98 (MQ=255)
                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGcac                                    <  2:1191325/129‑1 (MQ=255)
                                                                                         tGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGcac                                    >  1:1191325/1‑129 (MQ=255)
                                                                                              cggcggCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCa                     <  2:1071106/139‑1 (MQ=255)
                                                                                              cggcggCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCa                     >  1:1071106/1‑139 (MQ=255)
                                                                                                 cggcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAAc                  >  1:14835/1‑139 (MQ=255)
                                                                                                   gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGAtgctg                                       >  2:933335/1‑117 (MQ=255)
                                                                                                   gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAACCCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGAtgctg                                       <  1:933335/117‑1 (MQ=255)
                                                                                                   gcTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAActg                >  1:373389/1‑139 (MQ=255)
                                                                                                      gTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAActgctg             >  2:1095883/1‑139 (MQ=255)
                                                                                                            gTCTGGTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAAt       <  2:318662/139‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGTTGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCaa         >  1:290847/1‑132 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGTTGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCaa         <  2:290847/132‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                         <  1:527511/37‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATggg                                                                                                         >  2:527511/1‑37 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                >  1:1062855/1‑46 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                <  2:368637/46‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                <  2:1062855/46‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATtata                                                                                                >  1:368637/1‑46 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTc                                                                               >  2:142903/1‑62 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTc                                                                               <  1:142903/62‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    <  1:933740/73‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    >  2:933740/1‑73 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATcc                                                      >  2:1091956/1‑87 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATcc                                                      <  1:1091956/87‑1 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga  >  2:353307/1‑139 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga  >  2:1208262/1‑139 (MQ=255)
                                                                                                                 gTGCCTTTGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGga  <  1:353307/139‑1 (MQ=255)
                                                                                                                 gTGCCTATGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGt                                                    >  2:1226755/1‑89 (MQ=255)
                                                                                                                 gTGCCTATGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGt                                                    <  1:1226755/89‑1 (MQ=255)
                                                                                                                       ttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    >  2:357253/1‑67 (MQ=255)
                                                                                                                       ttGTCGAACCTTGTAAAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGc                                                                    <  1:357253/67‑1 (MQ=255)
                                                                                                                                      |                                                                                                                      
CCCGTGCAAGGTCGCCATATGGTGATTGTGGATGATGTCGTTACCACCGGAAGTACCGTCGCAGAGATTGCGCAGTTGCTTTTACGCAATGGTGCGGCGGCTGTCCAGGTCTGGTGCCTTTGTCGAACCTTGTAGAGCCTCGATGATGGGCGTATTATAA‑CCAACTAAAATAGTCAACTATTAGGCCATTACTATGATCCGTATTTCCGATGCTGCACAAGCGCACTTTGCCAAACTGCTGGCAAATCAGGA  >  NC_000913/3545431‑3545682

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: ‑