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.5% 288.2 / 9.1 112intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/5);  new base C (6/101);  total (6/106)
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.49e-01
*NC_0009133,957,0912.C95.6% 290.0 / 9.0 113intergenic (+25/+59)ilvC/ppiCketol‑acid reductoisomerase/peptidyl‑prolyl cis‑trans isomerase C
Reads supporting (aligned to +/- strand):  ref base . (0/5);  new base C (6/102);  total (6/107)
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.52e-01

Only 100 of 123 total aligned reads displayed.

AGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957002‑3957173
                                                                                              ||                                                                                       
 gCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCC‑‑CCCCC‑‑A‑‑cccc                                                                                     >  2:1964581/1‑90 (MQ=255)
     gtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑C‑‑cctctcc                                                                               <  2:2528387/93‑4 (MQ=255)
         gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCcaca                                                                             <  1:1028594/93‑1 (MQ=38)
         gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCcaca                                                                             <  2:5555980/93‑1 (MQ=38)
          taAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACag                                                                            <  1:5296332/93‑1 (MQ=37)
           aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagg                                                                           <  2:1667688/93‑1 (MQ=37)
           aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagg                                                                           <  1:4814810/93‑1 (MQ=37)
           aaGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagg                                                                           <  2:559559/93‑1 (MQ=37)
            aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACaggg                                                                          <  2:1651523/93‑1 (MQ=37)
             gAACCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCC‑‑CCCCGGCC‑‑C‑‑TCTCCCGCaggga                                                                         <  2:4288244/93‑1 (MQ=25)
             gAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACaggga                                                                         <  2:3196311/93‑1 (MQ=37)
               aaCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagggaga                                                                       <  2:711556/93‑1 (MQ=37)
                aCTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagggagag                                                                      <  1:651080/93‑1 (MQ=37)
                aCTACGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagggagag                                                                      <  1:2707590/93‑1 (MQ=25)
                  tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagggagaggg                                                                    <  1:6298810/93‑1 (MQ=25)
                  tGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagggagaggg                                                                    <  2:4103744/93‑1 (MQ=25)
                   gcgcGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACagggagaggga                                                                   <  1:2328671/93‑1 (MQ=25)
                      cgGCTATATGACAGATATGAAGCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGaa                                                                <  1:1634713/93‑1 (MQ=21)
                      cgGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGaa                                                                <  2:2792994/93‑1 (MQ=25)
                        gCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGaaaa                                                              <  1:6002934/93‑1 (MQ=21)
                         cTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAac                                                             <  1:4512987/93‑1 (MQ=21)
                          tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCC‑‑CCCCC‑‑ACCC‑‑CCTCCCCCCAGGAGAGGGGGaaaaaa                                                            >  1:5894655/1‑91 (MQ=21)
                          tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAaca                                                            <  1:2881568/93‑1 (MQ=21)
                           atatGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAacac                                                           <  2:692652/93‑1 (MQ=21)
                             atGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTc                                                        <  1:2920941/93‑1 (MQ=21)
                               gACAGATATGACACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  2:4598999/93‑1 (MQ=21)
                               gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGGCCTCACCCCCCC‑‑T‑‑C‑‑TCTCCCCCATGGAGGGGGAGAGAACACTaaa                                                     >  2:5140659/1‑93 (MQ=255)
                               gACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                     <  2:4549254/93‑1 (MQ=21)
                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:700721/93‑1 (MQ=21)
                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:1503394/93‑1 (MQ=21)
                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:2502993/93‑1 (MQ=21)
                                aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                   <  2:141965/93‑1 (MQ=21)
                                  aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                 <  1:3237166/93‑3 (MQ=21)
                                  aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAggc                                                 <  2:5546462/93‑3 (MQ=21)
                                   gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCT‑‑ACCC‑‑TCTTCCACATGGGGGGGGAGAAAACACTCaacgcc                                                 >  1:3186964/1‑89 (MQ=21)
                                   gATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                <  2:4692227/93‑4 (MQ=21)
                                   aatatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAggcc                                                <  1:3879865/92‑4 (MQ=18)
                                    atatGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAggcct                                               <  2:3063748/93‑5 (MQ=21)
                                      atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:5065353/93‑1 (MQ=18)
                                      atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTtc                                              <  2:5698581/93‑1 (MQ=18)
                                        gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCC‑‑ACCC‑‑TCTCCCCCAGGGAGAGGGAGAAAACACTCAAAGCCTtctc                                            >  2:2023207/1‑93 (MQ=18)
                                         aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                           <  1:2458931/93‑1 (MQ=18)
                                          aaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCt                                          <  2:4530529/93‑1 (MQ=18)
                                           aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  1:5029976/93‑1 (MQ=18)
                                           aCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTg                                         <  1:4757539/93‑1 (MQ=18)
                                             gTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGga                                       <  1:569997/93‑1 (MQ=21)
                                               aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                     <  2:4218555/93‑1 (MQ=18)
                                               aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                     <  1:1870230/93‑1 (MQ=18)
                                               aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGgaga                                     <  1:4922320/93‑1 (MQ=18)
                                               aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCCCA‑‑CCCCGCCC‑‑C‑‑CCTCCCCCCGGGGGGGGGGGGAAAAACTACAAGCCTTCTCCTggggg                                     >  2:4141417/1‑90 (MQ=12)
                                                ttGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  1:4267159/93‑1 (MQ=18)
                                                ttGCTGTAGGGGGTTAAGTGCGCGCTGATCCCCTCC‑‑CCCCG‑‑ACCC‑‑CCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                    <  2:5302792/93‑1 (MQ=18)
                                                 tGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                   <  1:650441/93‑1 (MQ=18)
                                                  gCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCCCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAgg                                  <  2:3574081/93‑1 (MQ=18)
                                                   cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                 <  2:3527743/93‑1 (MQ=18)
                                                   cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                 <  2:3143398/93‑1 (MQ=18)
                                                     tttGCGGGTTAAGTGCTCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGAGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                               <  2:6271825/92‑1 (MQ=18)
                                                      ttGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  2:6101429/93‑1 (MQ=18)
                                                      ttGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  1:2729665/93‑1 (MQ=18)
                                                      gtGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCtt                              <  1:605021/92‑1 (MQ=18)
                                                        gCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGGGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            >  1:3193077/1‑93 (MQ=18)
                                                        gCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  1:1570314/93‑1 (MQ=18)
                                                        gCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGc                            <  2:3180618/93‑1 (MQ=18)
                                                             ttAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  2:293587/93‑1 (MQ=18)
                                                             ttAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  1:4844363/93‑1 (MQ=18)
                                                             ttAAGTGCGCGCTGATGACCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTa                       <  1:2441000/93‑1 (MQ=18)
                                                              tAAGTGCGCGCTGATGCCCTCG‑‑CCCCGCCC‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:4723319/93‑1 (MQ=18)
                                                              tAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAg                      <  1:4031086/93‑1 (MQ=18)
                                                               aaGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGt                     <  1:970136/93‑1 (MQ=18)
                                                                aGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  1:3317814/93‑1 (MQ=18)
                                                                aGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGtt                    <  2:4489672/93‑1 (MQ=18)
                                                                  tGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  1:2513177/93‑1 (MQ=21)
                                                                  tGCGCGCTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  1:6160169/93‑1 (MQ=21)
                                                                  ggctcgcTGATGCCTTCG‑‑CCCGG‑‑A‑‑CCGTCTCCCGCAGGGAGAGGGAGAAAACTCTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGc                  <  2:1459393/89‑1 (MQ=19)
                                                                   gcgcgcTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  2:5962141/93‑1 (MQ=21)
                                                                   gcgcgcTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCg                 <  2:2213932/93‑1 (MQ=21)
                                                                    cgcgcTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCgg                <  1:6154697/93‑1 (MQ=21)
                                                                     gcgcTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt               <  2:4504717/93‑1 (MQ=21)
                                                                     gcgcTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGt               <  1:489710/93‑1 (MQ=21)
                                                                       gcTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTAc             <  2:23900/93‑1 (MQ=21)
                                                                        cggATGCCCTCC‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACGCTCAAGGCCTCCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa            <  2:3487822/91‑1 (MQ=21)
                                                                        cTGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACa            <  2:316080/93‑1 (MQ=21)
                                                                         tGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg           <  2:1126514/93‑1 (MQ=25)
                                                                         tGATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAg           <  1:462216/93‑1 (MQ=25)
                                                                          gATGCCCTCC‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  1:93596/93‑1 (MQ=25)
                                                                          gATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  1:931863/93‑1 (MQ=25)
                                                                          gATGCCCTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGc          <  2:508328/93‑1 (MQ=25)
                                                                              cccTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  1:6237745/93‑1 (MQ=25)
                                                                              cccTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  2:3325020/93‑1 (MQ=25)
                                                                              cccGCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCt      <  1:4880619/93‑1 (MQ=25)
                                                                               ccTCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCtt     <  2:6017890/93‑1 (MQ=25)
                                                                                 tCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   <  2:903249/93‑1 (MQ=37)
                                                                                 tCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  2:5958298/1‑93 (MQ=37)
                                                                                 tCA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTaa   >  2:4866222/1‑93 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:3932169/93‑1 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:2849179/93‑1 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  2:4788886/93‑1 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:6144739/93‑1 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  <  1:1871355/93‑1 (MQ=37)
                                                                                  cA‑‑CCCCG‑‑ACCC‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAat  >  2:4190781/1‑93 (MQ=37)
                                                                                              ||                                                                                       
AGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCA‑‑CCCCG‑‑A‑‑C‑‑TCTCCCACAGGGAGAGGGAGAAAACACTCAAGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAAT  >  NC_000913/3957002‑3957173

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