Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,812,482 T→C F198L (TTC→CTC)  katE → catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,812,4820TC100.0% 514.8 / NA 156F198L (TTC→CTC) katEcatalase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (66/90);  total (66/90)

Only 100 of 159 total aligned reads displayed.

ACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACTTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGG  >  NC_000913/1812392‑1812571
                                                                                          |                                                                                         
  cGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGt                                                                                       <  1:2419730/93‑1 (MQ=255)
   ggtCAGGGTCGTGGTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATGCCGAAGAGGGTATTTTTGACCGCGtt                                                                                      <  2:882421/91‑1 (MQ=255)
   gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                                      <  2:3563710/93‑1 (MQ=255)
   gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                                      <  1:2752476/93‑1 (MQ=255)
   gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                                      <  2:3513360/93‑1 (MQ=255)
   gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                                      <  2:2371179/93‑1 (MQ=255)
   gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                                      <  1:4676641/93‑1 (MQ=255)
   gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                                      <  1:4402495/93‑1 (MQ=255)
        gggTGGTGCTGGCTCTGCTGATACCGTGCGTGATATGCGTGGCTGTGTCACCAAGTTCTATACCGAAGAGTGTATTTTTGACCTCGTTGGcga                                                                                 >  2:633599/1‑91 (MQ=255)
        gggTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGGATTTTTGACCTCGTTGGCaa                                                                                 >  1:4268375/1‑93 (MQ=255)
          gtggtgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATa                                                                               >  2:2280901/1‑93 (MQ=255)
           tggtgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATaa                                                                              <  1:2756571/93‑1 (MQ=255)
            ggtgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAac                                                                             <  1:1148911/93‑1 (MQ=255)
            ggtgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAac                                                                             >  1:1571124/1‑93 (MQ=255)
                cTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGc                                                                         >  1:3550379/1‑93 (MQ=255)
                  ggCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCa                                                                       >  2:927833/1‑93 (MQ=255)
                  ggCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCa                                                                       >  2:355825/1‑93 (MQ=255)
                   gCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCaa                                                                      <  2:3058129/93‑1 (MQ=255)
                   gCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCaa                                                                      <  2:10239/93‑1 (MQ=255)
                   gCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCaa                                                                      >  1:3451107/1‑93 (MQ=255)
                     tctGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtc                                                                    <  1:2314443/93‑1 (MQ=255)
                      ctgctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtct                                                                   <  1:1700918/93‑1 (MQ=255)
                       tgctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtctt                                                                  >  1:894174/1‑93 (MQ=255)
                        tctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtcttc                                                                 <  1:2935301/92‑1 (MQ=255)
                        gctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtcttc                                                                 <  2:1182475/93‑1 (MQ=255)
                        gctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtcttc                                                                 >  1:2029299/1‑93 (MQ=255)
                            atgccGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTa                                                             <  1:1759567/90‑1 (MQ=255)
                            aTACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTa                                                             >  2:761649/1‑93 (MQ=255)
                            aTACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTa                                                             <  1:3719813/93‑1 (MQ=255)
                            aTACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTa                                                             <  1:3719815/93‑1 (MQ=255)
                             tACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTAt                                                            >  1:1838975/1‑93 (MQ=255)
                               ccGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATcc                                                          <  2:1516883/93‑1 (MQ=255)
                               ccGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATcc                                                          <  1:2504502/93‑1 (MQ=255)
                                cgtgcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCttt                                                              <  2:4136183/88‑1 (MQ=255)
                                cgtgcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCttt                                                              >  1:4136183/1‑88 (MQ=255)
                                cgtgcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCa                                                         <  1:1271364/93‑1 (MQ=255)
                                cgtgcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCa                                                         <  1:328660/93‑1 (MQ=255)
                                 gtgcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAg                                                        <  1:3937343/93‑1 (MQ=255)
                                   gcgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGa                                                      <  2:1858351/93‑1 (MQ=255)
                                     gtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATg                                                    >  2:1355692/1‑93 (MQ=255)
                                     gtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATg                                                    <  2:648094/93‑1 (MQ=255)
                                      tgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgc                                                   <  2:2455940/93‑1 (MQ=255)
                                       gATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcg                                                  <  1:2037365/93‑1 (MQ=255)
                                       gATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcg                                                  >  2:3371664/1‑93 (MQ=255)
                                         tatCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCa                                                >  2:3229122/1‑93 (MQ=255)
                                         tatCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCa                                                >  1:3427845/1‑93 (MQ=255)
                                         tatCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCa                                                >  1:4232567/1‑93 (MQ=255)
                                           tCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATa                                              >  2:3152569/1‑93 (MQ=255)
                                           tCCGGGGCTTTGCCACCAAGTGCGATACCGAAGAGGGGAGGCGTGACCTCGCCGGCAACAACACGCCAAGCTTCGGTAGCCAGGATGCGCaca                                              >  1:3607009/1‑91 (MQ=255)
                                            ccGTTGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATaa                                             <  1:4043438/93‑1 (MQ=255)
                                            ccGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATaa                                             <  1:4333309/93‑1 (MQ=255)
                                             cGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATaaa                                            >  2:2873365/1‑93 (MQ=255)
                                             cGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATaaa                                            >  1:4385492/1‑93 (MQ=255)
                                              gTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGGAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAAt                                           <  2:48343/93‑1 (MQ=255)
                                              gTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAAt                                           >  2:3523299/1‑93 (MQ=255)
                                               tGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAAtt                                          >  2:3023870/1‑93 (MQ=255)
                                                ggCTTTGCCTCCAAGTTCTATACCGAAGAGGGTATTTTTGCCCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTc                                         <  2:4381640/93‑1 (MQ=255)
                                                 gCTTTGAAACCAAGTTCTATACAGCAGAGGGCATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcc                                        >  2:4233523/1‑93 (MQ=255)
                                                  ctctGCCACCAAGTTCTATCCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAGCTTCTTTATCCAGGATGCGCATAAATTccc                                       <  2:3077859/90‑1 (MQ=255)
                                                  cTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTccc                                       <  1:310776/93‑1 (MQ=255)
                                                  cTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTccc                                       <  1:3666060/93‑1 (MQ=255)
                                                   tttGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATATATTcccc                                      >  1:2128507/1‑93 (MQ=255)
                                                   tttGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcccc                                      <  2:1206184/93‑1 (MQ=255)
                                                     tGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGa                                    <  2:2266723/93‑1 (MQ=255)
                                                      gCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGAt                                   <  2:1304584/93‑1 (MQ=255)
                                                       ccaccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGAtt                                  <  1:337924/93‑1 (MQ=255)
                                                         accaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGAtttt                                >  2:4552399/1‑93 (MQ=255)
                                                         accaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGAtttt                                >  2:4552400/1‑93 (MQ=255)
                                                          ccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTg                               >  2:1860033/1‑93 (MQ=255)
                                                            aaGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGtt                             <  1:2951442/93‑1 (MQ=255)
                                                            aaGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGtt                             >  2:3581590/1‑93 (MQ=255)
                                                            aaGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGtt                             >  2:2979298/1‑93 (MQ=255)
                                                             aGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCACGATTTTGTTc                            >  2:42197/1‑93 (MQ=255)
                                                              gTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCa                           >  1:4798/1‑93 (MQ=255)
                                                              gTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCa                           <  1:4559622/93‑1 (MQ=255)
                                                              gTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCa                           <  2:1155808/93‑1 (MQ=255)
                                                              gTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCa                           <  2:4495147/93‑1 (MQ=255)
                                                               ttCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCAt                          <  1:1576208/93‑1 (MQ=255)
                                                                tCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATg                         <  1:4491199/93‑1 (MQ=255)
                                                                  tataCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCg                       >  2:3639605/1‑93 (MQ=255)
                                                                    taCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGt                     >  2:1156555/1‑93 (MQ=255)
                                                                     aCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTa                    <  1:970698/93‑1 (MQ=255)
                                                                     aCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTa                    <  2:4665982/93‑1 (MQ=255)
                                                                      ccGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTaa                   >  1:4569939/1‑93 (MQ=255)
                                                                        gAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTaaaa                 <  1:412458/93‑1 (MQ=255)
                                                                           gagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCa              >  1:2355181/1‑93 (MQ=255)
                                                                           gagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCa              >  2:105171/1‑93 (MQ=255)
                                                                            agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAg             <  1:2419796/93‑1 (MQ=255)
                                                                            agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAg             <  2:445188/93‑1 (MQ=255)
                                                                            agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAg             <  2:597470/93‑1 (MQ=255)
                                                                            agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAg             <  1:1029852/93‑1 (MQ=255)
                                                                              ggTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGaa           <  2:1073883/93‑1 (MQ=255)
                                                                               gTATTTTTGACCTCGTTGGCAATAACCCGCCAATCTGCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAAc          <  1:4289168/93‑1 (MQ=255)
                                                                               gTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAAc          <  2:3292123/93‑1 (MQ=255)
                                                                               gTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAAc          <  2:2135804/93‑1 (MQ=255)
                                                                                 aTTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCg        >  1:1016398/1‑93 (MQ=255)
                                                                                  tttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGc       <  2:4122179/93‑1 (MQ=255)
                                                                                  tttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGc       <  2:2684670/93‑1 (MQ=255)
                                                                                   ttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCa      <  1:3996776/93‑1 (MQ=255)
                                                                                     ttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACt    >  1:541117/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACTTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGG  >  NC_000913/1812392‑1812571

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