Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,419,010 (C)6→8 coding (370/648 nt) yfcG → putative S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,419,0041.C100.0% 321.8 / NA 105A122P (GCC→CCC) yfcGputative S‑transferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (40/65);  total (40/65)
*NC_0009132,419,0042.C99.1% 312.3 / ‑3.6 105coding (364/648 nt)yfcGputative S‑transferase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (40/64);  total (40/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

Only 100 of 117 total aligned reads displayed.

GCGTGAGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑CC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAACAAGCGGCTGG  >  NC_000913/2418913‑2419085
                                                                                            ||                                                                                   
gcgTGAGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑c                                                                                    <  2:3376762/93‑1 (MQ=255)
  gTGAGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑cc‑‑c                                                                                >  2:2438053/1‑93 (MQ=255)
     aGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑cc‑‑cccc                                                                             >  2:3197388/1‑93 (MQ=255)
      gcgcgcCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑cc‑‑ccccc                                                                            >  1:55782/1‑92 (MQ=255)
        gcgcCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGAGTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑cc‑‑cccccca                                                                          <  2:4315217/93‑4 (MQ=255)
         cgccgcCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑cc‑‑ccccccaa                                                                         <  1:813121/93‑5 (MQ=255)
                acaTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTc                                                                    <  1:1039731/93‑1 (MQ=255)
                 caTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTcc                                                                   <  1:2164423/93‑1 (MQ=255)
                 caTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTcc                                                                   <  1:1854827/93‑1 (MQ=255)
                 caTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTcc                                                                   >  2:4522312/1‑93 (MQ=255)
                 caTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTcc                                                                   <  2:563413/93‑1 (MQ=255)
                  aTTACAGTGTTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGACAAAAGCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCt                                                                  <  1:1115365/93‑1 (MQ=255)
                  aTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCt                                                                  >  2:4463294/1‑93 (MQ=255)
                  aTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCt                                                                  <  2:1253090/93‑1 (MQ=255)
                  aTTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCt                                                                  <  2:3905493/93‑1 (MQ=255)
                   ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCtt                                                                 <  2:3683678/93‑1 (MQ=255)
                   ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCtt                                                                 >  1:2205084/1‑93 (MQ=255)
                   ttACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCtt                                                                 <  2:711109/93‑1 (MQ=255)
                    tccggtgtTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAACCCCGCAGCCCC‑‑CCCCAAACCATTCCTTa                                                                <  2:6673/89‑1 (MQ=255)
                    tACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTa                                                                <  2:274345/93‑1 (MQ=255)
                     aCAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTAc                                                               <  2:1698426/93‑1 (MQ=255)
                     aCAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTAc                                                               <  1:825229/93‑1 (MQ=255)
                       aGTGTTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGc                                                             <  2:4421127/93‑1 (MQ=255)
                       aGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGc                                                             <  1:1269561/93‑1 (MQ=255)
                        gTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCt                                                            >  1:2836233/1‑93 (MQ=255)
                         tGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTa                                                           <  1:731398/93‑1 (MQ=255)
                          ggTTATTCTGGCAGGTAGGTGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGAAG‑‑CCCACCCCCAACCATTCCGTACGCTAt                                                          >  2:784428/1‑93 (MQ=255)
                          ggTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTAt                                                          <  1:1199800/93‑1 (MQ=255)
                           gTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTAtt                                                         <  2:4483893/93‑1 (MQ=255)
                            ttATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTg                                                        >  1:588065/1‑93 (MQ=255)
                             tATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGa                                                       >  2:3333459/1‑93 (MQ=255)
                             tATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGa                                                       <  1:3506125/93‑1 (MQ=255)
                             tATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGa                                                       >  2:4146318/1‑93 (MQ=255)
                              aTTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGaa                                                      <  1:143689/93‑1 (MQ=255)
                               ttCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAAc                                                     >  2:411970/1‑93 (MQ=255)
                               ttCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAAc                                                     <  2:3541875/93‑1 (MQ=255)
                               ttCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAAc                                                     >  1:2722312/1‑93 (MQ=255)
                                   ggCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTa                                                 >  1:4068769/1‑93 (MQ=255)
                                    gCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTAt                                                <  1:3123244/93‑1 (MQ=255)
                                     cAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATc                                               <  2:3079418/93‑1 (MQ=255)
                                      aGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCa                                              >  2:3208697/1‑93 (MQ=255)
                                       ggTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAg                                             <  2:223143/93‑1 (MQ=255)
                                       ggTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAg                                             >  1:1810997/1‑93 (MQ=255)
                                       ggTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAg                                             <  2:1664493/93‑1 (MQ=255)
                                       ggTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAg                                             >  1:987716/1‑93 (MQ=255)
                                       ggTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAg                                             <  1:4488522/93‑1 (MQ=255)
                                         tAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGt                                           <  1:1223269/93‑1 (MQ=255)
                                          aGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGtt                                          <  1:2171276/93‑1 (MQ=255)
                                              ggACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGaaa                                      >  1:3927381/1‑93 (MQ=255)
                                               gACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAAc                                     >  2:2587183/1‑93 (MQ=255)
                                               gACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAAc                                     >  1:1451368/1‑93 (MQ=255)
                                                gcTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACt                                    <  1:2789869/92‑1 (MQ=255)
                                                aCTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACt                                    >  2:3538417/1‑93 (MQ=255)
                                                  tGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCa                                  <  1:698791/93‑1 (MQ=255)
                                                   ggggCCTATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAg                                 <  2:1036893/93‑1 (MQ=255)
                                                   ggggCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAg                                 <  1:2857357/93‑1 (MQ=255)
                                                    gggCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGc                                <  2:4433304/93‑1 (MQ=255)
                                                      gCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGt                              <  1:3022124/93‑1 (MQ=255)
                                                      gCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGt                              >  2:1870391/1‑93 (MQ=255)
                                                       ccGATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTc                             <  2:1094181/93‑1 (MQ=255)
                                                         gATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCtt                           <  2:1690196/93‑1 (MQ=255)
                                                         gATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCtt                           <  1:4582741/93‑1 (MQ=255)
                                                         gATGCTTGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCtt                           <  1:2633485/93‑1 (MQ=255)
                                                          aTGCTTGGGAAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACAATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCttt                          <  2:1707596/93‑1 (MQ=255)
                                                              ttGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCa                      >  2:1281388/1‑93 (MQ=255)
                                                               tGGGCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCAt                     <  1:3993136/93‑1 (MQ=255)
                                                                gggCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATg                    <  2:2755324/93‑1 (MQ=255)
                                                                 ggCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGt                   <  2:1679745/93‑1 (MQ=255)
                                                                  gCAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTa                  <  1:3629936/93‑1 (MQ=255)
                                                                   cAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTAc                 <  1:1301598/93‑1 (MQ=255)
                                                                   cAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTAc                 <  1:3836102/93‑1 (MQ=255)
                                                                   cAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTAc                 <  1:1777661/93‑1 (MQ=255)
                                                                   cAAAATCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTAc                 <  2:3920032/93‑1 (MQ=255)
                                                                    aaaaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACt                >  1:3158896/1‑93 (MQ=255)
                                                                     aaaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTg               <  2:28429/93‑1 (MQ=255)
                                                                     aaaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTg               >  2:2903041/1‑93 (MQ=255)
                                                                     aaaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTg               >  2:4010241/1‑93 (MQ=255)
                                                                     aaaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTg               >  2:229873/1‑93 (MQ=255)
                                                                     aaaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTg               >  1:4196043/1‑93 (MQ=255)
                                                                      aaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGa              <  2:2563155/93‑1 (MQ=255)
                                                                      aaTCATCATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGa              <  1:491126/93‑1 (MQ=255)
                                                                       atcatcATTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGag             >  1:778039/1‑92 (MQ=255)
                                                                        tcatcaTTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTTCGCTATTGAACGTTATCAGGGTGAAACTCAGCGTCTTTACCATGTACTGagc            >  2:728210/1‑91 (MQ=255)
                                                                        tcatcaTTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGagc            >  2:3654943/1‑91 (MQ=255)
                                                                        tcatcaTTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGagc            >  2:1986530/1‑91 (MQ=255)
                                                                          atcatTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCaa          <  1:2782520/93‑1 (MQ=255)
                                                                          atcatTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCaa          <  2:2511102/93‑1 (MQ=255)
                                                                          atcatTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCaa          <  2:1689971/93‑1 (MQ=255)
                                                                           tcatTTTAATCACGCCGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAg         <  2:1822157/93‑1 (MQ=255)
                                                                           tcatTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAg         >  1:2313634/1‑93 (MQ=255)
                                                                           tcatTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAg         >  1:4654891/1‑93 (MQ=255)
                                                                            catTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGc        <  1:723773/93‑1 (MQ=255)
                                                                             atTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCg       >  1:1717984/1‑93 (MQ=255)
                                                                             atTTTAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCg       <  1:2835805/93‑1 (MQ=255)
                                                                              ttttAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCgg      >  2:4508250/1‑93 (MQ=255)
                                                                              ttttAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCACGTTGACACTCAGCGTCTTTACCATGTACTGAGCAAGCgg      >  2:2025812/1‑93 (MQ=255)
                                                                                ttAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCt    >  1:3883191/1‑93 (MQ=255)
                                                                                ttAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCt    >  1:2138226/1‑93 (MQ=255)
                                                                                 tAATCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTg   <  2:3564497/93‑1 (MQ=255)
                                                                                  aaTCACGCAGCCCC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTgg  <  1:4152770/93‑1 (MQ=255)
                                                                                            ||                                                                                   
GCGTGAGCGCGCCGCCACATTACAGTGGTTATTCTGGCAGGTAGGCGGACTGGGGCCGATGCTTGGGCAAAATCATCATTTTAATCACGCAG‑‑CC‑‑CCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAACAAGCGGCTGG  >  NC_000913/2418913‑2419085

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