Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,007,718 2 bp→CG coding (943‑944/1023 nt) pldB → lysophospholipase L(2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,007,7180GC99.0% 266.9 / ‑5.6 102A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (40/61);  minor base A (1/0);  total (41/61)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
*NC_0009134,007,7190CG99.0% 268.9 / NA 102A315R (GCC→CGC) pldBlysophospholipase L(2)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (41/60);  minor base T (0/1);  total (41/61)
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.99e-01

Only 100 of 102 total aligned reads displayed.

TCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTA  >  NC_000913/4007650‑4007787
                                                                    ||                                                                    
tCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGt                                                                >  1:4363265/1‑76 (MQ=255)
 cGAAGGGGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                               <  1:2520997/76‑1 (MQ=255)
 cGAAGGAGGACGGCCGTTGGTAATGAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTc                                                               <  1:2766281/76‑1 (MQ=255)
  gAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcg                                                              <  1:3530069/76‑1 (MQ=255)
  gAAGGAGGACGGCCGTTGGTAAGTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcg                                                              <  2:156443/76‑1 (MQ=255)
   aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcgc                                                             <  1:3461184/76‑1 (MQ=255)
   aaGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTcgc                                                             <  2:1845072/76‑1 (MQ=255)
     ggaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCgcgc                                                           <  2:3057614/76‑1 (MQ=255)
     ggaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCgcgc                                                           <  2:568194/76‑1 (MQ=255)
     ggaggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCgcgc                                                           >  2:3607904/1‑76 (MQ=255)
       aggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTc                                                         <  1:1609460/76‑1 (MQ=255)
        ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTcc                                                        <  1:4535360/76‑1 (MQ=255)
        ggaCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTcc                                                        >  1:3890941/1‑76 (MQ=255)
          aCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                      <  1:3659872/76‑1 (MQ=255)
          aCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                      >  1:3025128/1‑76 (MQ=255)
          aCGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCAc                                                      >  2:3936881/1‑76 (MQ=255)
           cGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACg                                                     <  1:3291673/76‑1 (MQ=255)
               cggtGGTAATTAAAGGGGCTTACCATGAGATCCGTTTTGAAAAGGACGCAATGAGCACAGTCGCGCTCCACGCTAt                                                 >  1:3933008/4‑76 (MQ=255)
               cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                 <  2:2818396/76‑1 (MQ=255)
               cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                 <  1:3612657/76‑1 (MQ=255)
               cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                 <  2:3744663/76‑1 (MQ=255)
               cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                 <  2:1599001/76‑1 (MQ=255)
               cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                 >  1:1933468/1‑76 (MQ=255)
               cGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCAt                                                 <  2:2283959/76‑1 (MQ=255)
                ggtGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATc                                                <  2:1208621/74‑1 (MQ=255)
                gTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATc                                                >  2:3508267/1‑76 (MQ=255)
                 ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCg                                               >  1:1168995/1‑76 (MQ=255)
                 ttGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCg                                               >  1:1150953/1‑76 (MQ=255)
                       aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCCTCGTTGAtt                                         >  2:2049261/1‑76 (MQ=255)
                       aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtt                                         >  1:2236013/1‑76 (MQ=255)
                       aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtt                                         >  1:2270285/1‑76 (MQ=255)
                       aTTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGCCGCGCTCCACGCCCTCGTTGAtt                                         >  2:1497409/1‑76 (MQ=255)
                        ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCCTCGTTGAttt                                        >  2:2866392/1‑76 (MQ=255)
                        ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttt                                        >  2:4058039/1‑76 (MQ=255)
                        ttAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttt                                        <  1:1682132/76‑1 (MQ=255)
                         tAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttt                                       <  1:360510/76‑1 (MQ=255)
                          aaaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAttttt                                      <  2:1401793/76‑1 (MQ=255)
                           aaGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAtttttt                                     >  2:4545853/1‑76 (MQ=255)
                            aGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTc                                    <  2:4865463/76‑1 (MQ=255)
                            aGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTc                                    <  2:4720992/76‑1 (MQ=255)
                             ggTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCa                                   >  2:3677176/1‑76 (MQ=255)
                              gTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCaa                                  <  2:1797469/76‑1 (MQ=255)
                              gTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCaa                                  >  1:3205242/1‑76 (MQ=255)
                                 cTTACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAg                               <  2:2563708/76‑1 (MQ=255)
                                  ttACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAgg                              >  1:4115760/1‑76 (MQ=255)
                                   tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                             <  1:2374414/76‑1 (MQ=255)
                                   tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                             >  1:1622718/1‑76 (MQ=255)
                                   tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                             <  2:3067102/76‑1 (MQ=255)
                                   tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                             >  1:1577243/1‑76 (MQ=255)
                                   tACCATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGc                             <  2:225350/76‑1 (MQ=255)
                                     ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                           <  2:2416138/76‑1 (MQ=255)
                                     ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                           <  2:460807/76‑1 (MQ=255)
                                     ccATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCAt                           >  2:849497/1‑76 (MQ=255)
                                      cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATa                          <  1:1391552/76‑1 (MQ=255)
                                      cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATa                          <  1:4058039/76‑1 (MQ=255)
                                      cATGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATa                          <  2:4356610/76‑1 (MQ=255)
                                        tGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCACCAGGCATAAc                        <  1:2472706/76‑1 (MQ=255)
                                        tGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                        >  2:1659261/1‑76 (MQ=255)
                                        tGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                        >  1:4069539/1‑76 (MQ=255)
                                        tGAGATCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAAc                        <  2:425241/76‑1 (MQ=255)
                                        tGAGAGCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTCGATTTTTTCAACAGGCATAAc                        >  2:4627371/1‑76 (MQ=255)
                                         gagaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACt                       <  1:266564/76‑1 (MQ=255)
                                          agaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTc                      >  2:585249/1‑76 (MQ=255)
                                           gaTCCTTTTTGAATAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACCGGCATAACTCa                     <  1:1394740/76‑1 (MQ=255)
                                           gaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGCTTTTTTCAACAGGCCTCACTcc                     >  2:1100767/1‑75 (MQ=255)
                                           gaTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCa                     >  2:748732/1‑76 (MQ=255)
                                            ctCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCTTTGATTTTTTCAACAGGCATAACTCAc                    <  1:3021220/75‑1 (MQ=255)
                                            aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                    >  2:567566/1‑76 (MQ=255)
                                            aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                    <  2:3844141/76‑1 (MQ=255)
                                            aTCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACACGCATAACTCAc                    >  2:3584855/1‑76 (MQ=255)
                                            aTCCGTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAc                    <  2:4331512/76‑1 (MQ=255)
                                             tCCTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAcc                   <  2:3205242/76‑1 (MQ=255)
                                              ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                  <  2:2270284/76‑1 (MQ=255)
                                              ccTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAccc                  <  2:2236012/76‑1 (MQ=255)
                                               cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                 >  1:2124131/1‑76 (MQ=255)
                                               cTTTTTGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCa                 <  2:1816834/76‑1 (MQ=255)
                                               cTTTTTGAAAAGGACGCAATGCGCTCAATCGCGCTACACGCCATCGTTGATGTTTTAAACAGGCTGAACTcacaca                 >  2:1301778/1‑76 (MQ=255)
                                                tttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAg                <  2:4294184/76‑1 (MQ=255)
                                                 ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc               >  1:814327/1‑76 (MQ=255)
                                                 ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc               >  2:2778339/1‑76 (MQ=255)
                                                 ttttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGc               <  2:2574681/76‑1 (MQ=255)
                                                 ttttGAAAAGGACGCAACGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTATCAACAGGCATAACTCACCCAGc               <  2:693774/76‑1 (MQ=255)
                                                  tttGAAAAGGACGCAACGCGCTCAGTACCGCGCCAAGCCATCGTTGATTGTTTCAACAGGCAGAACGCACCCAgag              >  1:1678135/1‑74 (MQ=255)
                                                   ttGAAAAGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCAcc                   >  1:1580666/1‑70 (MQ=255)
                                                        aatgACGCAATGCTCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc        <  2:3253458/73‑1 (MQ=255)
                                                        aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc        <  2:2790252/76‑1 (MQ=255)
                                                        aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc        <  2:445003/76‑1 (MQ=255)
                                                        aaGGACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGAGTTTTTCAACAGGCATAACTCACCCAGCGGAAAcc        <  2:3429358/76‑1 (MQ=255)
                                                          ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGc      <  2:4777539/76‑1 (MQ=255)
                                                          ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGc      >  2:1578168/1‑76 (MQ=255)
                                                          ggACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGc      <  2:3615649/76‑1 (MQ=255)
                                                          ggACGCAATGCGCTCAGTCGCGCTCCAAGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCga      >  1:957375/1‑75 (MQ=255)
                                                           gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGct     <  2:1296186/76‑1 (MQ=255)
                                                           gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGct     >  1:3793415/1‑76 (MQ=255)
                                                           gACGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAAGCATAACTCACCCAGCGGAAACCGct     <  2:1160132/76‑1 (MQ=255)
                                                            gcGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc    <  2:268782/75‑1 (MQ=255)
                                                            aCGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc    >  1:2620120/1‑76 (MQ=255)
                                                            aCGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctc    <  1:3254332/76‑1 (MQ=255)
                                                             cGCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGctct   <  2:2363286/76‑1 (MQ=255)
                                                              gCAATGCGCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTa  <  1:3221961/76‑1 (MQ=255)
                                                                    ||                                                                    
TCGAAGGAGGACGGCCGTTGGTAATTAAAGGTGCTTACCATGAGATCCTTTTTGAAAAGGACGCAATGGCCTCAGTCGCGCTCCACGCCATCGTTGATTTTTTCAACAGGCATAACTCACCCAGCGGAAACCGCTCTA  >  NC_000913/4007650‑4007787

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