Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,430,453 +C coding (745/747 nt) smf ← putative Rossmann‑fold nucleotide‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,430,4531.C98.2% 295.5 / ‑1.0 108coding (745/747 nt)smfputative Rossmann‑fold nucleotide‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (61/45);  total (62/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTC  >  NC_000913/3430378‑3430522
                                                                             |                                                                     
gTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAgg                                                                        >  1:5283571/1‑76 (MQ=255)
gTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAgg                                                                        <  1:3994824/76‑1 (MQ=255)
gTCACAAGAATAGCTCCTTGTTTTATTAACCAGTGAGGGCCTTCGCTTCGCGGATTCCCTATTGGACCTGGCAAgg                                                                        <  2:2103329/76‑1 (MQ=255)
 tCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAggc                                                                       >  2:4281977/1‑75 (MQ=255)
 tCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCGATTGGACCTGGCAAggc                                                                       >  2:681476/1‑75 (MQ=255)
  cacaAGAATCGCGCCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAggca                                                                      <  2:2278585/76‑3 (MQ=255)
    caAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGcaaa                                                                   >  2:1210039/1‑76 (MQ=255)
     aaGAATCGCACCTTGTTTTAGTAACCAGTGAGGCCCTTCGCTTCCCGGATTACCTATTGGACCCGGCAATGcaaaa                                                                  >  2:517608/1‑76 (MQ=255)
       gAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAAc                                                                 >  2:4514614/1‑76 (MQ=255)
       gAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAAc                                                                 <  1:3889106/76‑1 (MQ=255)
        aaTCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACt                                                                <  2:917543/76‑1 (MQ=255)
        aaTCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACt                                                                >  2:4236772/1‑76 (MQ=255)
         aTCGCACCTTGTTTTATTAACCATTGAGGCCCTTCGCTTCGCTGATTCCCTATTGGAGCTGGCAAGGCAACACttt                                                              <  1:4729393/76‑2 (MQ=255)
         aTCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACtt                                                               <  2:3456204/76‑1 (MQ=255)
           cGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtct                                                             <  1:425952/76‑1 (MQ=255)
            gCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                            >  1:3856346/1‑76 (MQ=255)
            gCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                            <  2:560946/76‑1 (MQ=255)
            gCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                            >  1:3256546/1‑76 (MQ=255)
             cgccTTGTTTTATTAGCCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGATCTGGCAAGGCAAATACTTCTCg                                                           <  2:3649188/74‑1 (MQ=255)
             cACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCg                                                           >  1:3766467/1‑76 (MQ=255)
              cccTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGc                                                          <  2:908286/75‑1 (MQ=255)
              aCCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGc                                                          <  2:5091390/76‑1 (MQ=255)
               ccTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGcc                                                         <  2:1052032/76‑1 (MQ=255)
               ccTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGcc                                                         <  1:2730185/76‑1 (MQ=255)
                  tGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAACCTTCTCGCCCCt                                                      <  1:1269183/76‑1 (MQ=255)
                  tGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCt                                                      <  2:2973616/76‑1 (MQ=255)
                   gTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTg                                                     <  2:5647958/76‑1 (MQ=255)
                   gTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTg                                                     <  2:2709633/76‑1 (MQ=255)
                    ttttATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCCGc                                                    >  2:5215622/1‑76 (MQ=255)
                     tttATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCt                                                   >  2:5396549/1‑76 (MQ=255)
                      ttattaACCAGTGAGGCCCTTCTCTTCACGGATTCCCTCTTGGACCCGGCAAGCCCCACCCTTCTCGCCCCTGcac                                                  >  2:1242195/1‑74 (MQ=255)
                      ttattaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTc                                                  >  2:1230767/1‑76 (MQ=255)
                       tattaACCAGTGAGGCCCTTCGCTTCCAGGCTTCCATATTGGACCTGGTAAGGCAAAAACGTCTCGCACCAGCtaa                                                 >  1:2883873/1‑74 (MQ=255)
                        attaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAACAACTTCTCGCCCCTGCTcca                                                >  2:5317414/1‑74 (MQ=255)
                        attaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCaa                                                <  2:2175693/76‑1 (MQ=255)
                        attaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCaa                                                >  2:4266701/1‑76 (MQ=255)
                          taACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCACgc                                              >  2:3031537/1‑76 (MQ=255)
                          taACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAgc                                              <  1:5463557/76‑1 (MQ=255)
                          taACCAGTGAGGCCCTTCGCGTCACGGATTCCCTATGGGACCTGGCAAGGCAAAAACTTCTCGCCCCGGCTCAAgc                                              <  2:4524192/76‑1 (MQ=255)
                           aaCCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAgcg                                             >  2:4078843/1‑76 (MQ=255)
                            aCCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCACgcgc                                            >  2:2216896/1‑76 (MQ=255)
                            aCCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGATCAAgcgc                                            >  1:853319/1‑76 (MQ=255)
                             ccAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGca                                           >  2:2020636/1‑76 (MQ=255)
                             ccAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGca                                           <  2:2347321/76‑1 (MQ=255)
                               aGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGcaca                                         >  2:2142279/1‑76 (MQ=255)
                               aGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGcaca                                         >  1:2758582/1‑76 (MQ=255)
                                 tGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCTCCCTCGCTCAAGCGCACAAc                                       >  2:1823785/1‑76 (MQ=255)
                                 tGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAAc                                       <  1:2319114/76‑1 (MQ=255)
                                 tGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAAc                                       >  2:2400059/1‑76 (MQ=255)
                                  gAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACg                                      <  1:4578239/76‑1 (MQ=255)
                                  gAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACg                                      >  1:3165291/1‑76 (MQ=255)
                                  gAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACg                                      <  1:2598797/76‑1 (MQ=255)
                                   aGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGt                                     >  1:551265/1‑76 (MQ=255)
                                      cccTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCt                                  >  1:5726328/1‑76 (MQ=255)
                                      cccTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGAGCACAACGTGCt                                  >  2:4288242/1‑76 (MQ=255)
                                       ccTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTAGCCCCTGCTCAAGCGCACAACGTGCTg                                 <  1:3079490/76‑1 (MQ=255)
                                        cTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGc                                   <  2:4458867/73‑1 (MQ=255)
                                        cTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGc                                   >  1:4458868/1‑73 (MQ=255)
                                         ttCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTc                               <  1:846505/76‑1 (MQ=255)
                                         ttCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGCGCTGTc                               >  2:1481093/1‑76 (MQ=255)
                                         ttCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACCTGCTGTc                               >  2:5516756/1‑76 (MQ=255)
                                           cGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCAc                             >  2:799640/1‑76 (MQ=255)
                                           cGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCAc                             <  1:5278545/76‑1 (MQ=255)
                                            gCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCAcc                            >  1:122438/1‑76 (MQ=255)
                                            gCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCAcc                            >  1:3900509/1‑76 (MQ=255)
                                              ttCCCGGCTTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAg                          >  2:1656076/1‑76 (MQ=255)
                                              ttCCCGGATTCCCTCTTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAg                          >  2:3940469/1‑76 (MQ=255)
                                                  cGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGac                      >  2:793575/1‑75 (MQ=255)
                                                  cGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGaa                      >  1:4341746/1‑76 (MQ=255)
                                                  cGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGaa                      >  1:2183957/1‑76 (MQ=255)
                                                   ggATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCa                           >  1:5225836/1‑70 (MQ=255)
                                                   ggATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCa                           <  2:5225835/70‑1 (MQ=255)
                                                   ggATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAAc                     <  2:3091157/76‑1 (MQ=255)
                                                    gATTCCCTATTGGACCTGGCAAGGCAAGAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAGcc                    <  2:360912/76‑1 (MQ=255)
                                                    gATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAAcc                    >  1:2810762/1‑76 (MQ=255)
                                                       tCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCAcc                 >  2:687624/1‑75 (MQ=255)
                                                       tCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACt                 <  1:1757330/76‑1 (MQ=255)
                                                       tCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACt                 >  2:4444494/1‑76 (MQ=255)
                                                        cccTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTa                >  1:1171006/1‑76 (MQ=255)
                                                         ccTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTAc               >  1:4226244/1‑76 (MQ=255)
                                                          cTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACg              <  1:4569290/76‑1 (MQ=255)
                                                           tATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCACGCGCACAACGTGCTGTCACCAGCGAACCACTACGc             >  2:5299861/1‑76 (MQ=255)
                                                           tATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTAccc             >  2:386862/1‑74 (MQ=255)
                                                            cttGGACCAGGCAAGTCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCa            <  1:1210039/75‑1 (MQ=255)
                                                            aTTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCa            >  1:4830113/1‑76 (MQ=255)
                                                            aTTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCa            >  2:1894114/1‑76 (MQ=255)
                                                              tGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCaaa          <  1:505750/76‑1 (MQ=255)
                                                               ggACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAg         >  1:3037120/1‑76 (MQ=255)
                                                               ggACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAg         >  1:1756509/1‑76 (MQ=255)
                                                               ggACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAg         <  1:1238257/76‑1 (MQ=255)
                                                               ggACCTGGCAAGGCAAAAACTTCTCGCCCATGCTCAAGCGCATAACGTGCTGACACCTGCGAACAACGACGTAAAg         >  2:1823709/1‑76 (MQ=255)
                                                                 aCCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGcc       <  1:4236773/76‑1 (MQ=255)
                                                                 aCCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGcc       <  1:4720668/76‑1 (MQ=255)
                                                                  ccTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCATCAGCGAAACACTCCGCAAAGCCg      >  1:4695507/1‑76 (MQ=255)
                                                                  ccTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCg      >  2:3494050/1‑76 (MQ=255)
                                                                   cTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGc     <  1:4651995/76‑1 (MQ=255)
                                                                    tGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGNGAACCACTACGCAAAGCCGCt    <  1:4078844/76‑1 (MQ=255)
                                                                    tGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCt    >  1:2233279/1‑76 (MQ=255)
                                                                      gCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTc  >  2:3981241/1‑76 (MQ=255)
                                                                      gCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTc  >  2:3775286/1‑76 (MQ=255)
                                                                             |                                                                     
GTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTC  >  NC_000913/3430378‑3430522

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