Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,957,092 +CC intergenic (+26/+58) ilvC → / ← ppiC ketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,957,0911.C95.4% 401.8 / 16.1 153intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/7);  new base C (5/141);  total (5/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01
*NC_0009133,957,0912.C95.4% 402.6 / 16.1 153intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/7);  new base C (5/141);  total (5/148)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

Only 100 of 167 total aligned reads displayed.

TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3956999‑3957173
                                                                                               ||                                                                                       
 tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑c                                                                                        <  2:3254375/93‑1 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                                      <  2:1551755/93‑3 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                                      <  2:3112562/93‑3 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc                                                                                      <  2:5164218/93‑3 (MQ=255)
    gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑ccc‑‑‑t                                                                                  <  1:703850/93‑1 (MQ=255)
     cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccc‑‑‑tc                                                                                 <  1:1021298/93‑1 (MQ=255)
     cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccc‑‑‑tc                                                                                 <  2:809937/93‑1 (MQ=255)
     cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑Ccc‑‑‑tc                                                                                 <  1:1555842/93‑1 (MQ=255)
            gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCcaca                                                                            <  2:408603/93‑1 (MQ=38)
             taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACag                                                                           <  1:2864458/93‑1 (MQ=37)
              aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagg                                                                          <  1:2420264/93‑1 (MQ=37)
               aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggg                                                                         <  1:3184175/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCGCCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  1:3927307/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:6031079/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:5567110/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:4771747/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACaggga                                                                        <  2:2121284/93‑1 (MQ=37)
                   aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagag                                                                     <  1:80414/93‑1 (MQ=37)
                    cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagagg                                                                    <  2:1259763/93‑1 (MQ=25)
                    cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagagg                                                                    <  1:4563853/93‑1 (MQ=25)
                     tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagaggg                                                                   <  1:152831/93‑1 (MQ=25)
                      gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagaggga                                                                  <  1:177466/93‑1 (MQ=25)
                       cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACagggagagggag                                                                 <  2:3759527/93‑1 (MQ=25)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGaaaa                                                             <  1:5135331/93‑1 (MQ=21)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCC‑‑ACCCCC‑‑‑TCCCACAGGGAGCGGGAGaaaa                                                             >  2:3279794/1‑93 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAaca                                                           <  1:5849000/93‑1 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCACCCCC‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAaca                                                           >  2:5469641/1‑93 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCAACCCGACA‑‑CTC‑‑‑TCCCCCCAGGGGAGGGGGAAAaca                                                           >  1:3343323/1‑93 (MQ=25)
                              atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAacac                                                          <  2:2783006/93‑1 (MQ=21)
                                atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTc                                                        <  2:1174069/93‑1 (MQ=21)
                                atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTc                                                        <  2:2793655/93‑1 (MQ=21)
                                  gACAGATATGAAGCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGCCC‑‑CGC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  2:3294301/93‑1 (MQ=18)
                                  gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  2:1806617/93‑1 (MQ=21)
                                  gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:2508504/93‑1 (MQ=21)
                                  gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:1904400/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  1:1973059/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:397507/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  1:3033432/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:675265/93‑1 (MQ=21)
                                   aCAGATATGAAACGTATTGCTGTTGCAGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  1:5979940/93‑1 (MQ=21)
                                    cAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAgg                                                  <  1:5549775/93‑2 (MQ=21)
                                     aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                 <  2:1443408/93‑3 (MQ=21)
                                      gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                <  1:9338/93‑4 (MQ=21)
                                       atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  2:5607446/93‑5 (MQ=21)
                                       atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAAGGAGAAAACACTCAAggcct                                               <  2:6649291/93‑5 (MQ=21)
                                         atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:2687964/93‑1 (MQ=18)
                                         atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:2281234/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:392906/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:4103406/93‑1 (MQ=18)
                                             accGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCCCCCCGCCC‑‑CTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  1:4575859/91‑1 (MQ=18)
                                             aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  2:4929641/93‑1 (MQ=18)
                                              aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:3706130/93‑1 (MQ=18)
                                              aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:3890355/93‑1 (MQ=21)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  1:2062873/93‑1 (MQ=18)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  1:5716829/93‑1 (MQ=18)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  2:624573/93‑1 (MQ=18)
                                                gTATTGCTGTTGCGGGTGAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:399765/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCTCACCGCGAGAGGGAGAAAAAACGCAAAGCCTTCTCCTGGAGaa                                    >  2:4808822/1‑93 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:3736839/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:5425655/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:2692190/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:4010176/93‑1 (MQ=18)
                                                    tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  2:888500/93‑1 (MQ=18)
                                                      cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                 <  2:5464697/93‑1 (MQ=18)
                                                       tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  1:1270206/93‑1 (MQ=18)
                                                       tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  2:4740879/93‑1 (MQ=18)
                                                        tttGCGGGTTAAGTTCGCGCTGATCCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGTGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:552755/92‑1 (MQ=18)
                                                        gTTGCGGTTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:6429211/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:5530893/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  2:4337576/93‑1 (MQ=18)
                                                         ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  1:4950168/93‑1 (MQ=18)
                                                         ttGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  2:5130233/93‑1 (MQ=18)
                                                          tGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTg                             <  2:1958098/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:657926/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:3748690/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:5063741/93‑1 (MQ=18)
                                                            cGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                           <  2:3034061/93‑1 (MQ=18)
                                                             gggTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                          <  1:4431487/93‑1 (MQ=18)
                                                              ggTTAAGTGCGCGCTGCTGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         <  1:398457/93‑1 (MQ=18)
                                                              ggTTAAGTGCGCGCTGATGCCCCCCCCCCC‑‑C‑‑CCCCCTTCCCCCGGGGGGGGGGGGAAAAAACCCAAGCCCTCCCCTGGGGAAGGCCCTGCttt                         >  2:4506482/1‑91 (MQ=255)
                                                               gTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                        <  2:2273293/93‑1 (MQ=18)
                                                                 tAAGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  2:5187312/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:3226459/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:1487841/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  2:2729588/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:5071597/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:5640282/93‑1 (MQ=18)
                                                                    gTGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  2:1970743/93‑1 (MQ=18)
                                                                     tGCGCGCTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  2:4398858/93‑1 (MQ=21)
                                                                      gcgcgcTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  1:2529276/93‑1 (MQ=21)
                                                                       cgcgcTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                <  2:6278852/93‑1 (MQ=21)
                                                                       cgcgcTGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCGCCTGGAGAAGGCCTTGCTATTAGTTGCgg                <  1:6403099/93‑1 (MQ=18)
                                                                            tGATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg           <  1:2980254/93‑1 (MQ=25)
                                                                             gATGCCCTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  2:3058079/93‑1 (MQ=25)
                                                                                  ccTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGGGAGGGAGAAAACACTCAAGGCCTTTTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     >  2:3315625/1‑93 (MQ=25)
                                                                                  ccTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     <  1:5960311/93‑1 (MQ=25)
                                                                                   cTCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa    >  1:6513542/1‑93 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  1:1019747/93‑1 (MQ=37)
                                                                                    tCACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  1:2986405/93‑1 (MQ=37)
                                                                                     cACCCCG‑‑ACCCTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:79103/93‑1 (MQ=37)
                                                                                               ||                                                                                       
TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCG‑‑A‑‑CTC‑‑‑TCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3956999‑3957173

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