Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,248,939 T→C intergenic (+117/‑30) mzrA → / → yqjC modulator of EnvZ/OmpR regulon/DUF1090 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,248,9390TC100.0% 284.4 / NA 78intergenic (+117/‑30)mzrA/yqjCmodulator of EnvZ/OmpR regulon/DUF1090 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/42);  total (36/42)

ATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGG  >  NC_000913/3248809‑3249068
                                                                                                                                  |                                                                                                                                 
aTCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGa                                                                                                                           >  2:360980/1‑139 (MQ=255)
aTCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGa                                                                                                                           >  2:444740/1‑139 (MQ=255)
 tCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGaa                                                                                                                          <  2:840328/139‑1 (MQ=255)
           tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                <  2:737203/139‑1 (MQ=255)
           tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                <  2:711417/139‑1 (MQ=255)
           tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                <  1:444740/139‑1 (MQ=255)
              cTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAgg                                                                                                             >  1:267162/1‑139 (MQ=255)
                   tGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                         >  1:675905/1‑122 (MQ=255)
                   tGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                         <  2:675905/122‑1 (MQ=255)
                   tGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                        <  2:25310/139‑1 (MQ=255)
                   tGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                        >  2:496850/1‑139 (MQ=255)
                        tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATg                                                                                                   >  2:978375/1‑139 (MQ=255)
                             aaCATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                <  1:925455/121‑1 (MQ=255)
                             aaCATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                >  2:925455/1‑121 (MQ=255)
                                  atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCa                                                                                         >  1:176809/1‑139 (MQ=255)
                                      ctcactTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCAt                                                                                        <  2:825561/136‑1 (MQ=255)
                                      ctcactTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCAt                                                                                        >  1:825561/1‑136 (MQ=255)
                                         actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                  >  2:1102341/1‑139 (MQ=255)
                                          ctTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATc                                                                                                                        <  2:1023850/100‑1 (MQ=255)
                                          ctTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATc                                                                                                                        >  1:1023850/1‑100 (MQ=255)
                                           tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                >  1:430934/1‑139 (MQ=255)
                                                tGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGtt                                                                           <  2:343903/139‑1 (MQ=255)
                                                     tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGa                                                                                                                           <  1:707629/86‑1 (MQ=255)
                                                     tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGa                                                                                                                           >  2:707629/1‑86 (MQ=255)
                                                     tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtctct                                                                      >  2:184462/1‑139 (MQ=255)
                                                       aCCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTctctct                                                                    <  1:978375/139‑1 (MQ=255)
                                                                 tCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAg                                                          <  2:664972/139‑1 (MQ=255)
                                                                  cTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGt                                                         <  2:367197/139‑1 (MQ=255)
                                                                     tGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCttt                                                                  <  2:204967/127‑1 (MQ=255)
                                                                     tGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCttt                                                                  >  1:204967/1‑127 (MQ=255)
                                                                     tGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGcc                                                      <  2:176809/139‑1 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATg                                                                                                                  >  2:850672/1‑73 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATg                                                                                                                  <  1:850672/73‑1 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                          <  2:1132811/81‑1 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                          >  1:1132811/1‑81 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGCAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGcc                                                      <  2:800257/133‑1 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGcc                                                      >  1:800257/1‑133 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                <  2:1116189/139‑1 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                >  2:907208/1‑139 (MQ=255)
                                                                                gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGtt                                                                                                           >  2:823489/1‑75 (MQ=255)
                                                                                gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGtt                                                                                                           <  1:823489/75‑1 (MQ=255)
                                                                                 cGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCtt                                                            >  2:444162/1‑121 (MQ=255)
                                                                                 cGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCtt                                                            <  1:444162/121‑1 (MQ=255)
                                                                                   gtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAg                                                          <  2:111207/121‑1 (MQ=255)
                                                                                   gtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAg                                                          >  1:111207/1‑121 (MQ=255)
                                                                                   gtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCAc                                        <  1:360980/139‑1 (MQ=255)
                                                                                    tggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACt                                       <  2:77512/139‑1 (MQ=255)
                                                                                      gtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAc                                     >  2:1019844/1‑139 (MQ=255)
                                                                                       ttAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                    <  1:184462/139‑1 (MQ=255)
                                                                                       ttAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                    <  1:907208/139‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGaa                                                                                                                          >  1:132054/1‑44 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGaa                                                                                                                          <  2:132054/44‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                  >  2:791891/1‑84 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                  <  1:791891/84‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCAc                                        >  2:487875/1‑126 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCAc                                        <  1:487875/126‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCAc                                        >  2:764711/1‑126 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCAc                                        <  1:764711/126‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                    >  1:665690/1‑130 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                    <  2:665690/130‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACGCTGTGTCAg                           >  2:715843/1‑139 (MQ=255)
                                                                                                   cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                         <  1:1051568/42‑1 (MQ=255)
                                                                                                   cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                         >  2:1051568/1‑42 (MQ=255)
                                                                                                     tCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaaa                      <  2:430934/139‑1 (MQ=255)
                                                                                                      cGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACGCTGTGTCAGGAAAAg                     <  1:715843/139‑1 (MQ=255)
                                                                                                        ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                          <  2:710505/52‑1 (MQ=255)
                                                                                                        ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                          >  1:710505/1‑52 (MQ=255)
                                                                                                            cTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGc                                                                                     <  2:531928/69‑1 (MQ=255)
                                                                                                            cTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGc                                                                                     >  1:531928/1‑69 (MQ=255)
                                                                                                                   aTCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATcc        <  2:897119/139‑1 (MQ=255)
                                                                                                                    tCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCt       >  2:834422/1‑139 (MQ=255)
                                                                                                                     cGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa             <  2:713505/132‑1 (MQ=255)
                                                                                                                     cGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa             >  1:713505/1‑132 (MQ=255)
                                                                                                                      gTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTa     <  1:1102341/139‑1 (MQ=255)
                                                                                                                      gTTCCTTTGTAACAGGTCCGAATCAAAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTa     <  2:90276/139‑1 (MQ=255)
                                                                                                                         ccTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAgg  >  1:591417/1‑139 (MQ=255)
                                                                                                                             tGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaa                       >  1:453333/1‑114 (MQ=255)
                                                                                                                             tGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaa                       <  2:453333/114‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                 
ATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGG  >  NC_000913/3248809‑3249068

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: