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.C92.8% 374.8 / 29.0 152intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/11);  new base C (8/133);  total (8/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01
*NC_0009133,957,0912.C92.8% 377.9 / 29.0 153intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/11);  new base C (8/134);  total (8/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 167 total aligned reads displayed.

AGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑CTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957002‑3957173
                                                                                              ||                                                                                      
aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑ccc                                                                                     <  1:5713765/93‑3 (MQ=255)
aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑ccc                                                                                     <  2:602621/93‑3 (MQ=255)
aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑ccc                                                                                     <  1:2461541/93‑3 (MQ=255)
 gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑ccct                                                                                    <  2:3752005/93‑1 (MQ=255)
  cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑Ccctc                                                                                   <  2:2178166/93‑1 (MQ=255)
  cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑Ccctc                                                                                   <  1:1370501/93‑1 (MQ=255)
  cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑Ccctc                                                                                   <  2:5042240/93‑1 (MQ=255)
    ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑Ccctctc                                                                                 <  2:3782884/93‑3 (MQ=255)
      taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑Ccctctccc                                                                               <  1:2949086/93‑5 (MQ=255)
        ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑Ccac                                                                             <  1:2920408/93‑1 (MQ=38)
          taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACag                                                                           <  1:4521467/93‑1 (MQ=37)
           aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagg                                                                          <  1:3506559/93‑1 (MQ=37)
           aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagg                                                                          <  1:3266285/93‑1 (MQ=37)
            aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACaggg                                                                         <  1:3934721/93‑1 (MQ=37)
            aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACaggg                                                                         <  2:829868/93‑1 (MQ=37)
            aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACaggg                                                                         <  1:1451057/93‑1 (MQ=37)
               aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagggaga                                                                      <  1:631773/93‑1 (MQ=37)
               aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagggaga                                                                      <  1:3706097/93‑1 (MQ=37)
               aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagggaga                                                                      <  2:865728/93‑1 (MQ=37)
                 cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagggagagg                                                                    <  1:2720606/93‑1 (MQ=25)
                  tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagggagaggg                                                                   <  2:6109133/93‑1 (MQ=25)
                   gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACagggagaggga                                                                  <  2:673026/93‑1 (MQ=25)
                       ggCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGaaa                                                             <  1:1639594/93‑1 (MQ=25)
                        gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCA‑‑CCCCC‑‑ACCCCCTC‑‑CCCCACGGAGAGGGGGaaaa                                                            >  1:6606902/1‑93 (MQ=21)
                          tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAaca                                                          <  2:5506543/93‑1 (MQ=21)
                           atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAacac                                                         <  1:2913276/93‑1 (MQ=21)
                            tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCCC‑‑C‑‑CTCTC‑‑CCCCCGGGAGAGGGGGAAAaaacc                                                        >  2:2005335/1‑92 (MQ=18)
                              tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCa                                                      <  2:2417895/93‑1 (MQ=21)
                              tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCa                                                      <  2:844228/93‑1 (MQ=21)
                              tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCa                                                      <  2:4523411/93‑1 (MQ=21)
                               gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCaa                                                     <  2:4440229/93‑1 (MQ=21)
                               gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:5125767/93‑1 (MQ=21)
                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCCCC‑‑C‑‑CCCTC‑‑CCCCAGGGGGGGGGGGAAAAAACACAAg                                                   >  2:969982/1‑93 (MQ=16)
                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:2996035/93‑1 (MQ=21)
                                   gATATGAACCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAggcc                                                <  1:2299866/93‑4 (MQ=18)
                                   gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCT‑‑CCCCG‑‑A‑‑CGCCCCGCCACCGGGGGGGGGAGAAAAAAGTCaaaggt                                                 >  1:2194852/1‑89 (MQ=19)
                                    atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  2:4014547/93‑5 (MQ=21)
                                    atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  2:3131501/93‑5 (MQ=21)
                                    atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  1:937749/93‑5 (MQ=21)
                                      atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCC‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:3809512/93‑1 (MQ=18)
                                      atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:5321362/93‑1 (MQ=18)
                                      atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGcccttc                                             <  2:3919462/93‑5 (MQ=21)
                                        gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:1458888/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCGGCC‑‑CTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:882053/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:4738124/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:746265/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:4713500/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:3111185/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:1582890/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:5252362/93‑1 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:4857127/93‑1 (MQ=18)
                                          aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  1:5491200/93‑1 (MQ=18)
                                           aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcccg                                         <  1:4200575/93‑3 (MQ=18)
                                             gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:4794914/93‑1 (MQ=18)
                                             gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:1357909/93‑1 (MQ=18)
                                             gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  1:1435846/93‑1 (MQ=18)
                                             gTATTGCTGTTGCGGGTTAAGTGCGCGCTAATGCCCTCA‑‑CCCCGCCC‑‑CCCTC‑‑CGACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  1:5770991/93‑1 (MQ=18)
                                              tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCC‑‑CCCCG‑‑ACCCTCTC‑‑CCCCCAGGAGAGGGAGGAAAAAATACAAGACTTCTCCTGgag                                      >  1:1058497/1‑93 (MQ=12)
                                              tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgag                                      <  1:4290702/93‑1 (MQ=18)
                                              tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgag                                      >  2:3392736/1‑93 (MQ=18)
                                               aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                     <  2:358088/93‑1 (MQ=18)
                                                ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:4538550/93‑1 (MQ=18)
                                                 tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  2:6287635/93‑1 (MQ=18)
                                                  gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                  <  2:4312446/93‑1 (MQ=18)
                                                    tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  2:3651621/93‑1 (MQ=18)
                                                    tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  2:4064021/93‑1 (MQ=18)
                                                     gTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  2:3400166/93‑1 (MQ=18)
                                                     gTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  2:1109236/93‑1 (MQ=18)
                                                     gTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:5193645/93‑1 (MQ=18)
                                                      ttGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  2:2828836/93‑1 (MQ=18)
                                                          gggTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                          <  2:2828260/93‑1 (MQ=18)
                                                          gggTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTa                          <  2:6342734/93‑1 (MQ=18)
                                                            gTTAAGTGCGCGCTGATGCCCTCG‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGGGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                        >  1:3128733/1‑93 (MQ=18)
                                                            gTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                        <  2:1409147/93‑1 (MQ=18)
                                                             ttAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  2:4852292/93‑1 (MQ=18)
                                                              tAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:6448866/93‑1 (MQ=18)
                                                               aaGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:639765/93‑1 (MQ=18)
                                                               aaGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:3669094/93‑1 (MQ=18)
                                                                aGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  2:5415298/93‑1 (MQ=18)
                                                                 gTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  1:2005227/93‑1 (MQ=18)
                                                                   gcgcgcTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  2:5802428/93‑1 (MQ=21)
                                                                    cgcgcTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                <  2:3691476/93‑1 (MQ=21)
                                                                     gcgcTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt               <  2:5975865/93‑1 (MQ=21)
                                                                      cgcTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTa              <  1:4677394/93‑1 (MQ=21)
                                                                       gcTGATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc             <  2:3586276/93‑1 (MQ=21)
                                                                         tGATGCCCTCA‑‑CCCCG‑‑ACCCTCTT‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg           >  1:127430/1‑93 (MQ=25)
                                                                          gATGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  1:1040366/93‑1 (MQ=25)
                                                                            tGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAACCACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc        <  1:1197763/93‑1 (MQ=25)
                                                                            tGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc        <  2:842131/93‑1 (MQ=25)
                                                                            tGCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAc        <  2:6625785/93‑1 (MQ=25)
                                                                             gCCCTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc       <  1:1713963/93‑1 (MQ=25)
                                                                             acccTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc       <  1:5117258/92‑1 (MQ=25)
                                                                              cccTCC‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  1:2575744/93‑1 (MQ=25)
                                                                              cccTCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  2:5656545/93‑1 (MQ=25)
                                                                                 tCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTATTTGCGGTACAGCACCTTaa   >  2:3220518/1‑93 (MQ=25)
                                                                                 tCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  2:1080743/93‑1 (MQ=37)
                                                                                 tCA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  2:2285145/1‑93 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  2:6578053/93‑1 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  >  1:1822179/1‑93 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCCTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:144886/93‑1 (MQ=37)
                                                                                              ||                                                                                      
AGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑CTCTC‑‑CCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957002‑3957173

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