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.6% 823.1 / ‑4.8 252intergenic (‑45/‑101)dapA/b2479dihydrodipicolinate synthase/GcvR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (119/132);  total (119/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.76e-01

Only 100 of 260 total aligned reads displayed.

TTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAG  >  NC_000913/2597735‑2597914
                                                                                          |                                                                                         
tttCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGtgtg                                                                                         >  2:5929606/1‑93 (MQ=255)
tttCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGtgtg                                                                                         >  1:703632/1‑93 (MQ=255)
tttCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGtgtg                                                                                         >  1:4730056/1‑93 (MQ=255)
  tCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGc                                                                                       >  1:4037252/1‑93 (MQ=255)
   catccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCa                                                                                      <  2:6182252/93‑1 (MQ=255)
   catccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCa                                                                                      >  2:3303906/1‑93 (MQ=255)
    atccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAt                                                                                     >  2:5987245/1‑93 (MQ=255)
    atccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAt                                                                                     >  2:5459031/1‑93 (MQ=255)
    atccatcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCAt                                                                                     >  2:2844945/1‑93 (MQ=255)
       catcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaa                                                                                  <  2:4751940/93‑1 (MQ=255)
       catcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaa                                                                                  >  1:6201227/1‑93 (MQ=255)
        atcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaa                                                                                 >  1:6120503/1‑93 (MQ=255)
         tcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaaa                                                                                >  1:109882/1‑93 (MQ=255)
         tcGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTaaaa                                                                                <  1:2360538/93‑1 (MQ=255)
            gAGTAACAATCGCAACAATACTGCCCGTGAACATGAACCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagca                                                                             <  1:5308361/93‑1 (MQ=255)
              gTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaag                                                                           <  1:5818184/93‑1 (MQ=255)
               tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaagc                                                                          <  2:510651/93‑1 (MQ=255)
               tAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAaagcaagc                                                                          <  2:1207400/93‑1 (MQ=255)
                aaCAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAagcaagca                                                                         <  1:2252852/93‑1 (MQ=255)
                aaCAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAagcaagca                                                                         <  1:1165854/93‑1 (MQ=255)
                aaCAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAagcaagca                                                                         <  1:2019297/93‑1 (MQ=255)
                 aCAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcag                                                                        >  2:5169258/1‑93 (MQ=255)
                  cAATCGCGACAATACTTCCAGTGAACATGGGCCATCCTCTGTGCAAACAAGTGGCTCAATGGTACGTTTGGTGAGGCATTAAaagcaagccgc                                                                       >  1:2854314/1‑90 (MQ=255)
                    aTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagaca                                                                     >  1:1186140/1‑93 (MQ=255)
                    aTCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagaca                                                                     >  1:1316694/1‑93 (MQ=255)
                      cgcgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagacaga                                                                   <  1:4191479/93‑1 (MQ=255)
                      cgcgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGcagacaga                                                                   >  2:1688774/1‑93 (MQ=255)
                       gcgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGaa                                                                  <  2:2797970/93‑1 (MQ=255)
                        cgACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAc                                                                 <  1:466157/93‑1 (MQ=255)
                          aCAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAAcct                                                               >  2:1771477/1‑92 (MQ=255)
                          aCAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCg                                                               <  2:463879/93‑1 (MQ=255)
                          aCAAGACTGCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCg                                                               <  1:6304950/93‑1 (MQ=255)
                           cAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGt                                                              <  1:5646084/93‑1 (MQ=255)
                           aaaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGt                                                              <  2:3986123/92‑1 (MQ=255)
                            aaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGtt                                                             <  2:4992050/93‑1 (MQ=255)
                            aaTACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGtt                                                             >  2:5840001/1‑93 (MQ=255)
                              tACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCt                                                           >  1:743216/1‑93 (MQ=255)
                              tACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCt                                                           >  1:1699793/1‑93 (MQ=255)
                              tACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCt                                                           <  2:3131913/93‑1 (MQ=255)
                               aCTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTg                                                          <  1:5960901/93‑1 (MQ=255)
                               aCTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTg                                                          <  1:805420/93‑1 (MQ=255)
                                 ttCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAt                                                        <  1:5461394/93‑1 (MQ=255)
                                 ttCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAt                                                        >  2:2560046/1‑93 (MQ=255)
                                  tCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAtt                                                       >  1:3738289/1‑93 (MQ=255)
                                   cccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttg                                                      <  2:6604173/93‑1 (MQ=255)
                                    ccGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgt                                                     >  1:3802900/1‑93 (MQ=255)
                                     cGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgtt                                                    >  2:6313857/1‑93 (MQ=255)
                                      gTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGAttgttg                                                   <  1:3645578/93‑1 (MQ=255)
                                       tGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATtgttgt                                                  >  2:3921888/1‑93 (MQ=255)
                                        gAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTa                                                 >  2:1399072/1‑93 (MQ=255)
                                        gAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTa                                                 <  2:4073369/93‑1 (MQ=255)
                                        gAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTa                                                 <  1:1230239/93‑1 (MQ=255)
                                        gAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTa                                                 <  2:919115/93‑1 (MQ=255)
                                          aCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATg                                               <  2:5712431/93‑1 (MQ=255)
                                          aCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATg                                               >  2:41899/1‑93 (MQ=255)
                                          aCATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATg                                               >  1:4353691/1‑93 (MQ=255)
                                          aCATGGGCCATCCTCTGTGCAAACAAGAGGCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACGGCTCTCCTTGTTGTATg                                               >  1:1030377/1‑93 (MQ=255)
                                             tGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCAt                                            >  1:6653146/1‑93 (MQ=255)
                                              gggCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATg                                           <  1:2117669/93‑1 (MQ=255)
                                              gggCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATg                                           <  2:3820430/93‑1 (MQ=255)
                                               ggCCATCCGCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGt                                          <  1:3720816/93‑1 (MQ=255)
                                                gCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                         >  1:1895485/1‑93 (MQ=255)
                                                gCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                         <  1:3363449/93‑1 (MQ=255)
                                                gCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                         <  2:2466154/93‑1 (MQ=255)
                                                  cATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                       <  1:4241227/93‑1 (MQ=255)
                                                  cATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                       >  1:4586900/1‑93 (MQ=255)
                                                   aTCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttt                                      >  2:3643029/1‑93 (MQ=255)
                                                    tCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttt                                     >  2:4474449/1‑93 (MQ=255)
                                                    tCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttt                                     >  2:726888/1‑93 (MQ=255)
                                                    tCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttt                                     >  2:1609484/1‑93 (MQ=255)
                                                    tCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttt                                     >  1:6543420/1‑93 (MQ=255)
                                                     ccTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttttttt                                    <  2:3753582/93‑1 (MQ=255)
                                                      ctctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttttt                                   >  1:6410229/1‑93 (MQ=255)
                                                       tctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTa                                  <  2:5610516/93‑1 (MQ=255)
                                                       tctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTa                                  <  2:4962357/93‑1 (MQ=255)
                                                        ctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTAt                                 <  1:6034990/93‑1 (MQ=255)
                                                        ctGTGCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTAt                                 <  2:192301/93‑1 (MQ=255)
                                                           tgCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCt                              <  1:2332022/93‑1 (MQ=255)
                                                            tcAACCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGCTGTATGCATGTTTTTTTTATGCtt                             <  1:3263530/92‑1 (MQ=255)
                                                            gCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                             <  1:1575072/93‑1 (MQ=255)
                                                            gCAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCtt                             >  1:4524039/1‑93 (MQ=255)
                                                             cAAACAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTGTATGCttt                            <  1:6528428/93‑1 (MQ=255)
                                                                aCAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCt                         <  2:2954685/93‑1 (MQ=255)
                                                                 cAAGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCtt                        >  2:2216070/1‑93 (MQ=255)
                                                                  aaGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTa                       <  1:2475093/93‑1 (MQ=255)
                                                                   aGTGTCTCAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTaa                      <  2:4162289/93‑1 (MQ=255)
                                                                         tcAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaa                >  1:3602084/1‑93 (MQ=255)
                                                                          cAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaac               <  2:1379087/93‑1 (MQ=255)
                                                                          cAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaac               <  1:5161553/93‑1 (MQ=255)
                                                                          cAATGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGaacaac               <  2:3530760/93‑1 (MQ=255)
                                                                             tGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCa            >  1:927668/1‑93 (MQ=255)
                                                                             tGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCa            <  2:2257796/93‑1 (MQ=255)
                                                                             tGGTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCa            <  1:4934023/93‑1 (MQ=255)
                                                                               gTACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAcc          >  1:1761856/1‑93 (MQ=255)
                                                                                tACGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCAccc         >  2:340640/1‑93 (MQ=255)
                                                                                  cGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCt       <  2:1966287/93‑1 (MQ=255)
                                                                                  cGTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCt       >  1:406821/1‑93 (MQ=255)
                                                                                   gTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCtt      <  1:5213878/93‑1 (MQ=255)
                                                                                   gTTTGGTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCtt      >  2:3496192/1‑93 (MQ=255)
                                                                                       ggTGTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAg  <  1:3094831/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTTCATCCATCGGAGTAACAATCGCGACAATACTTCCCGTGAACATGGGCCATCCTCTGTGCAAACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAG  >  NC_000913/2597735‑2597914

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