Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,853,675 A→T 6.1% intergenic (+33/‑247) tisB → / → emrD toxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,853,6750AT 6.1% 689.6 / 21.9 210intergenic (+33/‑247)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (111/84);  new base T (5/8);  total (117/93)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

TGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTA  >  NC_000913/3853528‑3853822
                                                                                                                                                   |                                                                                                                                                   
tGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGa                                                                                                                                                     >  8:31845/1‑148 (MQ=255)
    ccTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCg                                                                                                                                                >  8:124879/1‑149 (MQ=255)
      tGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                   <  8:364545/144‑1 (MQ=255)
           cacaGGAGACGCGTATGAACCTGGTGGATATAGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCt                                                                                                                                         >  3:427411/1‑149 (MQ=255)
           cacaGGAGACGCGTATGAACCTGGTGGATATAGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCt                                                                                                                                         >  3:434137/1‑149 (MQ=255)
            acaGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCtt                                                                                                                                        >  2:76412/1‑149 (MQ=255)
            acaGGAGACGCGTATGAACCTGGTGGATATAGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCtt                                                                                                                                        <  4:434137/149‑1 (MQ=255)
             caGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCt                                                                                                                                         <  3:384435/147‑1 (MQ=255)
              aGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCtttg                                                                                                                                      <  7:77114/149‑2 (MQ=255)
              aGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCtttt                                                                                                                                      <  6:489041/149‑1 (MQ=255)
                   aCGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATAtt                                                                                                                                 <  5:116785/149‑1 (MQ=255)
                     gcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTcc                                                                                                                                            >  4:96199/1‑136 (MQ=255)
                     gcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATAtt                                                                                                                                 >  1:402077/1‑147 (MQ=255)
                        tATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCg                                                                                                                            >  2:269213/1‑149 (MQ=255)
                         aTGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCtt                                                                                                                                        >  4:48516/1‑136 (MQ=255)
                          tGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTg                                                                                                                          >  4:122176/1‑149 (MQ=255)
                           gAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGc                                                                                                                         >  6:202240/1‑149 (MQ=255)
                             aCCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCgg                                                                                                                       >  1:392167/1‑149 (MQ=255)
                              ccTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCggg                                                                                                                      >  4:103423/1‑149 (MQ=255)
                                  gtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTgg                                                                                                                  >  2:149076/1‑149 (MQ=255)
                                  gtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTgg                                                                                                                  <  6:165806/149‑1 (MQ=255)
                                  gtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTgg                                                                                                                  >  8:103836/1‑149 (MQ=255)
                                  gtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTg                                                                                                                   <  7:124879/148‑1 (MQ=255)
                                       tatCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCg                                                                                                             <  8:482/149‑1 (MQ=255)
                                         tCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGga                                                                                                           <  1:293052/149‑1 (MQ=255)
                                           gCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgat                                                                                                          <  2:212594/148‑1 (MQ=255)
                                                 cgatctcTCATACTCATTGTTTCATGAGTGCATCTGGTTGATGCTGTTCTGATATACCTGAAGTAATTCAGATTCTAGTCGCACCTAAGGGTTGCGGGTAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgcgatgc                                                                                                   <  7:179432/143‑1 (MQ=255)
                                                   tATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCgg                                                                                                                       <  2:259892/127‑1 (MQ=255)
                                                   tATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCgg                                                                                                                       >  1:259892/1‑127 (MQ=255)
                                                     tCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAATGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGa                                                                                                <  8:53600/148‑1 (MQ=255)
                                                       cTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAGACCGCTCCCCTTTATATTAAGCGTGCCGGTTGGTGTCGGATTCGATGCTGACGct                                                                                            >  1:206147/1‑148 (MQ=255)
                                                        tCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGCTTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCa                                                                                            >  3:30626/1‑149 (MQ=255)
                                                        tCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCa                                                                                            >  7:227371/1‑149 (MQ=255)
                                                        tCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCa                                                                                            >  6:59662/1‑149 (MQ=255)
                                                        tCAAACACATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAACGGGAGCGGGAAACCGCTCCCCTTATATATTTAGCGTCCCGGTTGGTGTCGgatgcgatgc                                                                                                   >  8:179432/1‑142 (MQ=255)
                                                          aaaCTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATc                                                                                          <  5:202240/149‑1 (MQ=255)
                                                          aaaCTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATc                                                                                          <  7:31845/149‑1 (MQ=255)
                                                          aaaCTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATc                                                                                          <  4:144960/149‑1 (MQ=255)
                                                           aaCTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCt                                                                                         >  2:161548/1‑149 (MQ=255)
                                                             cTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGAc                                                                                               <  7:209126/141‑1 (MQ=255)
                                                             cTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGAc                                                                                               >  8:209126/1‑141 (MQ=255)
                                                                aTTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGGTGGTGTCGGATGCGATGCTGACGCATCTTATc                                                                                     <  7:237500/148‑1 (MQ=255)
                                                                   gttgCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGc                                                                                  <  5:317269/148‑1 (MQ=255)
                                                                   gttgCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGc                                                                                  <  4:252213/148‑1 (MQ=255)
                                                                       cagcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTAc                                                                             >  2:185376/1‑149 (MQ=255)
                                                                        agcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATc                                                                                     <  1:435879/140‑1 (MQ=255)
                                                                        agcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATc                                                                                     >  2:435879/1‑140 (MQ=255)
                                                                        agcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTAc                                                                             >  4:427188/1‑148 (MQ=255)
                                                                             tGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCGTCTTATCCGCCCTACCAtctc                                                                       >  1:103717/1‑149 (MQ=255)
                                                                             tGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGAGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtct                                                                        >  7:246772/1‑148 (MQ=255)
                                                                               cAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtctctc                                                                     >  3:125539/1‑149 (MQ=255)
                                                                                aaCTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTcc                                                                    >  8:379091/1‑149 (MQ=255)
                                                                                     cTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGc                                                                <  3:92516/148‑1 (MQ=255)
                                                                                      ttGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCaa                                                              >  4:80627/1‑149 (MQ=255)
                                                                                      ttGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCaa                                                              >  3:11367/1‑149 (MQ=255)
                                                                                         aTGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACAt                                                           >  8:253573/1‑149 (MQ=255)
                                                                                           gCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACAtt                                                          >  1:26841/1‑148 (MQ=255)
                                                                                            cTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACAttt                                                         <  2:402077/148‑1 (MQ=255)
                                                                                            cTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTa                                                        <  5:59662/149‑1 (MQ=255)
                                                                                            cTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTa                                                        <  2:392167/149‑1 (MQ=255)
                                                                                               ttCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTg                                                     >  7:220138/1‑149 (MQ=255)
                                                                                                 cTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGcc                                                   >  4:52132/1‑149 (MQ=255)
                                                                                                         ccTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGt                                           <  6:286171/149‑1 (MQ=255)
                                                                                                          cTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGtt                                          >  5:3549/1‑149 (MQ=255)
                                                                                                           tGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGtt                                          >  8:237587/1‑148 (MQ=255)
                                                                                                           tGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGtt                                          >  8:213691/1‑148 (MQ=255)
                                                                                                           tGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGtt                                          >  1:383483/1‑148 (MQ=255)
                                                                                                            gAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCgg                                                                                                                       >  3:157984/1‑70 (MQ=255)
                                                                                                             aaGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTa                                        <  1:149076/148‑1 (MQ=255)
                                                                                                             aaGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTa                                        <  5:196029/148‑1 (MQ=255)
                                                                                                                  aTTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACGATCTCTCACGGCAACATTTATTGCCGCTTTTGTTTACATAtt                                  <  8:246772/149‑1 (MQ=255)
                                                                                                                  aTTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATAtt                                  <  8:401137/149‑1 (MQ=255)
                                                                                                                   ttCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTg                                            >  6:238423/1‑138 (MQ=255)
                                                                                                                   ttCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACata                                    >  5:94397/1‑146 (MQ=255)
                                                                                                                   ttCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACata                                    <  6:94397/146‑1 (MQ=255)
                                                                                                                   ttCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTc                                 >  5:151737/1‑149 (MQ=255)
                                                                                                                   ttCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTc                                 >  2:311514/1‑149 (MQ=255)
                                                                                                                      aGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGc                              <  1:154945/149‑1 (MQ=255)
                                                                                                                      aGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGc                              >  2:463636/1‑149 (MQ=255)
                                                                                                                      aGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGc                              <  5:274421/149‑1 (MQ=255)
                                                                                                                        aTTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCg                            >  3:239216/1‑149 (MQ=255)
                                                                                                                         ttCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGc                           <  1:18022/149‑1 (MQ=255)
                                                                                                                          tCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGc                           >  1:471673/1‑148 (MQ=255)
                                                                                                                            aaGTCGCTCCAAAGGGGATCTGGAAACCGCTCCCCTTTATATTTAGCGTGCGGGTTTGTGTGGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAAGATTTATTGCCGCTTTTGTTTACATATTCTGGCGCTaaa                       <  2:206147/149‑1 (MQ=255)
                                                                                                                             aGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTaaa                       <  6:350159/149‑1 (MQ=255)
                                                                                                                             aGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTaaa                       <  3:52132/149‑1 (MQ=255)
                                                                                                                              gTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGttt                                         <  6:367333/130‑1 (MQ=255)
                                                                                                                              gTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGttt                                         >  5:367333/1‑130 (MQ=255)
                                                                                                                                  tataaaaGGGGAGCGGTTTCCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGcgg                            >  1:494434/5‑137 (MQ=255)
                                                                                                                                   aCCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTc                 <  8:152775/149‑1 (MQ=255)
                                                                                                                                    taaaaGGGGAGCGGTTTCCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGc                              <  7:409165/133‑1 (MQ=255)
                                                                                                                                    taaaaGGGGAGCGGTTTCCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGc                              >  8:409165/3‑135 (MQ=255)
                                                                                                                                     cAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTccc               >  7:451974/1‑149 (MQ=255)
                                                                                                                                        aGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCAt            >  1:230055/1‑149 (MQ=255)
                                                                                                                                         ggggAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAATAAATGTTGCCggga                                               <  4:189571/113‑4 (MQ=38)
                                                                                                                                         ggggAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAATAAATGTTGCCggga                                               >  3:189571/1‑110 (MQ=38)
                                                                                                                                            gAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCt        >  4:337326/1‑149 (MQ=255)
                                                                                                                                            gAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCt        >  7:63762/1‑149 (MQ=255)
                                                                                                                                               cGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGc     <  2:383483/149‑1 (MQ=255)
                                                                                                                                                gggAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTc                 >  8:438118/1‑136 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTA  >  NC_000913/3853528‑3853822

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