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,9390TC97.9% 327.7 / ‑3.2 93intergenic (+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 (48/43);  total (49/44)
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.43e-01

ATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGGAG  >  NC_000913/3248809‑3249070
                                                                                                                                  |                                                                                                                                   
aTCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGa                                                                                                                             >  1:1096495/1‑139 (MQ=255)
     ttCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATcac                                                                                                                        >  2:212081/1‑139 (MQ=255)
         gATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATg                                                                                                                    <  2:363508/139‑1 (MQ=255)
                tCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGtt                                                                                                             <  1:159892/139‑1 (MQ=255)
                 ccTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                            <  1:224002/139‑1 (MQ=255)
                   tGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  >  1:828284/1‑131 (MQ=255)
                   tGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  <  2:828284/131‑1 (MQ=255)
                                  atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGa                                                                                                                  <  2:1108875/116‑1 (MQ=255)
                                  atCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGa                                                                                                                  >  1:1108875/1‑116 (MQ=255)
                                     actcactTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCg                                                                                        >  2:665752/1‑139 (MQ=255)
                                         actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          >  1:51638/1‑117 (MQ=255)
                                         actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          <  2:51638/117‑1 (MQ=255)
                                         actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                    >  1:1229685/1‑139 (MQ=255)
                                         actTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCttt                                                                                    >  2:905257/1‑139 (MQ=255)
                                           tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                  >  1:385304/1‑139 (MQ=255)
                                           tttGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                  <  2:385304/139‑1 (MQ=255)
                                              ggTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGt                                                                                                                                 <  1:384700/89‑1 (MQ=255)
                                              ggTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGt                                                                                                                                 >  2:384700/1‑89 (MQ=255)
                                               gTGATTTCACCGTAACTGTCTATGATTAATGCGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTTCTTTGTAACAGGACCGAATCTCAATGGAAGGTTCAAGAATGTAATACCGAATCGCTTTAGCTGt                                                                              >  2:604500/1‑139 (MQ=255)
                                                   tttCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTGACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTtct                                                                          >  2:817942/1‑139 (MQ=255)
                                                    ttCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  <  1:1083020/98‑1 (MQ=255)
                                                    ttCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  >  2:1083020/1‑98 (MQ=255)
                                                     tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCGTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGc                                                                                            >  2:523091/1‑119 (MQ=255)
                                                     tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCGTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGc                                                                                            <  1:523091/119‑1 (MQ=255)
                                                     tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  <  1:327959/97‑1 (MQ=255)
                                                     tCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  >  2:327959/1‑97 (MQ=255)
                                                          gTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTg                                                                   <  1:1139829/139‑1 (MQ=255)
                                                                 tCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAg                                                            >  2:519219/1‑139 (MQ=255)
                                                                   tATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATca                                                                                                                         <  2:1174456/76‑1 (MQ=255)
                                                                   tATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATca                                                                                                                         >  1:1174456/1‑76 (MQ=255)
                                                                   tATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                            <  1:314924/89‑1 (MQ=255)
                                                                   tATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                            >  2:314924/1‑89 (MQ=255)
                                                                       aTTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCgg                                                      >  1:1130071/1‑139 (MQ=255)
                                                                         tAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTa                                                    >  1:466622/1‑139 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          <  1:1259447/83‑1 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          >  2:1259447/1‑83 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCa                                                                                           >  2:834183/1‑98 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCa                                                                                           <  1:834183/98‑1 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                  >  2:121670/1‑139 (MQ=255)
                                                                           aTGAGGCGGTGGTTAACACTACCGCTTCGTTGTATGGACAATCGTTCCTTTGTAACAGGACCGAATCAAATTGGAAGGTTAAAGACTGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                  >  2:689565/1‑139 (MQ=255)
                                                                                gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGaa                                                                                                       >  2:230184/1‑81 (MQ=255)
                                                                                gCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGaa                                                                                                       <  1:230184/81‑1 (MQ=255)
                                                                                 cGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGcc                                            <  2:1229685/139‑1 (MQ=255)
                                                                                  ggtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  >  2:802322/1‑68 (MQ=255)
                                                                                  ggtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          >  1:653617/1‑76 (MQ=255)
                                                                                  ggtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          <  2:653617/76‑1 (MQ=255)
                                                                                  ggtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTACTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  <  1:802322/68‑1 (MQ=255)
                                                                                   gtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  >  2:444680/1‑67 (MQ=255)
                                                                                   gtggtTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGa                                                                                                                  <  1:444680/67‑1 (MQ=255)
                                                                                       ttAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                      <  2:895575/139‑1 (MQ=255)
                                                                                         aaCACTACCGCTTCGTTGTCTGGACAATCGATCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCt                                    <  1:604500/139‑1 (MQ=255)
                                                                                              tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATAc                                                                                               >  2:58/1‑75 (MQ=255)
                                                                                              tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATAc                                                                                               <  1:58/75‑1 (MQ=255)
                                                                                              tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTc                               >  2:269360/1‑139 (MQ=255)
                                                                                              tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTc                               >  1:279206/1‑139 (MQ=255)
                                                                                              tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTc                               >  1:334722/1‑139 (MQ=255)
                                                                                              tACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTc                               >  1:466327/1‑139 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTc                                                                                                                                <  1:194254/40‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTc                                                                                                                                >  2:194254/1‑40 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATca                                                                                                                         <  2:97519/47‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATca                                                                                                                         >  1:97519/1‑47 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTTGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                      >  1:369439/1‑130 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTTGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTAcc                                      <  2:369439/130‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                  >  2:310604/1‑86 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAg                                                                                  <  1:310604/86‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                  <  1:949139/118‑1 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGt                                                  >  2:949139/1‑118 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCt                                    >  1:953923/1‑132 (MQ=255)
                                                                                                ccGCTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCt                                    <  2:953923/132‑1 (MQ=255)
                                                                                                   cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAgg                                                                                                               >  1:477529/1‑54 (MQ=255)
                                                                                                   cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAgg                                                                                                               <  2:477529/54‑1 (MQ=255)
                                                                                                   cTTCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaa                          >  1:1085962/1‑139 (MQ=255)
                                                                                                     tCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                            >  2:784191/1‑55 (MQ=255)
                                                                                                     tCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                            <  1:784191/55‑1 (MQ=255)
                                                                                                     tCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaaa                        >  1:1205632/1‑139 (MQ=255)
                                                                                                     tCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaaa                        >  1:1079015/1‑139 (MQ=255)
                                                                                                     tCGTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGaaaa                        >  2:86457/1‑139 (MQ=255)
                                                                                                       gTTGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAgg                      <  2:1085962/139‑1 (MQ=255)
                                                                                                        ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          <  1:1142377/54‑1 (MQ=255)
                                                                                                        ttGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCaa                                                                                                          >  2:1142377/1‑54 (MQ=255)
                                                                                                         tGTCTGGACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAg                    <  2:466327/139‑1 (MQ=255)
                                                                                                              ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa               <  2:1205632/139‑1 (MQ=255)
                                                                                                              ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa               <  1:269360/139‑1 (MQ=255)
                                                                                                              ggACAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCaaaa               <  1:121670/139‑1 (MQ=255)
                                                                                                                 cAATCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAatat            <  2:285827/139‑1 (MQ=255)
                                                                                                                   aTCGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATcc          <  1:665752/139‑1 (MQ=255)
                                                                                                                     cGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCtt                                                                                     >  1:465809/1‑62 (MQ=255)
                                                                                                                     cGTTCCTTTGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCtt                                                                                     <  2:465809/62‑1 (MQ=255)
                                                                                                                         ccTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                            >  1:267342/1‑35 (MQ=255)
                                                                                                                         ccTTTGTAACAGGTCCGAATCACAATGGAAGGTTc                                                                                                            <  2:267342/35‑1 (MQ=255)
                                                                                                                           tttGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGgag  <  1:86457/139‑1 (MQ=255)
                                                                                                                            ttGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAat              >  1:526075/1‑126 (MQ=255)
                                                                                                                            ttGTAACAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAat              <  2:526075/126‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
ATCGCTTCGGATAACTTCCTGTTATCCGAAACATATCACTCACTTTGGTGATTTCACCGTAACTGTCTATGATTAATGAGGCGGTGGTTAACACTACCGCTTCGTTGTCTGGACAATCGTTCCTTTGTAATAGGTCCGAATCACAATGGAAGGTTCAAGAATGAAATACCGCATCGCTTTAGCTGTTTCTCTCTTTGCTCTTAGTGCCGGTAGTTATGCCACTACCCTGTGTCAGGAAAAGGAGCAAAATATCCTTAAGGAG  >  NC_000913/3248809‑3249070

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: