Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,224 C→G 16.9% intergenic (+17/+14) 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,2240CG16.9% 294.3 / 52.0 118intergenic (+17/+14)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (39/59);  new base G (10/10);  total (49/69)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-02

Only 100 of 120 total aligned reads displayed.

CGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAA  >  NC_000913/1620125‑1620321
                                                                                                   |                                                                                                 
cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAc                                                                                                   <  1:3783527/100‑1 (MQ=255)
cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAc                                                                                                   <  1:3329684/100‑1 (MQ=255)
 gTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACa                                                                                                  >  1:3220531/1‑100 (MQ=255)
    tGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCATATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACAtct                                                                                               >  1:2034728/1‑100 (MQ=255)
     gctttAAATCACATGGGTACTGGCAGTGATAAGTCGGATACGCTCGGAGTGCCCTATTATAATAAACACGCTAAGTAGTTAGTACTGGCACAGACa                                                                                                  >  2:2112456/3‑96 (MQ=255)
         tAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCgggg                                                                                          <  1:2603531/100‑1 (MQ=255)
          aaccaacaTGGGTAATGGCAGTGATAAGTAGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTTTTCTGGCCCCGACATCTCGGGGc                                                                                         <  2:2945382/95‑1 (MQ=255)
          aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                         <  2:3811151/100‑1 (MQ=255)
          aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCACCGACATCTCGGGGc                                                                                         <  1:1297114/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:1371515/100‑5 (MQ=255)
              cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccaga                                                                                     <  2:1372085/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: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)
                  tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                 <  1:1632971/100‑1 (MQ=255)
                   gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                <  2:2851725/100‑1 (MQ=255)
                   gagTACTGGCAGTGATAAGTCGGATGCACTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                <  1:2548648/98‑1 (MQ=255)
                      tACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                             <  1:3830955/100‑1 (MQ=255)
                       aCTGGCAGTGATAAGTCGGATGCGCTCGGCGTGTCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTccc                                                                            <  1:39454/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)
                                    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)
                                              gcTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCg                                                     <  2:2714731/100‑1 (MQ=255)
                                                   gCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGc                                                <  2:641884/100‑1 (MQ=255)
                                                    cGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                               <  1:3703181/100‑1 (MQ=255)
                                                    cGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACACCACCGCAAGcc                                               <  2:1770354/100‑1 (MQ=255)
                                                      tGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACATCCCACCTTTACCGCTTTACGCCACCGCAAGCCaa                                             >  1:1601432/1‑100 (MQ=255)
                                                       gCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGATTTACGCCACCGCAAGCCaaa                                            >  1:722798/1‑100 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                           <  1:2158649/100‑1 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                           <  1:2791152/100‑1 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                           >  2:1577621/1‑100 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                           >  1:1110623/1‑100 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACAGCTTTACGCCACCGCAAGCCAAAt                                           >  2:3200979/1‑100 (MQ=255)
                                                         ccTATTATAATCAACACGCTATGTAGTTTGTTCTGTCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                          <  2:2347760/100‑1 (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)
                                                          cTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATAc                                         >  1:2529868/1‑100 (MQ=255)
                                                          cTATAATAATCAACACGCTATGTAGATTGATCTAGCCCCGACATCTCGCGGCTTATTAACCTCCCACCTTCACCGCTTTAAGCCCCCGCAAACCAAATAc                                         >  2:2344870/1‑100 (MQ=255)
                                                             ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                      <  1:3779494/100‑1 (MQ=255)
                                                             ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                      >  1:1830596/1‑100 (MQ=255)
                                                             ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                      <  1:149578/100‑1 (MQ=255)
                                                             ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                      <  1:3807281/100‑1 (MQ=255)
                                                              tataATGAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                      <  2:393293/99‑1 (MQ=255)
                                                              tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGTCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTg                                     <  1:161036/100‑1 (MQ=255)
                                                                taATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGat                                   <  1:2793061/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)
                                                                             aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCgg                        >  1:1389764/1‑98 (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)
                                                                                        ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCACATACATTGATATACAGCCCGGTCATAATGAgc             >  2:3090491/5‑98 (MQ=255)
                                                                                        ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAAACATAGATATACAGCCCGGTCATAATGAgcac           >  1:2141615/5‑100 (MQ=255)
                                                                                          aagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccg         >  1:3888370/3‑100 (MQ=255)
                                                                                           ggCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc        <  1:1450587/100‑1 (MQ=255)
                                                                                           agCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc        >  2:2458034/2‑100 (MQ=255)
                                                                                           agCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc        >  1:1169366/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)
                                                                                             ccccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTAATAATGAgcaccgcac      >  1:851099/1‑100 (MQ=255)
                                                                                              cccGAGATGTCGGGGCCTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcacc     >  2:2018315/1‑100 (MQ=255)
                                                                                              cccGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgcacc     <  2:1606277/100‑1 (MQ=255)
                                                                                                cGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTa   <  1:2249139/100‑1 (MQ=255)
                                                                                                 gACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTaa  >  1:2240296/1‑100 (MQ=255)
                                                                                                 gACATCTCGAGGCTTATTCCCTTCCCACCTTTACCGCTTTATGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTa   <  2:2672149/99‑1 (MQ=255)
                                                                                                   |                                                                                                 
CGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAA  >  NC_000913/1620125‑1620321

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: