Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,597,825 A→G intergenic (‑45/‑101) dapA ← / → b2479 dihydrodipicolinate synthase/GcvR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,597,8250AG99.4% 559.6 / ‑4.7 176intergenic (‑45/‑101)dapA/b2479dihydrodipicolinate synthase/GcvR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (90/85);  total (90/86)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 176 total aligned reads displayed.

CATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAG  >  NC_000913/2597738‑2597914
                                                                                       |                                                                                         
 atccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAt                                                                                     >  1:855619/1‑93 (MQ=255)
 atccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAt                                                                                     <  1:4304917/93‑1 (MQ=255)
 atccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAt                                                                                     >  2:372664/1‑93 (MQ=255)
  tccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAtt                                                                                    >  2:4121443/1‑93 (MQ=255)
  tccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAtt                                                                                    >  2:132022/1‑93 (MQ=255)
  tccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAtt                                                                                    >  1:1863223/1‑93 (MQ=255)
   ccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTa                                                                                   <  2:2607311/93‑1 (MQ=255)
    catcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaa                                                                                  <  1:79397/93‑1 (MQ=255)
     atcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaa                                                                                 >  2:2459961/1‑93 (MQ=255)
      tcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaaa                                                                                >  1:3931445/1‑93 (MQ=255)
        ggAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagc                                                                              >  2:3040978/1‑93 (MQ=255)
         gAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagca                                                                             >  1:4106163/1‑93 (MQ=255)
          aGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaa                                                                            >  1:2860728/1‑93 (MQ=255)
           gTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaag                                                                           <  1:1010597/93‑1 (MQ=255)
           gTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaag                                                                           >  2:1080809/1‑93 (MQ=255)
           gTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaag                                                                           >  2:3489937/1‑93 (MQ=255)
            tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaagc                                                                          <  1:1968658/93‑1 (MQ=255)
            tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaagc                                                                          >  2:1241875/1‑93 (MQ=255)
            tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaagc                                                                          >  2:1616079/1‑93 (MQ=255)
             aaCAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAagcaagca                                                                         <  2:3738172/93‑1 (MQ=255)
              aCAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcag                                                                        >  2:3192702/1‑93 (MQ=255)
               cAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcaga                                                                       <  1:4049815/93‑1 (MQ=255)
               cAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcaga                                                                       <  2:4426269/93‑1 (MQ=255)
               cAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGAAAGcaga                                                                       >  2:1933001/1‑93 (MQ=255)
                aaTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagac                                                                      <  1:996119/93‑1 (MQ=255)
                aaTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagac                                                                      >  1:2185121/1‑93 (MQ=255)
                aaTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagaa                                                                      >  1:2098973/1‑92 (MQ=255)
                 aTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagaca                                                                     <  1:1909731/93‑1 (MQ=255)
                 aTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagaca                                                                     >  1:1260778/1‑93 (MQ=255)
                    gcgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGaa                                                                  <  2:1214164/93‑1 (MQ=255)
                     cgACAATACTTCCCGTGAGCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                                 <  1:2369824/93‑1 (MQ=255)
                     cgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                                 >  1:232610/1‑93 (MQ=255)
                      gACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAcc                                                                >  2:2101903/1‑93 (MQ=255)
                        cAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGt                                                              <  1:1856447/93‑1 (MQ=255)
                         aaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGtt                                                             >  1:1658880/1‑93 (MQ=255)
                          aTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTc                                                            >  2:4258769/1‑93 (MQ=255)
                          aTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTc                                                            <  1:1060416/93‑1 (MQ=255)
                              ttCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAt                                                        <  1:3186825/93‑1 (MQ=255)
                              ttCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAt                                                        <  1:534425/93‑1 (MQ=255)
                               tCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAtt                                                       <  1:1106122/93‑1 (MQ=255)
                               tCCCGTGAACAGGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTGTGGTGTGGCATTAAAAGCAAGCAGACGGAACCGTGCTGAtt                                                       <  2:2997411/93‑1 (MQ=255)
                                 ccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgt                                                     >  2:4296247/1‑93 (MQ=255)
                                   gTGAACTTGGGCCAGCCTCTGTTCAAGCGAGTGGCGCAAGGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgttg                                                   <  1:2374058/93‑1 (MQ=255)
                                   gTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgttg                                                   >  2:896004/1‑93 (MQ=255)
                                   gTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgttg                                                   <  2:3596220/93‑1 (MQ=255)
                                   gTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgttg                                                   >  1:4388283/1‑93 (MQ=255)
                                   gTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgttg                                                   <  1:1502698/93‑1 (MQ=255)
                                     gAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTa                                                 <  1:4493010/93‑1 (MQ=255)
                                      aaCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTAt                                                <  1:2485182/93‑1 (MQ=255)
                                      aaCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTAt                                                >  2:4012356/1‑93 (MQ=255)
                                       aCATGGGCCATCCTCTGTGCAAACAAGGGTCTCAATGGTACGGTTGGTGTGGCATTAATAGCAAGCAGACAGAACCGTTCTGATTGTTGTATg                                               >  2:4657543/1‑93 (MQ=255)
                                          tGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCAt                                            >  2:270826/1‑93 (MQ=255)
                                          tGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCAt                                            >  2:602566/1‑93 (MQ=255)
                                           gggCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATg                                           <  1:4309077/93‑1 (MQ=255)
                                           gggCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATg                                           <  2:1676588/93‑1 (MQ=255)
                                             gCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                         >  2:4551773/1‑93 (MQ=255)
                                             gCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                         >  1:4556742/1‑93 (MQ=255)
                                               cATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                       <  2:265805/93‑1 (MQ=255)
                                               cATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                       >  1:72540/1‑93 (MQ=255)
                                               cATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                       <  1:920226/93‑1 (MQ=255)
                                               cATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                       >  2:4489133/1‑93 (MQ=255)
                                                aTCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttt                                      <  2:1643014/93‑1 (MQ=255)
                                                  ccTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttttt                                    <  2:3063092/93‑1 (MQ=255)
                                                   ctctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttttt                                   >  1:365342/1‑93 (MQ=255)
                                                   ctctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttttt                                   >  2:1039480/1‑93 (MQ=255)
                                                    tctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTa                                  <  2:1752288/93‑1 (MQ=255)
                                                    tctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTa                                  <  1:4492790/93‑1 (MQ=255)
                                                     ctGTGCAAACAAGTGTCTAAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTAt                                 <  2:405791/93‑1 (MQ=255)
                                                      tgtgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATg                                <  2:2138753/93‑1 (MQ=255)
                                                       gtgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGc                               >  2:2192984/1‑93 (MQ=255)
                                                       gtgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGc                               >  1:4328388/1‑93 (MQ=255)
                                                        tgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCt                              <  2:948285/93‑1 (MQ=255)
                                                        tgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCt                              <  2:1447865/93‑1 (MQ=255)
                                                         gCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                             <  1:3888481/93‑1 (MQ=255)
                                                         gCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                             >  2:1975013/1‑93 (MQ=255)
                                                          cAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCttt                            <  1:503931/93‑1 (MQ=255)
                                                            aaCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTcc                          >  1:4613081/1‑93 (MQ=255)
                                                            aaCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTcc                          >  2:2570865/1‑93 (MQ=255)
                                                             aCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCt                         >  2:3890150/1‑93 (MQ=255)
                                                              cAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCtt                        >  1:1182567/1‑93 (MQ=255)
                                                                aGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTaa                      <  2:2863026/93‑1 (MQ=255)
                                                                aGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTaa                      <  2:29368/93‑1 (MQ=255)
                                                                  tgtCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGa                    <  2:1288432/93‑1 (MQ=255)
                                                                   gtCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaa                   >  1:3547868/1‑93 (MQ=255)
                                                                    tctcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaac                  <  1:3324385/93‑1 (MQ=255)
                                                                    tctcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaac                  <  2:1339162/93‑1 (MQ=255)
                                                                     ctcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaaca                 >  2:4430733/1‑93 (MQ=255)
                                                                      tcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaa                >  1:4357384/1‑93 (MQ=255)
                                                                        aaTGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACt              >  2:2758805/1‑93 (MQ=255)
                                                                         aTGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTc             >  2:910891/1‑93 (MQ=255)
                                                                           ggTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAc           >  2:3349279/1‑93 (MQ=255)
                                                                            gTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAcc          <  1:2782888/93‑1 (MQ=255)
                                                                             tACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAccc         >  2:1327281/1‑93 (MQ=255)
                                                                                gTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCtt      >  1:2078332/1‑93 (MQ=255)
                                                                                 tttGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTa     >  2:2503108/1‑93 (MQ=255)
                                                                                  ttGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTaa    <  1:3771927/93‑1 (MQ=255)
                                                                                   tGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTaaa   <  2:2354009/93‑1 (MQ=255)
                                                                                   tGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTaaa   <  1:3089646/93‑1 (MQ=255)
                                                                                    ggTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAg  >  1:1711899/1‑93 (MQ=255)
                                                                                    ggTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAg  >  1:123909/1‑93 (MQ=255)
                                                                                       |                                                                                         
CATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAG  >  NC_000913/2597738‑2597914

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