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,9390TC98.2% 382.0 / ‑3.4 108intergenic (+117/‑30)mzrA/yqjCmodulator of EnvZ/OmpR regulon/DUF1090 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (53/53);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

Only 100 of 108 total aligned reads displayed.

ACTCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGGAGATCA  >  NC_000913/3248803‑3249074
                                                                                                                                         |                                                                                                                                       
aCTCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAg                                                                                                                                       >  1:1736275/1‑139 (MQ=255)
     cATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCg                                                                                                                                  <  2:115672/139‑1 (MQ=255)
     cATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCg                                                                                                                                  <  1:2489/139‑1 (MQ=255)
        cGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:1736275/139‑1 (MQ=255)
        cGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               <  1:804852/139‑1 (MQ=255)
         gCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATc                                                                                                                              <  1:1674818/139‑1 (MQ=255)
         gCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATc                                                                                                                              <  1:17878/139‑1 (MQ=255)
                 tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:909415/139‑1 (MQ=255)
                 tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:1595177/139‑1 (MQ=255)
                    cTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATg                                                                                                                        <  1:141908/134‑1 (MQ=255)
                    cTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATg                                                                                                                        >  2:141908/1‑134 (MQ=255)
                       ccTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                                >  2:240997/1‑139 (MQ=255)
                           ttATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGa                                                                                                            <  1:240997/139‑1 (MQ=255)
                              tCCGAATCATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:655429/117‑1 (MQ=255)
                              tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               >  1:655429/1‑117 (MQ=255)
                              tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               >  1:193343/1‑117 (MQ=255)
                              tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:193343/117‑1 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:1493386/107‑1 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               >  1:1493386/1‑107 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCa                                                                                               >  2:1315996/1‑139 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTCGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGaa                                                                                                           >  2:1155750/1‑127 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTCGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGaa                                                                                                           <  1:1155750/127‑1 (MQ=255)
                                            ctcactTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGc                                                                                           >  2:140148/1‑139 (MQ=255)
                                            ctcactTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGc                                                                                           >  1:1422711/1‑139 (MQ=255)
                                             tcactTTGGTGATTTCACCGTAACTGTCCATGATTAATGATGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACATGTCCGAATCCCAATGGAGGGTTCAAGAATGAAATACCGCATCGCt                                                                                          >  1:538632/1‑139 (MQ=255)
                                               actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                        >  1:892368/1‑139 (MQ=255)
                                               actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                        >  2:1448717/1‑139 (MQ=255)
                                                 tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      >  2:326710/1‑139 (MQ=255)
                                                 tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      >  1:1636053/1‑139 (MQ=255)
                                                   tGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATAc                                                                                                   <  2:316883/124‑1 (MQ=255)
                                                   tGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATAc                                                                                                   >  1:316883/1‑124 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:772282/1‑103 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:1450981/103‑1 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:1450981/1‑103 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:772282/103‑1 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGt                                                                                  <  2:1748832/139‑1 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGt                                                                                  >  1:438918/1‑139 (MQ=255)
                                                        aTTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGGCAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtc                                                                               <  1:1315996/139‑1 (MQ=255)
                                                         tttCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtct                                                                              >  1:458196/1‑139 (MQ=255)
                                                           tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtctct                                                                            >  2:1093423/1‑139 (MQ=255)
                                                            cACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtctctc                                                                           <  2:438918/139‑1 (MQ=255)
                                                              ccGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCtt                                                                         <  1:1755067/139‑1 (MQ=255)
                                                              ccGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCtt                                                                         <  1:326710/139‑1 (MQ=255)
                                                                         tATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTg                                                              >  1:206631/1‑139 (MQ=255)
                                                                             aTTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCgg                                                          <  1:1093423/139‑1 (MQ=255)
                                                                              ttAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGt                                                         <  2:1075739/139‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAg                                                       >  1:169132/1‑139 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGt                                                                                                                                     <  2:395210/60‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               >  2:185017/1‑66 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              >  1:924268/1‑83 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              <  2:555916/83‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              <  2:924268/83‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              >  1:555916/1‑83 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGttt                                                                                >  1:851399/1‑113 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGttt                                                                                <  2:851399/113‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                      >  2:193778/1‑139 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                      <  2:458196/139‑1 (MQ=255)
                                                                                     ggCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTg                                                                       <  2:3575/118‑1 (MQ=255)
                                                                                     ggCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTg                                                                       >  1:3575/1‑118 (MQ=255)
                                                                                      gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 <  1:193778/139‑1 (MQ=255)
                                                                                      gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 >  2:1817165/1‑139 (MQ=255)
                                                                                            gtTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAc                                           <  2:1422711/139‑1 (MQ=255)
                                                                                             ttAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                          <  2:892368/139‑1 (MQ=255)
                                                                                               aaCACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCt                                        >  1:641124/1‑139 (MQ=255)
                                                                                                  acTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCtgtg                                     <  1:1817165/139‑1 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTcc                                                                                                                                   >  2:311800/1‑43 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTcc                                                                                                                                   <  1:311800/43‑1 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCtt                                                                  >  2:1068964/1‑108 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCtt                                                                  <  1:1068964/108‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTTTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:323017/1‑54 (MQ=255)
                                                                                                      ccGCTTCGTTGTTTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:323017/54‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAg                                                                                                                                       >  1:336587/1‑37 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:1737357/54‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:1737357/1‑54 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              >  2:1079299/1‑62 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 >  2:356695/1‑123 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 <  1:356695/123‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAg                                 >  2:1270708/1‑139 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAg                                 >  2:539275/1‑139 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCAAATCACAATGGAAGGTTCaa                                                                                                              <  1:1079299/62‑1 (MQ=255)
                                                                                                          ttCGTTGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaa                             <  2:641124/139‑1 (MQ=255)
                                                                                                              ttGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               <  1:1241255/37‑1 (MQ=255)
                                                                                                              ttGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAAt                                                                                                                               >  2:1241255/1‑37 (MQ=255)
                                                                                                              ttGTCTGGACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                         <  2:1636053/139‑1 (MQ=255)
                                                                                                                    ggACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGtt                                                                                                                 <  1:215708/45‑1 (MQ=255)
                                                                                                                    ggACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGtt                                                                                                                 >  2:215708/1‑45 (MQ=255)
                                                                                                                    ggACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa                   <  1:539275/139‑1 (MQ=255)
                                                                                                                    ggACAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa                   <  1:1270708/139‑1 (MQ=255)
                                                                                                                      tcAATCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      >  1:1144376/2‑70 (MQ=255)
                                                                                                                      acgaTCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      <  2:1144376/67‑1 (MQ=255)
                                                                                                                        aaTCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                         <  2:684043/129‑1 (MQ=255)
                                                                                                                        aaTCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                         >  1:684043/1‑129 (MQ=255)
                                                                                                                          tCGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCt             >  1:98831/1‑139 (MQ=255)
                                                                                                                           cGTT‑CCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      <  1:1364127/65‑1 (MQ=255)
                                                                                                                                ccTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAg                                                                <  1:1345964/83‑1 (MQ=255)
                                                                                                                                ccTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAg                                                                >  2:1345964/1‑83 (MQ=255)
                                                                                                                                ccTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTa                                                    <  1:1651675/95‑1 (MQ=255)
                                                                                                                                ccTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTa                                                    >  2:1651675/1‑95 (MQ=255)
                                                                                                                                   ttGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGgaga     >  1:1672790/1‑139 (MQ=255)
                                                                                                                                      tAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGGAGATCa  >  1:1117450/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                       
ACTCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTT‑CCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGGAGATCA  >  NC_000913/3248803‑3249074

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