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.

ACTCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGGAGATCA  >  NC_000913/3248803‑3249074
                                                                                                                                         |                                                                                                                                       
aCTCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAg                                                                                                                                       >  1:1736275/1‑139 (MQ=255)
     cATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCg                                                                                                                                  <  2:115672/139‑1 (MQ=255)
     cATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCg                                                                                                                                  <  1:2489/139‑1 (MQ=255)
        cGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:1736275/139‑1 (MQ=255)
        cGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               <  1:804852/139‑1 (MQ=255)
         gCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATc                                                                                                                              <  1:1674818/139‑1 (MQ=255)
         gCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATc                                                                                                                              <  1:17878/139‑1 (MQ=255)
                 tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:909415/139‑1 (MQ=255)
                 tAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:1595177/139‑1 (MQ=255)
                    cTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATg                                                                                                                        <  1:141908/134‑1 (MQ=255)
                    cTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATg                                                                                                                        >  2:141908/1‑134 (MQ=255)
                       ccTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                                >  2:240997/1‑139 (MQ=255)
                           ttATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGa                                                                                                            <  1:240997/139‑1 (MQ=255)
                              tCCGAATCATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:655429/117‑1 (MQ=255)
                              tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               >  1:655429/1‑117 (MQ=255)
                              tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               >  1:193343/1‑117 (MQ=255)
                              tCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:193343/117‑1 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               <  2:1493386/107‑1 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               >  1:1493386/1‑107 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCa                                                                                               >  2:1315996/1‑139 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTCGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGaa                                                                                                           >  2:1155750/1‑127 (MQ=255)
                                        atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTCGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGaa                                                                                                           <  1:1155750/127‑1 (MQ=255)
                                            ctcactTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGc                                                                                           >  2:140148/1‑139 (MQ=255)
                                            ctcactTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGc                                                                                           >  1:1422711/1‑139 (MQ=255)
                                             tcactTTGGTGATTTCACCGTAACTGTCCATGATTAATGATGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACATGTCCGAATCCCAATGGAGGGTTCAAGAATGAAATACCGCATCGCt                                                                                          >  1:538632/1‑139 (MQ=255)
                                               actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                        >  1:892368/1‑139 (MQ=255)
                                               actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                        >  2:1448717/1‑139 (MQ=255)
                                                 tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      >  2:326710/1‑139 (MQ=255)
                                                 tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      >  1:1636053/1‑139 (MQ=255)
                                                   tGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATAc                                                                                                   <  2:316883/124‑1 (MQ=255)
                                                   tGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATAc                                                                                                   >  1:316883/1‑124 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:772282/1‑103 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:1450981/103‑1 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:1450981/1‑103 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:772282/103‑1 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGt                                                                                  <  2:1748832/139‑1 (MQ=255)
                                                     gTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGt                                                                                  >  1:438918/1‑139 (MQ=255)
                                                        aTTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGGCAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtc                                                                               <  1:1315996/139‑1 (MQ=255)
                                                         tttCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtct                                                                              >  1:458196/1‑139 (MQ=255)
                                                           tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtctct                                                                            >  2:1093423/1‑139 (MQ=255)
                                                            cACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtctctc                                                                           <  2:438918/139‑1 (MQ=255)
                                                              ccGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCtt                                                                         <  1:1755067/139‑1 (MQ=255)
                                                              ccGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCtt                                                                         <  1:326710/139‑1 (MQ=255)
                                                                         tATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTg                                                              >  1:206631/1‑139 (MQ=255)
                                                                             aTTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCgg                                                          <  1:1093423/139‑1 (MQ=255)
                                                                              ttAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGt                                                         <  2:1075739/139‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAg                                                       >  1:169132/1‑139 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGt                                                                                                                                     <  2:395210/60‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               >  2:185017/1‑66 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              >  1:924268/1‑83 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              <  2:555916/83‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              <  2:924268/83‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              >  1:555916/1‑83 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGttt                                                                                >  1:851399/1‑113 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGttt                                                                                <  2:851399/113‑1 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                      >  2:193778/1‑139 (MQ=255)
                                                                                 aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                      <  2:458196/139‑1 (MQ=255)
                                                                                     ggCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTg                                                                       <  2:3575/118‑1 (MQ=255)
                                                                                     ggCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTg                                                                       >  1:3575/1‑118 (MQ=255)
                                                                                      gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 <  1:193778/139‑1 (MQ=255)
                                                                                      gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 >  2:1817165/1‑139 (MQ=255)
                                                                                            gtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAc                                           <  2:1422711/139‑1 (MQ=255)
                                                                                             ttAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                          <  2:892368/139‑1 (MQ=255)
                                                                                               aaCACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCt                                        >  1:641124/1‑139 (MQ=255)
                                                                                                  acTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCtgtg                                     <  1:1817165/139‑1 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTcc                                                                                                                                   >  2:311800/1‑43 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTcc                                                                                                                                   <  1:311800/43‑1 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCtt                                                                  >  2:1068964/1‑108 (MQ=255)
                                                                                                    tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCtt                                                                  <  1:1068964/108‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTTTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:323017/1‑54 (MQ=255)
                                                                                                      ccGCTTCGTTGTTTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:323017/54‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAg                                                                                                                                       >  1:336587/1‑37 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      <  2:1737357/54‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                      >  1:1737357/1‑54 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                              >  2:1079299/1‑62 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 >  2:356695/1‑123 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGc                                                 <  1:356695/123‑1 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAg                                 >  2:1270708/1‑139 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAg                                 >  2:539275/1‑139 (MQ=255)
                                                                                                      ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCAAATCACAATGGAAGGTTCaa                                                                                                              <  1:1079299/62‑1 (MQ=255)
                                                                                                          ttCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaa                             <  2:641124/139‑1 (MQ=255)
                                                                                                              ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               <  1:1241255/37‑1 (MQ=255)
                                                                                                              ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAAt                                                                                                                               >  2:1241255/1‑37 (MQ=255)
                                                                                                              ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                         <  2:1636053/139‑1 (MQ=255)
                                                                                                                    ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGtt                                                                                                                 <  1:215708/45‑1 (MQ=255)
                                                                                                                    ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGtt                                                                                                                 >  2:215708/1‑45 (MQ=255)
                                                                                                                    ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa                   <  1:539275/139‑1 (MQ=255)
                                                                                                                    ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa                   <  1:1270708/139‑1 (MQ=255)
                                                                                                                      tcAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      >  1:1144376/2‑70 (MQ=255)
                                                                                                                      acgaTCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      <  2:1144376/67‑1 (MQ=255)
                                                                                                                        aaTCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                         <  2:684043/129‑1 (MQ=255)
                                                                                                                        aaTCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGa                         >  1:684043/1‑129 (MQ=255)
                                                                                                                          tCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCt             >  1:98831/1‑139 (MQ=255)
                                                                                                                           cGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                      <  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)
                                                                                                                                         |                                                                                                                                       
ACTCTCATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGGAGATCA  >  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: