Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,227 C→G 18.1% intergenic (+20/+11) marB → / ← eamA periplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,620,2270CG18.1% 262.8 / 56.3 113intergenic (+20/+11)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (33/58);  new base G (10/10);  total (44/69)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

Only 100 of 117 total aligned reads displayed.

TGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATT  >  NC_000913/1620129‑1620326
                                                                                                  |                                                                                                   
 gATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtct                                                                                                    >  1:1813398/1‑99 (MQ=255)
    ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCggg                                                                                                >  1:2142746/1‑100 (MQ=255)
    ttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCggg                                                                                                <  1:3873435/100‑1 (MQ=255)
     tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCgggg                                                                                               <  1:2603531/100‑1 (MQ=255)
      aaccaacaTGGGTAATGGCAGTGATAAGTAGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTTTTCTGGCCCCGACATCTCGGGGc                                                                                              <  2:2945382/95‑1 (MQ=255)
      aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                              <  2:3811151/100‑1 (MQ=255)
       aaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCt                                                                                             <  1:2937174/100‑1 (MQ=255)
        aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGCGATGTCGGGGcca                                                                                            <  2:2413297/100‑3 (MQ=255)
        aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTAGAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCt                                                                                             <  2:2352770/99‑1 (MQ=255)
         tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                           <  2:40727/100‑4 (MQ=255)
          cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccaga                                                                                          <  2:1506540/100‑5 (MQ=255)
          cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccaga                                                                                          <  2:1371515/100‑5 (MQ=255)
          cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                          >  2:418152/1‑100 (MQ=255)
           acaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                         <  2:2412460/100‑1 (MQ=255)
            caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCTTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGCCCCGAGATGTCGGGGccagaaca                                                                                       <  1:3258477/100‑8 (MQ=255)
            caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccagaac                                                                                        <  1:2602124/100‑7 (MQ=255)
            caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccagaac                                                                                        <  1:243825/100‑7 (MQ=255)
            caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                        <  2:3035179/100‑1 (MQ=255)
            caTGGGTACTGGCAGTGATAAGTCGGATGCGCTAGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                         >  2:2585936/1‑99 (MQ=255)
             aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTaa                                                                                       >  1:2471282/1‑100 (MQ=255)
              tGGGTACTGGCAGTGATAAGTCGTATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                      <  1:1534337/100‑1 (MQ=255)
              tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                      <  1:1632971/100‑1 (MQ=255)
               gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                     <  2:2851725/100‑1 (MQ=255)
               gagTACTGGCAGTGATAAGTCGGATGCACTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                     <  1:2548648/98‑1 (MQ=255)
                 gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTaa                                                                                       >  2:2343952/1‑96 (MQ=255)
                  tACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                  <  1:3830955/100‑1 (MQ=255)
                   aCTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                                 <  1:1163238/100‑1 (MQ=255)
                    cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                <  1:2352840/100‑1 (MQ=255)
                    cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                >  1:600840/1‑100 (MQ=255)
                      ggCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAcc                                                                              <  2:3563361/100‑1 (MQ=255)
                      ggCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAcc                                                                              >  1:2714755/1‑100 (MQ=255)
                            gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcn                                                                        <  1:2093176/100‑2 (MQ=255)
                             aTAAGTCGGATGCGCTCGACGTGCCCTATTATAATCAACACGCTATTTAGTTTGTTCTGGCCCTGACATCTCGGTTCTTATTAACTTCCCACCttaa                                                                          <  2:1754257/97‑3 (MQ=255)
                              tAAGTCGGATGCGCTCGGCGTGCCTTATTATAATCAACAGGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGc                                                                      <  1:677336/100‑1 (MQ=255)
                              tAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGc                                                                      >  2:1371681/1‑100 (MQ=255)
                               aaGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCt                                                                     <  1:669168/100‑1 (MQ=255)
                                aGTCGGATGCGCTCGGCGTGACCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATATCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                    >  2:636929/1‑100 (MQ=255)
                                  tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTa                                                                  <  1:3136206/100‑1 (MQ=255)
                                   cGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGTCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                 <  2:1924364/100‑1 (MQ=255)
                                   cGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                                 >  1:247635/1‑100 (MQ=255)
                                    ggATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                                <  1:2002222/100‑1 (MQ=255)
                                     gATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGc                                                               >  1:2815051/1‑100 (MQ=255)
                                      aTGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGcc                                                              >  2:1357367/1‑100 (MQ=255)
                                         cgcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCAcc                                                           <  2:1155167/100‑1 (MQ=255)
                                          gcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCg                                                          <  2:2714731/100‑1 (MQ=255)
                                          gcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCg                                                          >  2:2793036/1‑100 (MQ=255)
                                               gCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGc                                                     <  2:641884/100‑1 (MQ=255)
                                                cGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                    <  1:3703181/100‑1 (MQ=255)
                                                cGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACACCACCGCAAGcc                                                    <  2:1770354/100‑1 (MQ=255)
                                                   gCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGATTTACGCCACCGCAAGCCaaa                                                 >  1:722798/1‑100 (MQ=255)
                                                    cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                                >  1:1110623/1‑100 (MQ=255)
                                                    cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                                <  1:2158649/100‑1 (MQ=255)
                                                    cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                                >  2:1577621/1‑100 (MQ=255)
                                                    cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                                <  1:2791152/100‑1 (MQ=255)
                                                    cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACAGCTTTACGCCACCGCAAGCCAAAt                                                >  2:3200979/1‑100 (MQ=255)
                                                     ccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                               <  1:3177418/100‑1 (MQ=255)
                                                      ctactaTAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAGCTCGAGGCTGATGAACTTGCCACCTTTACCGCTTTACGCCACCGCAAGCTAAATa                                               <  2:1818562/99‑1 (MQ=255)
                                                      cTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGTCATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                               <  1:694405/99‑1 (MQ=255)
                                                      cTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATAc                                              <  2:124942/100‑1 (MQ=255)
                                                      cTATAATAATCAACACGCTATGTAGATTGATCTAGCCCCGACATCTCGCGGCTTATTAACCTCCCACCTTCACCGCTTTAAGCCCCCGCAAACCAAATAc                                              >  2:2344870/1‑100 (MQ=255)
                                                         ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                           <  1:3807281/100‑1 (MQ=255)
                                                         ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                           <  1:149578/100‑1 (MQ=255)
                                                         ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                           >  1:1830596/1‑100 (MQ=255)
                                                          tataATGAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                           <  2:393293/99‑1 (MQ=255)
                                                          tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGTCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                          <  1:161036/100‑1 (MQ=255)
                                                          tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                          >  1:3666478/1‑100 (MQ=255)
                                                            taATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGat                                        <  1:1253358/100‑1 (MQ=255)
                                                                  acacCCTATGTAGTTTGTTCTGGCCCCGACATCTAGCGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAg                                  <  2:1913569/100‑1 (MQ=255)
                                                                    acGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                                 >  1:1429532/1‑99 (MQ=255)
                                                                       cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                             <  2:2790422/100‑1 (MQ=255)
                                                                       cTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                             <  1:2844051/100‑1 (MQ=255)
                                                                         aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                             >  1:1389764/1‑98 (MQ=255)
                                                                         aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                           >  1:2435015/1‑100 (MQ=255)
                                                                          tGTAGTTTGTTCTGGCCCCGACATCTCGGTGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCa                          <  1:3645516/100‑1 (MQ=255)
                                                                           gTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCAt                         >  1:1577118/1‑100 (MQ=255)
                                                                            tAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAACCAAATACATTGATATACAGCCCGGTCATa                        >  1:1606292/1‑100 (MQ=255)
                                                                             aGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATaa                       <  1:3259889/100‑1 (MQ=255)
                                                                              gTTTGTTCTGGCCCCGACATCTCGGGGCTTATTGACTTCCCACCTTTACCGCTTTACGCCTCCGCAAGCCAAATACTTTGATATACAGCCCGGTCATAAt                      <  2:3344987/100‑1 (MQ=255)
                                                                              gTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAAt                      >  1:3141981/1‑100 (MQ=255)
                                                                                gttaataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAACCAAATACATTAATATACAGCCCGGTCATAATGa                    >  1:622785/9‑100 (MQ=38)
                                                                                 ttaataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAg                   >  1:3910724/8‑100 (MQ=255)
                                                                                 tGTTCTGGCCCCTACCTCTCGGTGCTTATTAACTTGCCACCTTGAACGCTTTACGCCACCGCAATCCAAATACATTTATATACAGCCCGGTCATAAt                      <  2:3506098/97‑1 (MQ=255)
                                                                                  naataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCAAAATACATTGATATACAGCCCGGTCATAATGAg                   >  1:3483647/7‑99 (MQ=255)
                                                                                  gTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCCCTTTACACCACCGCAAGCCAAATACATAAATATACAGCACGGTCATAATAAg                   >  1:862137/1‑99 (MQ=255)
                                                                                   ttCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca                 <  2:3183841/100‑1 (MQ=255)
                                                                                    ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAAACATAGATATACAGCCCGGTCATAATGAgcac                >  1:2141615/5‑100 (MQ=255)
                                                                                     cTGGCCCCGACATCTCGTGGCTTATTAACATCCCACCTTTACCGCTTTACGCCACCGCAAGCAAAATACATTGTGATACAGCCCGGTCATAATGAgcacc               >  2:2437275/1‑100 (MQ=255)
                                                                                      aagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccg              >  1:3888370/3‑100 (MQ=255)
                                                                                       ggCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc             <  1:1450587/100‑1 (MQ=255)
                                                                                       agCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc             >  1:1169366/2‑100 (MQ=255)
                                                                                       agCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc             >  2:2458034/2‑100 (MQ=255)
                                                                                       agCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCACTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc             >  1:2331958/2‑100 (MQ=255)
                                                                                        gCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgca            <  1:1122070/100‑1 (MQ=255)
                                                                                         ccccGACATCTCGGGGCTTATTAACTTTCCTCCTTTTCCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac           <  1:1219478/100‑1 (MQ=255)
                                                                                         ccccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcac           >  2:2325299/1‑100 (MQ=255)
                                                                                          cccGAGATGTCGGGGCCTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcacc          >  2:2018315/1‑100 (MQ=255)
                                                                                          cccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcacc          <  2:1606277/100‑1 (MQ=255)
                                                                                             gACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTaa       >  1:2240296/1‑100 (MQ=255)
                                                                                                  ctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAtt  >  1:2240917/1‑100 (MQ=255)
                                                                                                  ctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAtt  >  1:2696872/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
TGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATT  >  NC_000913/1620129‑1620326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: