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.C91.2% 808.6 / 77.4 321intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (9/19);  new base C (33/259);  total (42/279)
Fisher's exact test for biased strand distribution p-value = 5.12e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,957,0912.C91.0% 808.5 / 80.9 321intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (10/19);  new base C (32/260);  total (42/279)
Fisher's exact test for biased strand distribution p-value = 1.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

TGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957000‑3957173
                                                                                                   ||                                                                                     
tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑c                                                                                      <  2:4644303/93‑1 (MQ=255)
tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑c                                                                                      <  1:562221/93‑1 (MQ=255)
   gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑ccct                                                                                   <  2:2273340/93‑1 (MQ=255)
   gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑ccct                                                                                   <  1:2973484/93‑1 (MQ=255)
    cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑Ccctc                                                                                  <  2:4864663/93‑1 (MQ=255)
    cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑Ccctc                                                                                  <  2:3886969/93‑1 (MQ=255)
    cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑Ccctc                                                                                  <  1:3202057/93‑1 (MQ=255)
    cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCAGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑Ccctc                                                                                  <  2:7217282/93‑1 (MQ=255)
       gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑GCCTCA‑‑CCCCC‑‑A‑‑Ccctctcc                                                                               >  2:3634167/1‑90 (MQ=255)
       gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCC‑‑CCCCG‑‑A‑‑Ccctctcc                                                                               <  1:2817074/93‑4 (MQ=255)
        taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑GCCTCA‑‑CCCCG‑‑A‑‑Ccctctccc                                                                              >  1:1125637/1‑89 (MQ=255)
          ggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCcac                                                                              <  2:6564739/93‑1 (MQ=38)
             aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACagg                                                                           <  2:1687157/93‑1 (MQ=37)
             aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACagg                                                                           <  1:4452697/93‑1 (MQ=37)
               gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACaggga                                                                         <  2:3097597/93‑1 (MQ=37)
                 aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACagggaga                                                                       <  2:819394/93‑1 (MQ=37)
                   cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACagggagagg                                                                     <  1:1391038/93‑1 (MQ=25)
                   cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACagggagagg                                                                     <  1:4609895/93‑1 (MQ=25)
                    tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACagggagaggg                                                                    <  2:3066985/93‑1 (MQ=25)
                     gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATCG‑‑CCCCCA‑‑CCCCG‑‑ACCCCCCCCCACagggtgaggg                                                                    >  1:6688712/1‑93 (MQ=255)
                      cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACagggagagggag                                                                  <  1:7408261/93‑1 (MQ=25)
                       gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCCCA‑‑CCCCC‑‑ACCCTCTCCCCCCgggaggggggga                                                                >  1:2323734/1‑93 (MQ=25)
                        cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaa                                                               <  2:7046330/93‑1 (MQ=25)
                        cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaa                                                               <  2:3721309/93‑1 (MQ=25)
                          gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGaaaa                                                             <  2:7160109/93‑1 (MQ=21)
                            tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAaca                                                           <  1:5462008/93‑1 (MQ=21)
                            tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAaca                                                           <  2:5663519/93‑1 (MQ=21)
                            tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCGCA‑‑CCCCGGCC‑‑CTCTCCCACAGGGAGAGGGGGgaaaaa                                                           >  2:4394582/1‑91 (MQ=18)
                             atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCC‑‑CCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAGAAAacac                                                          <  1:4994329/93‑1 (MQ=21)
                             atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCCCC‑‑C‑‑CTCTCCCACCGGGAGGGGGGGGAAacac                                                          >  1:4956730/1‑93 (MQ=21)
                             atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAacac                                                          <  1:3811649/93‑1 (MQ=21)
                              tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACt                                                         <  2:6459125/93‑1 (MQ=21)
                              tatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCCCC‑‑CCCCC‑‑ACCCCCTCCCCCCGGGAGGGGGAGgaaaaaaa                                                         >  2:3456809/1‑89 (MQ=19)
                               atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTc                                                        <  1:3277699/93‑1 (MQ=21)
                                tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCa                                                      >  2:4986826/1‑93 (MQ=21)
                                 gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:4222819/93‑1 (MQ=21)
                                 gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  1:2296873/93‑1 (MQ=21)
                                  aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  1:999425/93‑1 (MQ=21)
                                    aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                 <  2:7447878/93‑3 (MQ=21)
                                     gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                <  2:4814619/93‑4 (MQ=21)
                                      atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCC‑‑CCCCGCCC‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  1:6137445/93‑5 (MQ=18)
                                      atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  2:5410152/93‑5 (MQ=21)
                                        atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:47468/93‑1 (MQ=18)
                                        atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:3573679/93‑1 (MQ=18)
                                         tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                             <  2:407865/93‑1 (MQ=21)
                                         tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                             <  2:4309308/93‑1 (MQ=18)
                                         tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtct                                             <  2:2343567/93‑1 (MQ=18)
                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:2106839/93‑1 (MQ=18)
                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:595887/93‑1 (MQ=18)
                                           aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:3420256/93‑1 (MQ=18)
                                            aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  2:6237566/93‑1 (MQ=18)
                                             aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:1786916/93‑1 (MQ=18)
                                             aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  1:6632012/93‑1 (MQ=18)
                                             aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  1:7073373/93‑1 (MQ=18)
                                             aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  1:2648754/93‑1 (MQ=18)
                                             aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  2:2742789/93‑1 (MQ=18)
                                                tATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  1:1208530/93‑1 (MQ=18)
                                                  ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  2:3256808/93‑1 (MQ=18)
                                                   tGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  2:6123468/93‑1 (MQ=18)
                                                    gCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                  <  2:5822223/93‑1 (MQ=18)
                                                    gCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                  <  1:788244/93‑1 (MQ=18)
                                                     cTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                 <  2:2890828/93‑1 (MQ=18)
                                                       gTTGCGGGTTAAGTGCGCGCTGATGCCCCCCCC‑‑CCCCC‑‑C‑‑CCCTCCCCCAGGGGGGGGGGGAAAAAACTCAAAGCCCTCTCCCTGGGAAgggcc                               >  1:3721389/1‑90 (MQ=9)
                                                          gCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:6452089/93‑1 (MQ=18)
                                                          gCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:525219/93‑1 (MQ=18)
                                                            gggTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCCTCCCCCACGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         >  2:473366/1‑93 (MQ=21)
                                                             ggTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAt                         <  2:5961346/93‑1 (MQ=18)
                                                               ttAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  2:2143146/93‑1 (MQ=18)
                                                                tAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:1856043/93‑1 (MQ=18)
                                                                tAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  2:5297196/93‑1 (MQ=18)
                                                                tAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:1985038/93‑1 (MQ=18)
                                                                tAAGTGCGCGCTGAAG‑‑CCCCCC‑‑CCCCGGCC‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:5504985/93‑1 (MQ=18)
                                                                 aaGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:19776/93‑1 (MQ=18)
                                                                 aaGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:7473852/93‑1 (MQ=18)
                                                                  aGTGCGCGCTGATG‑‑CCCCCA‑‑CCCCC‑‑ACCCCCTCCCCCCGGGGGGGGGAGAAAAAAACCAAAGGCTTCTCCTCGAGAAGGGCTTGCTATTAGtt                    >  1:1032302/1‑93 (MQ=9)
                                                                   nngcgcgcTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  2:1351512/91‑1 (MQ=18)
                                                                   gTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  2:2171043/93‑1 (MQ=18)
                                                                   gTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTg                   <  2:6876363/93‑1 (MQ=18)
                                                                    tGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  2:7322482/93‑1 (MQ=21)
                                                                     gcgcgcTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  1:7341242/93‑1 (MQ=21)
                                                                       gcgcTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt               <  1:1480774/93‑1 (MQ=21)
                                                                         gcTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc             <  2:2843395/93‑1 (MQ=21)
                                                                          cTGATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa            <  2:301339/93‑1 (MQ=21)
                                                                            gATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  2:1128213/93‑1 (MQ=25)
                                                                            gATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  2:4260743/93‑1 (MQ=25)
                                                                            gATG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  2:2560913/93‑1 (MQ=25)
                                                                             aTG‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCa         <  1:1168425/93‑1 (MQ=25)
                                                                                g‑‑CCCTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCAcc       <  2:6318590/93‑1 (MQ=25)
                                                                                   cccTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  1:1558610/93‑1 (MQ=25)
                                                                                   cccTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  2:6784023/93‑1 (MQ=25)
                                                                                   cccTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  1:1791919/93‑1 (MQ=25)
                                                                                   cccTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  2:1860029/93‑1 (MQ=25)
                                                                                     cTCA‑‑CCCCG‑‑ACCCTCTCCCCCAGGGGGAGGGGGAAAAAACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa    >  2:2480432/1‑93 (MQ=255)
                                                                                     cTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa    <  1:5082802/93‑1 (MQ=37)
                                                                                     cTCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTa    <  1:6424387/93‑1 (MQ=37)
                                                                                      tCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  1:4340204/93‑1 (MQ=37)
                                                                                      tCA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  1:2478433/93‑1 (MQ=37)
                                                                                       cA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:4383812/93‑1 (MQ=37)
                                                                                       cA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:6055371/93‑1 (MQ=37)
                                                                                       cA‑‑CCCCG‑‑ACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:3694200/93‑1 (MQ=37)
                                                                                                   ||                                                                                     
TGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCA‑‑CCCCG‑‑A‑‑CTCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957000‑3957173

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