Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,658,592 C→T 100% A175V (GCC→GTC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,5920CT100.0% 267.3 / ‑8.1 88A175V (GCC→GTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (30/55);  minor base A (0/2);  total (30/58)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAC  >  NC_000913/1658495‑1658688
                                                                                                 |                                                                                                
cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCt                                                                                                >  2:551253/1‑100 (MQ=255)
cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCt                                                                                                >  1:593681/1‑100 (MQ=255)
  gaagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcg                                                                                              <  2:289970/100‑1 (MQ=255)
    agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggc                                                                                            <  1:2497758/100‑1 (MQ=255)
      aaaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcgg                                                                                          >  1:1798818/1‑100 (MQ=255)
       aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGt                                                                                         <  1:2498869/100‑1 (MQ=255)
        acccGTCGAACAATATGGCAATGAGGCTGTATATACTCAGTACTCTTCGGGAATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcgg                                                                                          <  2:1216164/96‑1 (MQ=255)
            gTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAAc                                                                                    >  2:533271/1‑100 (MQ=255)
            gTCGAACAATATGGCAATGAGGCGGTAAATATTCAGGCCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAAc                                                                                    <  2:593674/100‑1 (MQ=255)
              cGAGCAGTATGGAGATGAGGGTGTAAATATACAGCACTGTTGGTGTATCTATTGCGGCAATATGACCCGCTGTTCGCCATCAGTCTCGGCGGTCaaa                                                                                     <  2:100159/97‑1 (MQ=255)
                aaCAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCt                                                                                >  2:542265/1‑100 (MQ=255)
                   aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgat                                                                             >  2:1452346/1‑100 (MQ=255)
                    atatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgatg                                                                            >  2:849261/1‑100 (MQ=255)
                    atatGGAGATGAGGCAGAACATATTCACAACTCTTCGGCGATCCTTTGCGGCTACAGGACCCGCTCTTCGCCGTCAGTCTCGGCGGTCAAACGCCtga                                                                              <  2:768588/98‑1 (MQ=255)
                      atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCATTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGaa                                                                          <  2:1416541/100‑1 (MQ=255)
                        ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACt                                                                        >  2:71906/1‑100 (MQ=255)
                        ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACt                                                                        >  1:1562208/1‑100 (MQ=255)
                         gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTg                                                                       >  1:583254/1‑100 (MQ=255)
                         gCAATGAGGCGGTATACATCCAATACTCTTCCGGGATCGTTGGCGGCAATAAGACCCGTTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgat                                                                             <  2:202187/94‑1 (MQ=255)
                         gCAATGAGGCGATATATGTTCATTACTCTTCCGATATCGTTGGCGGCAATAAGACCCGTTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgat                                                                             <  2:632835/94‑1 (MQ=255)
                          cAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGc                                                                      <  1:962440/100‑1 (MQ=255)
                          cAATGAGCATCTATATATTCAGCACTCTCCGGGGATCCTTGGCGGCAATATGACCCGCTCTTCGCCGTCAGTCTCGGCGGTCAAACGCCTGATGAACTGc                                                                      <  2:466074/100‑1 (MQ=255)
                           aaggAGGCTGAATATATTCATCCATCATTGTGGAACGGTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAAc                                                                         <  2:1771438/93‑1 (MQ=255)
                            aTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTa                                                                    <  1:1605747/100‑1 (MQ=255)
                                ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                                <  1:1336009/100‑1 (MQ=255)
                                    gTATATATTCAGTACTCTTCGGGGATCGTGGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcgg                                                              <  2:1075068/98‑1 (MQ=255)
                                     tatatcTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                           <  1:1014895/100‑1 (MQ=255)
                                      atatatCTAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                           <  2:1562197/99‑1 (MQ=255)
                                       tatatTCAGTACTCTTCGTGGATCGTTGGCGGCAATATGACCCGCTTTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGATACGACGGTTCAc                                                         >  2:958607/1‑100 (MQ=255)
                                        atatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACt                                                        <  1:1427418/100‑1 (MQ=255)
                                            tCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAc                                                    <  2:1448459/100‑1 (MQ=255)
                                            acagccctctTCTGAGCTCGGTAGCGGCAATAGGACCCGCTCTTCGCCATCAGTCTCGCCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                         <  2:709930/89‑1 (MQ=255)
                                             cAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAcc                                                   <  1:2075579/100‑1 (MQ=255)
                                             cAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAcc                                                   <  1:1561604/100‑1 (MQ=255)
                                               gTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAg                                                 >  2:625551/1‑100 (MQ=255)
                                                tACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGt                                                <  1:579318/100‑1 (MQ=255)
                                                  ctctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                               <  2:1510134/99‑1 (MQ=255)
                                                  ctctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTAt                                              >  1:633388/1‑100 (MQ=255)
                                                   tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                             <  2:1313371/100‑1 (MQ=255)
                                                   tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                             <  1:405231/100‑1 (MQ=255)
                                                   tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                             <  2:1313277/100‑1 (MQ=255)
                                                   tcgtCGGGGATCGTTAGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                             <  2:2058950/97‑1 (MQ=255)
                                                    ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                            <  1:1081662/100‑1 (MQ=255)
                                                    ctTCGGGGATCGTTAGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                            <  2:436011/100‑1 (MQ=255)
                                                       cGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTc                                         >  1:237692/1‑100 (MQ=255)
                                                        ggggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                        >  2:1337793/1‑100 (MQ=255)
                                                        ggggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                        <  2:1606853/100‑1 (MQ=255)
                                                        gcggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                        <  2:1552964/98‑1 (MQ=255)
                                                           gATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAc                                     <  1:1252119/100‑1 (MQ=255)
                                                             tCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                   <  1:399800/100‑1 (MQ=255)
                                                              cGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                  <  2:2427317/100‑1 (MQ=255)
                                                               gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                 >  2:170881/1‑100 (MQ=255)
                                                               gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                 <  1:542272/100‑1 (MQ=255)
                                                                ttGGCGGCAATATGACCCGCTCTTCGCCATCAGTCGCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCAc                                <  2:1890960/100‑1 (MQ=255)
                                                                  gggggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                    <  2:1097669/91‑1 (MQ=255)
                                                                  ggcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTg                              >  1:1100276/1‑100 (MQ=255)
                                                                   gcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                             >  2:2281697/1‑100 (MQ=255)
                                                                   gcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                             <  1:170883/100‑1 (MQ=255)
                                                                      gcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCa                          >  1:1205207/1‑100 (MQ=255)
                                                                      gcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCa                          <  1:1539320/100‑1 (MQ=255)
                                                                       cAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaa                         >  2:551694/1‑100 (MQ=255)
                                                                         atatGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAt                       >  1:1920209/1‑100 (MQ=255)
                                                                          tatGACGCGCCCTGAGCCAACCGACCCGGCGGTCAAACGCTTGATGCGCTGCTACGGCGGTTCACTCAACGAGTATGGCTGCTACAGCACTGCGCAAAtt                      <  2:176247/100‑1 (MQ=255)
                                                                          tatGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                      >  1:1066276/1‑100 (MQ=255)
                                                                            tGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTc                    <  1:1360695/100‑1 (MQ=255)
                                                                            tGACACGCTTTTCGCAATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCaaa                        <  2:841528/96‑1 (MQ=255)
                                                                             gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                   >  1:566282/1‑100 (MQ=255)
                                                                              aCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                  >  2:2213563/1‑100 (MQ=255)
                                                                               cccGCTCTTCGCCAACAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                  <  2:2071645/99‑1 (MQ=255)
                                                                                 cGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg               <  1:134146/100‑1 (MQ=255)
                                                                                 cGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg               <  1:1429132/100‑1 (MQ=255)
                                                                                 cGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgt                >  1:1954442/1‑99 (MQ=255)
                                                                                  gcgctTCGACATTAGACTCGGCGGTCAAACGCCTGAGGAACTGCTACGGCTGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc              <  2:2364437/97‑1 (MQ=255)
                                                                                  gCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCTTACAGCACTGCGCAAATTTCCTGTGc              <  1:1556270/100‑1 (MQ=255)
                                                                                   ctctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGcc             <  1:1321553/100‑1 (MQ=255)
                                                                                    tctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCa            <  1:20085/100‑1 (MQ=255)
                                                                                       tCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc         <  2:79292/100‑1 (MQ=255)
                                                                                       tCGCCATCAGTCTCGGCGGGCAAACGCCTGATGAACTGCTTCGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGCGCCATGc         <  2:129794/100‑1 (MQ=255)
                                                                                       tCGCCATCAGTCTAGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc         >  2:1588986/1‑100 (MQ=255)
                                                                                         gCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc       <  2:583247/100‑1 (MQ=255)
                                                                                         gCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAAAGCACTGCGCAAATTTCCTGTACCATGccc       >  2:536920/1‑100 (MQ=255)
                                                                                           cATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa     >  1:209324/1‑100 (MQ=255)
                                                                                           cATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa     >  1:1384767/1‑100 (MQ=255)
                                                                                             tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca   <  1:2087043/100‑1 (MQ=255)
                                                                                             tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca   <  1:2306631/100‑1 (MQ=255)
                                                                                             tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca   <  1:529956/100‑1 (MQ=255)
                                                                                             tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca   >  2:1820464/1‑100 (MQ=255)
                                                                                              cAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTGACTCAACCAGTGTGGCTCCTACAGCACTGCTCATATTTCCAGTGCCATGCTCTacac  <  2:95357/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAC  >  NC_000913/1658495‑1658688

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: