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.C93.9% 324.3 / 18.0 132intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/8);  new base C (2/122);  total (2/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
*NC_0009133,957,0912.C94.0% 326.6 / 18.0 133intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/8);  new base C (2/123);  total (2/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 158 total aligned reads displayed.

TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3956999‑3957173
                                                                                                  ||                                                                                      
ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑a                                                                                          <  1:2515700/93‑1 (MQ=255)
 tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑c                                                                                       <  1:1783410/93‑1 (MQ=255)
 tGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑c                                                                                       <  2:1986456/93‑1 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑ccc                                                                                     <  1:4590245/93‑3 (MQ=255)
   aGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑ccc                                                                                     <  1:3727938/93‑3 (MQ=255)
     cAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑C‑‑cctc                                                                                 <  2:502238/93‑1 (MQ=255)
         taggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGTTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑C‑‑cctctccc                                                                             <  2:4327998/93‑5 (MQ=38)
              aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagg                                                                          <  1:66379/93‑1 (MQ=37)
               aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACaggg                                                                         <  2:1011315/93‑1 (MQ=37)
               aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACaggg                                                                         <  1:3123040/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACaggga                                                                        <  2:2283446/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACaggga                                                                        <  2:4022670/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACaggga                                                                        <  2:1193931/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACaggga                                                                        <  2:2025620/93‑1 (MQ=37)
                gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACaggga                                                                        <  1:1973708/93‑1 (MQ=37)
                  aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggaga                                                                      <  1:2454025/93‑1 (MQ=37)
                  aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggaga                                                                      <  1:2716087/93‑1 (MQ=37)
                  aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggaga                                                                      <  2:1616798/93‑1 (MQ=37)
                   aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggagag                                                                     <  2:837989/93‑1 (MQ=37)
                   aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggagag                                                                     <  1:4013836/93‑1 (MQ=37)
                   aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggagag                                                                     <  2:656650/93‑1 (MQ=37)
                    cTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggagagg                                                                    <  2:3526764/93‑1 (MQ=25)
                     tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggagaggg                                                                   <  2:1063633/93‑1 (MQ=25)
                       cgcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACagggagagggag                                                                 <  1:3791417/93‑1 (MQ=25)
                        gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAgggagagggaga                                                                <  1:4264171/93‑1 (MQ=25)
                        gcgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑CCTCCCACAcggagcggggga                                                                >  1:4623196/1‑93 (MQ=25)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGaaaa                                                             <  2:476662/93‑1 (MQ=21)
                           gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGaaaa                                                             <  2:4507182/93‑1 (MQ=21)
                            cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAac                                                            <  2:2715881/93‑1 (MQ=21)
                            cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAac                                                            <  1:2579504/93‑1 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑GCCCCCCCCCC‑‑ACCC‑‑CCTCCCCCCCGGAGAGGGAGgaaaaa                                                           >  2:4620486/1‑91 (MQ=255)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAaca                                                           <  1:3395608/93‑1 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAaca                                                           <  1:937925/93‑1 (MQ=21)
                             tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAaca                                                           <  2:837859/93‑1 (MQ=21)
                              atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAacac                                                          <  2:570964/93‑1 (MQ=21)
                                atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTc                                                        <  1:2725439/93‑1 (MQ=21)
                                 tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCCCACTCCCCC‑‑C‑‑C‑‑CCTCCCACCCCGCGAAGGAGAAAAcacaca                                                      >  2:4168807/1‑93 (MQ=255)
                                 tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCCCACCCCC‑‑ACCC‑‑CCTCCCCCCGGGGGGGGGGGGAAACACTCa                                                      >  1:1454056/1‑93 (MQ=255)
                                      gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACCGGGAGAGGGAGAAAACACTCAAggcc                                                >  2:705210/1‑90 (MQ=21)
                                       atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  1:1365027/93‑5 (MQ=21)
                                       atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  1:2156561/93‑5 (MQ=21)
                                         atGAAGCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:4422825/93‑1 (MQ=18)
                                         atGAAGCGTATTGCTGTTGCGGGTTAAGTGCGCCCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:4218278/93‑1 (MQ=18)
                                         atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  1:2604471/93‑1 (MQ=18)
                                           gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:3497452/93‑1 (MQ=18)
                                           gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtctc                                            <  2:1551671/93‑1 (MQ=18)
                                            aaaCGTATTGCTTTTGCGGTTTAAGTGCGCGCTGAGG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:497037/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:1021635/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  2:2060299/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:191149/93‑1 (MQ=18)
                                            aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:749614/93‑1 (MQ=18)
                                               cGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTgg                                        <  2:1454695/93‑1 (MQ=18)
                                                gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:2636398/93‑1 (MQ=18)
                                                gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  2:1276195/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:162991/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:4495315/93‑1 (MQ=18)
                                                   ttGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  2:3295827/93‑1 (MQ=18)
                                                    tGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  2:1612940/93‑1 (MQ=18)
                                                      cTGTTGCAGGTTAAGCGTGCGCTGATA‑‑CCACTACTCCG‑‑A‑‑CCGTCTCCTACAGGGAGAGGGAGAAAAAACTCGAGGCATGCTCCTGGAGAAGGc                                 <  2:1596728/93‑1 (MQ=12)
                                                       tGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGcc                                <  1:4195381/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  2:2772129/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  1:4620925/93‑1 (MQ=18)
                                                        gTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCGTCTCCTGGAGAAGGCCt                               <  2:760563/93‑1 (MQ=18)
                                                         ttGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  1:411107/93‑1 (MQ=18)
                                                           tcGGGTTACGCGCGCTCTCATG‑‑CCCCCACCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACCCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:1112027/92‑1 (MQ=12)
                                                           gCGGGTTAAGTGCGCGTTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:28612/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGTCTTCTCCTGGAGAAGGCCTTGc                            <  1:3569633/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:1026310/93‑1 (MQ=18)
                                                           gCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:3414100/93‑1 (MQ=18)
                                                            cGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCt                           <  1:3827185/93‑1 (MQ=18)
                                                               gTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTAtt                        <  1:2287765/93‑1 (MQ=18)
                                                                ttAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  2:3801561/93‑1 (MQ=18)
                                                                ttAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  1:1872305/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:2939972/93‑1 (MQ=18)
                                                                  aaGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTACCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:1702174/93‑1 (MQ=18)
                                                                  aaGGGCGCGCTGATG‑‑CCCTCACCCCGGCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:936676/93‑1 (MQ=17)
                                                                  aaGGGCGCGCTGATG‑‑CCCTCACCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  2:2847661/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGGTG‑‑TCCTCACCGTCG‑‑ACC‑‑TCTCCGCCAGGGAGAGGGAGAAAACACTCAAGGCCTTCCCCGGGAGAGGCCCTTGCTATTAgct                    <  2:734400/93‑3 (MQ=12)
                                                                   aGTGCGCGCTGATG‑‑CCCTCCCCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACCCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  2:1228243/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:376387/93‑1 (MQ=18)
                                                                   aGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:504056/93‑1 (MQ=18)
                                                                   aGTGCGCGCGGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  2:2958789/93‑1 (MQ=18)
                                                                     tGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCGAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  2:380806/93‑1 (MQ=18)
                                                                     tGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  2:3366499/93‑1 (MQ=21)
                                                                     tGCGCGCTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  1:3772101/93‑1 (MQ=21)
                                                                     gccgcgcAGATG‑‑CCCTCCCCCCG‑‑ACCC‑‑TCTCCGACAGGGATAGAGAGAAAACACACAAGGCCTTCTCCTGGAGAAGTCCTTGCTATTAGTTGc                  <  2:2475745/91‑1 (MQ=255)
                                                                       cgcgcTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                <  1:3637716/93‑1 (MQ=21)
                                                                          gcTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc             <  1:2040369/93‑1 (MQ=21)
                                                                          gcTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc             <  2:1408232/93‑1 (MQ=21)
                                                                           cTGATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa            <  1:940577/93‑1 (MQ=21)
                                                                             gATG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAgt          <  1:3639552/93‑2 (MQ=25)
                                                                              aTG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAgta         <  2:4281242/93‑3 (MQ=25)
                                                                               tG‑‑CCCTCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTAc        <  1:3936794/93‑1 (MQ=25)
                                                                                   cccTCACCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCt      <  2:1689805/93‑1 (MQ=21)
                                                                                      tCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTaa   <  2:3089731/93‑1 (MQ=25)
                                                                                      tCACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTaa   <  2:2362953/93‑1 (MQ=25)
                                                                                       cACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAat  <  1:2426297/93‑1 (MQ=25)
                                                                                       cACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAat  <  1:814439/93‑1 (MQ=25)
                                                                                       cACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAat  <  1:1476726/93‑1 (MQ=25)
                                                                                       cACCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGTACCTTAat  <  1:1326695/93‑1 (MQ=25)
                                                                                                  ||                                                                                      
TTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATG‑‑CCCTCACCCCG‑‑A‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3956999‑3957173

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